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CHAPTER.  1 
INTRODUCTION 

In  high  rise  buildings,  it  is  important  to  provide  the 
lateral  resistance  of  the  frames  to  the  loads  due  to  winds, 
earthquakes  or  blasts.  These  forces  can  develop  high 
stresses,  and  yield  lateral  movements  ,  thereby  causing 
discomfort  to  the  occupants.  In  addition  to  stiffening  the 
frames  due  to  lateral  loads,  the  shear  walls  support  vertical 
loads.  Although  reinforced  concrete  walls  were  often  used  to 
provide  the  necessary  resistance  to  horizontal  forces,  modern 
concept  of  shear  walls  designed  as  vertical  slender 
cantilevers  was  first  used  in  1948  in  housing  projects  in  New 
York  City  and  in  Chicago  in  buildings  designed  for  wind 
forces . 

In  nonearthquake  areas  reinforced  concrete  buildings 
with  more  than  IS  to  20  stories  are  usually  designed  with 
shear  walls  to  improve  their  stiffness.  Economically,  the 
inclusion  of  shear  walls  seems  to  be  the  least  expensive  way 
to  increase  the  overall  rigidity  of  building  structures. 
Shear  walls  serve  the  functions  as  follow:  a)  to  support 
gravity  loads  b)  to  provide  lateral  resistance  ,  and  c)  to 
function  as  a  wall  . 

Usually  most  of  the  recent  high  rise  reinforced  concrete 
buildings  were  built  without  any  additional  cost  for  making 
the  design  stresses  in  the  range  of  allowable  stresses  when 
lateral   loads   are  considered  by  the  use  of  the  shear   walls. 


The   high   lateral  rigidity  was  also  achieved  as  one   of   the 

results  of  the  shear  wall-frame  interaction. 

OBJECTIVES 

The  merits  of  shear  walls  for  earthquake  resistance,  and 
examinations  of  their  performance  in  earthquakes  were  made  by 
many  researchers. 

In  this  study,  the  main  objectives  are    as  follows: 

a)  to  suggest  the  use  of  the  finite  element  method,  a 
convenient  numerical  analysis,  by  which  an  accurate  estimate 
may  be  made. 

b)  to  obtain  experience  in  developing  programs  using 
the  finite  element  method. 

c)  to  learn  the  stress  distributions  of  the  shear  walls 
under  the  application  of  design  loads. 

d)  to  improve  the  technology  in  shear  wall  design. 

e)  how  to  reduce  the  data  to  be  manageable  by  the  micro- 
computer. 


CHAPTER  2 
REVIEW  OF  PAST  WORK 

The  design  of  high  rise  buildings  has  developed  from  the 
bearing  wall  types  to  rigid  steel  frames,  thus  causing 
curtain  wall  buildings. 

The  design  of  shear  walls  has  been  the  subject  of 
considerable  discussion  since  1950' s.  The  simplest  way  of 
designing  shear  walls  is  to  consider  the  shear  wall  as  an 
independent  member  and  design  it  as  a  vertical  cantilever. 
While  this  method  has  the  advantage  of  simplicity,  it  is  not 
universally  valid,  depending  on  the  nature  of  the  column  and 
girder  frame  provided.  The  design  of  the  shear  wall  to  resist 
all  lateral  loads  as  a  cantilever  is  very  conservative  with 
respect  to  the  shear  wall,  while  it  may  be  unconservative 
with  respect  to  the  frame. 

Gould (  1  >   suggested   a   method   of   proportioning    the 
external  load  between  the  shear  walls  and  the  frames  so  as  to 
satisfy   the  requirements  of  deformation  compatibility  with  a 
minimum  of  mathematical  complication. 
2.1  SHEAR  CONNECTION  TECHNIQUE 
2.1.1  ELASTIC  ANALYSIS 

Coull  and  Choudhury<  2  )  developed  the  curves  to  enable 
design  calculations  to  be  performed  rapidly  by  making  use  of 
the  basic  assumptions  of  the  shear  connection  technique.  This 
method  assumes  both  beams  and  walls  in  coupled  shear  walls 
deflect  equally  with  a  point  of  contraf lexure  at  midspan  of 
the  connecting  beams,   and  the  discrete  system  of   connecting 


beans  is  replaced  by  a  continuous  medium  of  equivalent 
stiffness.  The  behavior  of  coupled  shear  wall  is  expressed  as 
a  second  order  differential  equation. 

Schwaighof er <  3  )  modified  Rosman's  theory.  Rosman's 
theory  is  based  on  the  following  assumptions.  a)  the  floor- 
to-floor  height  is  constant,  b)  the  modulus  of  elasticity 
of  the  walls  and  the  beams  is  constant.  c)  the  cross- 
sectional  dimensions  of  the  shear  walls  and  beams  are 
constant,  and  d)  the  points  of  deflection  of  all  beams  are 
located  at  midspan.  Replacing  the  discrete  system  of  shear 
walls  and  beams  by  an  equivalent  continuous  system  and  using 
the  shear  force  along  the  center  of  the  lamina  as  the 
redundant.  Rosman's  approach  can  be  used  but  appears  to  be 
too  lengthy  for  practical  purposes.  Schwaighof er  therefore 
introduced  a  simplification  in  Rosman's  theory.  The 
characteristic  curve  is  replaced  by  a  straight  line.  This 
results  in  considerable  savings  of  time  in  the  analysis. 

CoullC  4  )  presented  the  method  for  the  analysis  of 
coupled  shear  walls  supported  on  elastic  foundations.  He 
indicated  that  the  general  case  of  a  foundation  conditions 
allowing  both  vertical  settlement  and  rotation  to  occur  could 
be  obtained  by  superposition  of  solution  worked  out. 

Because  the  solutions  allowing  for  vertical  settlement  of 
the  foundation  alone  and  for  rotation  of  the  foundation  alone 
are  only  of  limited  application  to  foundations,  Tso  and 
Chant  5  )  presented  the  problem  of  lateral  loading  of  coupled 
shear  walls  resting  on  flexible  foundation. 


ArvidssonC  6  )  presented  a  method  which  took  into 
consideration  the  effect  of  door  openings  near  the  edge  of 
the  wall.  He  included  corrective  terms  in  mathematical  model, 
thus  removing  the  limitations  of  applicabilities.  Also  he<  7 
)  presented  the  method  for  the  analysis  of  shear  walls  with 
two  rows  of  opening.  He  reduced  from  the  high  degree  of 
static  indeterminacy  of  the  system  to  a  fourth  order 
differential  equation. 

ArvidssonC  8  )  analysed  the  interaction  between  coupled 
shear  walls  and  frames.  Secondary  design  effects,  such  as 
column  shortening  in  the  frame  and  elastic  foundation 
conditions,  were  also  taken  into  account. 

Coull  and  Low<  9  )  presented  for  the  analysis  of  non- 
plane  coupled  shear  walls  subjected  to  lateral  loads  which 
produce  combined  bending  and  torsional  deformations. 

Tso  and  BiswasC  10,  11  )  presented  the  analysis  of 
nonplanar  coupled  shear  walls  using  continuum  method.  Both 
bending  and  torsional  deformation  were  taken  into  account. 
The  shear  walls  were  treated  as  thin-walled  beams  of  open 
section  and  the  connecting  beams  between  shear  walls  were 
replaced  by  an  elastic  medium  of  known  stiffness  properties. 
The  total  stiffness  of  the  coupled  shear  wall  system  was  then 
equal  to  the  sum  of  the  stiffness  of  the  shear  walls  and  the 
elastic  medium. 
2.1.2  PLASTIC  ANALYSIS 

Winokur  and  Gluck(  12  )  presented  an  ultimate  strength 
method  to  analyze  coupled  shear  walls.  The  collapse  mechanism 


was  assumed  to  have  plastic  hinges  at  the  points  of  contact 
of  the  connecting  beams  and  the  shear  wall  and  at  the  bottom 
of  the  latter. 

Paulay<  13  )  emphasized  that  cracking  might  have 
considerable  effect  on  the  distribution  of  internal  actions. 
The  elasto-plastic  behavior  was  assessed  by  an  approximate 
analysis,  which  also  yielded  the  plastic  rotations  of  the 
coupling  beams.  He  emphasized  ductility  in  the  shear  wall 
component.  On  shear  walls,  which  carry  a  low  intensity  of 
gravity  loading,  the  lateral  load  may  be  large  enough  to 
cause  cracking.  This  is  certainly  the  case  when  the  structure 
is  exposed  to  a  moderate  earthquake.  In  building  frames, 
cracking  under  service  conditions  can  be  expected  to  cause  a 
considerable  reduction  of  stiffness.  However,  relatively 
little  change  in  relative  stiffness  occurs.  Hence,  cracking 
need  not  significantly  affect  the  bending  moment  pattern  for 
a  frame.  In  shear  walls  on  the  other  hand  the  structural 
components  react  in  markedly  different  ways  to  cracking.  By 
means  of  a  step  by  step  procedure  the  postelastic  performance 
of  structure  can  be  traced . 

Coull  and  Choo<  14  )  presented  an  approximate  technique. 
based  on  the  continuum  representation  of  the  system  of 
connecting  beams,  for  the  analysis  of  the  elasto-plast ic 
behavior  of  laterally  loaded  coupled  shear  walls.  The 
assumption  is  made  that  plastic  hinges  form  at  the  ends  of  a 
connecting  beam,  and  the  substituted  continuous  medium,  when 
the  end  moments  reach  a  limiting  value. 


Gluck<  15  )  presented  an  elasto-plastic  procedure  of 
lateral  load  analysis  of  coupled  shear  walls  based  on  the 
continuum  approach.  Due  to  limitations  on  the  rotational 
ductility  factor,  full  plastif ication  with  height  of  the 
lamina  system  is  very  rarely  possible.  The  ultimate  state  of 
the  structure  being  associated  with  an  incomplete  collapse 
mechanism,  a  plastif ication  on  part  of  the  height  was  to  be 
considered . 
2.2  FRAME  METHOD 
2.2.1  ELASTIC  ANALYSIS 

MacLeod (  16  )  showed  how  three-dimensional  behavior  of  a 
shear  wall  structure  could  be  considered  in  analysis  without 
using  a  full  space  frame  approach.  A  plane  frame  analysis 
involving  line  elements  with  rigid  ends  was  used  to  idealize 
the  various  wall  conditions. 

MacLeod (  17  )  developed  the  complete  stiffness  matrix, 
showing  how  the  shear,  bending  and  axial  deformation  of  the 
flexible  part,  rotational  springs  and  the  rigid  parts  are 
taken  into  account.  He  also  described  a  method  for  carrying 
out  the  back  substitution  for  member  actions. 

Heidebrecht  and  Swift <  18  )  presented  a  method  for 
analysing  asymmetrical  coupled  shear  walls  which  takes  into 
consideration  the  warping  of  the  shear  wall  piers  and  the 
coupling  of  the  piers  by  floor  beams  and  slabs. 

MacLeod  and  HosnyC  19  )  presented  frame  analysis  of  shear 
wall  cores*  He  dealt  with  any  type  of  rectangular  layout, 
e.g.,   unsymmetrical   cores,   multicell  cores,   variation   of 


properties   with  height,  etc..   He  took   into   consideration, 
torsional,   lateral,  and  vertical  load.  Supports  can  be  fixed 
or  elastically  restrained. 
2.2.2  PLASTIC  ANALYSIS 

Nayar  and  CoullC  20  >  presented  a  linear  stepwise 
elastic-plastic  analysis  based  on  the  wide-column  frame 
analogy.  Incremental  loads  are  applied  until  a  plastic  hinge 
forms  at  the  most  heavily  stressed  section.  The  assemblage  of 
the  elastoplastic  stiffness  matrix  is  organized. 

2.3  FINITE  ELEMENT  METHOD 

MacLeod  (  21  )  ,  Giri javallabhan  C  22  )  ,  Al-Mahoidi  and 
NilsonC  23  >  presented  the  finite  element  method  in 
determining  stress  distributions  in  shear  walls  with  openings 
and  in  predicting  more  accurate  values  for  the  bending 
moment,  and  axial  and  shear  forces  which  act  upon  the  beam 
when  the  shear  walls  are  subjected  to  boundary  loads.  MacLeod 
used  the  polynomials  of  second  order  as  displacement 
functions  which  were  chosen  by  a  trial  and  error  process  on 
the  basis  of  the  boundary  compatibility  criterion. 
Researchers  mentioned  above  experimented  only  the  homogeneous 
materials  and  didn't  analyse  the  shear  walls  with 
reinforcement . 

2.4  DYNAMIC  ANALYSIS  AND  SIMULATED  SEISMIC  ANALYSIS 

PaulayC  24  )  presented  the  behavior  of  coupling  beams 
of  reinforced  concrete  shear  walls.  For  a  coupled  shear  wall 
structure  to  develop  its  ultimate  strength  with  modest 
overall   ductility,   the   coupling  beams  need   to   possess   a 


considerable  ability  to  deform  plastically.  To  ensure 
maximum  ductility,  diagonal  tension  failure  in  the  beams  must 
be  avoided.  When  diagonal  tension  failure  in  the  beams  was 
prevented  by  adequate  stirrup  reinforcement,  the  ultimate 
flexural  strength  of  the  coupling  beams  could  be  predicted. 
To  ensure  the  utilization  of  the  maximum  available  ductility, 
the  shear  strength  must  not  be  permitted  to  govern  the 
design .  The  loss  of  stiffness,  after  the  cracking  of  the  test 
beams,  was  considerably  more  than  the  loss  encountered  with 
normally  proportioned  reinforced  concrete  flexural  members. 

PaulayC  25  )  presented  the  analysis  of  spandrel  beam 
under  the  simulated  seismic  loading.  In  deep  spandrel  beam 
(span  to  depth  ratio  is  less  than  2) ,  after  cyclic  loading  of 
high  intensity,  enough  to  cause  the  onset  of  yield,  only  85  X 
of  the  flexural  strength,  as  predicted  by  a  conventional 
ultimate  load  analysis,  can  be  expected  to  develop.  To  ensure 
the  maximum  ductility  obtainable,  a  shear  failure,  by 
separation  along  main  diagonal ,  must  be  suppressed . 

PaulayC  26  )  presented  the  ductile  behavior  of  coupled 
shear  walls.  In  coupling  beams,  it  was  evident  that  the 
energy  absorbed  by  the  diagonal  reinforcement  was 
considerably  larger  particularly  at  low  loads.  He  showed  that 
the  extra  energy  absorbed  at  the  instrumented  lengths  of  the 
diagonal  reinforcement  was  approximately  50  *e  of  the 
corresponding  total  energy  absorbed  in  usually  reinfored 
beams.  The  use  of  ductile  coupling  beams  enabled  a 
considerable  portion  of  the  total  energy  to  be  dissipated   in 
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the  coupling  syatem.  This  dissipation  was  dispersed  over  the 
full  height  of  the  structure  and  it  resulted  in  a  relief  of 
the  plastic  hinges  at  the  base  of  the  walls. 

Paulay,  Priestly,  and  5ynge(  27  )  analysed  the  ductility 
in  earthquake  resisting  squat  shear  walls.  They  experimented 
squat  shear  walls  with  a  height  to  length  ratio  of  less  than 
2.  They  examined  the  possibilities  of  achieving  an  acceptable 
level  of  energy  dissipation  in  squat  shear  walls,  mainly  by 
flexural  yielding.  They  showed  sliding  shear  along  the  base 
is  the  most  significant  cause  of  loss  of  stiffness  and 
strength  and  consequent  reduction  in  ability  to  dissipate 
energy.  Crossed  diagonal  reinforcement  in  two  directions, 
effective  in  both  tension  and  compression,  was  found  to 
considerably  improve  seismic  response  of  squat  shear  walls. 
They  presented  how  to  determine  a  desirable  amount  of 
diagonal  shear  reinforcement  to  control  sliding  and  necessary 
horizontal  shear  reinforcement  to  prevent  diagonal  tension 
failure. 
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CHAPTER  3 
METHOD  OF  ANALYSIS 

The  use  of  finite  element  method  as  a  tool  of  stress 
analysis  is  a  computer  based  technique  for  establishing  and 
solving  the  governing  equations  for  complex  problem  in  a  very 
effective  way.  In  this  method,  the  actual  continuum  or 
structure  is  represented  as  an  assemblage  of  finite 
elements.  These  elements  are  interconnected  at  nodes.  Since 
the  functions  of  displacements,  strains  and  stresses  are 
unknowns  ,  the  displacement  components  within  an  element  are 
approximately  characterized  by  the  unknown  parameters  of  the 
corresponding  displacement  components  and/or  their  partial 
derivatives  at  the  selected  nodal  points  of  the  element  by 
the  use  of  interpolation  functions,  the  shape  functions'  28, 
29,30,31  >.  Using  these  expressions,  one  can  calculate  the 
strain  components  and  use  the  obtained  results  to  find  the 
desired  stresses  and  the  corresponding  strain  energy  by  the 
method  of  numerical  integrations. 

By  the  use  of  the  principle  of  virtual  work,  or  the 
stationary  conditions  of  the  total  energy  in  each  element  of 
the  structural  system,  the  relationships  between  the 
generalized  nodal  displacement  components  and  the  generalized 
force  components  are  obtained.  When  the  elements  are 
assembled  together  to  form  the  entire  structure  of  interest, 
a  matrix  equation  for  finding  the  nodal  displacement 
parameters  is  obtained.  Also  appropriate  boundary  conditions 
should  be  used  in  the  establishment  of  the  matrix  equation  to 
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ensure  the  existence  of  the  unique  solution  of  the  problem. 
A  list  of  the  standard  steps  of  the  finite  element  analysis 
is  : 

i>  Discretization  of  the  structure 
ii)  Selection  of  proper  interpolation  functions 
iii)  Calculation  of  element  stiffness  matrices  and  load 
vectors  by  the  use  of  numerical  integration 
iv)  Assemblage  of  element  stiffness  matrices  to  obtain 
overall  global  stiffness  matrix  and  loading  vector 
v)  Enforcement  of  appropriate  boundary  condition 
vi)  Solution  for  the  unknown  nodal  displacements 
vii)  Computation   of  strains  and  stresses  at  the  points 
of  interest 
Because  of  the  limitation  of  the  memory  space  of   Zenith 
100   micro-computer,   the  matrix  equation  was  further  reduced 
by   the  use  of  antisymmetric  property  of  the  displacement   of 
the   particular  problem  of  interest.   The  detailed   steps   of 
these  operations  will  be  reported  later. 
SHAPE  FUNCTIONS  AND  DEGREES  OF  FREEDOM 

The  shear  wall  of  interest  was  treated  as  a  plane  stress 
problem  since  the  stress  component  in  the  direction  of 
thickness  is  equal  to  zero  almost  everywhere. (  30,31  ).  There 
are  many  shape  functions  for  solving  2-dimensional  elasticity 
problems  available  such  as  constant  strain  element,  element 
with  quadratic  interpolations,  Hermitian  interpolation  and  ra- 
th order  Hermitian  interpolating  functions, etc . .  Through  the 
numerical    experiments   during   the   investigation   of    the 
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properties  of  the  Hermitian  interpolating  function,  the 
author  discovered  that  the  strain  energy  density  function 
obtained  by  the  use  of  Hermitian  interpoation  functions  to 
represent  the  displacement  function  of  elements  under  the 
application  of  rigid  body  displacement  is  zero  only  when  the 
domain  of  element  is  a  rectangle  (  See  Fig.l-Fig.lb  and 
program  "ENERGY"  in  APPENDIX.  Note  that  this  plot  shows  that 
strain  energy  value  is  zero  everywhere  when  the  given  nodal 
displacements  are  specified  according  to  a  rigid  body 
displacement  ).  Otherwise.  for  example  an  arbitrary 
quadrilateral  element,  the  obtained  strain  energy  density 
functions  for  rigid  body  displacements  are    not  zero. 

In  order  to  reduce  the  number  of  unknowns,  it  is  not 
suitable  to  use  constant  strain  element.  The  use  of  the  first 
order  Hermitian  element  can  provide  quadratic  stress 
expressions  over  each  element  and  the  obtained  stress 
expression  over  the  entire  continuum  will  be  continuous.  It 
is  natural  to  use  rectangular  shaped  domain  for  the  analysis 
of  shear  wall -typed  structures  which  will  produce  consistent 
strain  distributions  according  to  theoretical  expressions  as 
well  as  numerical  experiment  performed  by  the  use  of  the 
Zenith  100. (  See  Fig.  2  ) 

The  Hermitian  shape  function  for  this  particular  case  can 

be   constructed   according  to  the  products  of   the   following 

functions: 

L-X     X(L-X) 

Fi<X)=  *  CL-2X)  (  1  ) 

L         L3 
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X      XCL-X) 

F2<X)=  *   (2X-L) 

L        L3 


F3(X>= 


X<L-X>2 
L2 
<L-X)X2 


F4<X)=  - 


<  3  ) 


(  4  ) 


Note  that  Fi<X)  and  F2<X)  are  exactly  the  influence  lines 
of  Ra  and  RD  under  the  application  of  a  moving  unit  load  on  a 
fixed-fixed  beam.  The  F3<X)  and  F4<X>  are  the  influence  lines 
of  the  fixed  end  moments  Ma<X)  and  Mfc>(X>  as  shown  in  Fig.  2  . 

Let  the  nodal  displacement  components  at  the  i-th  node  be 
arranged  according  to  the  following  vectors. 


(  Ui  } 


Ou/3x>i 
<9v/©x>i 
C8u/9y>i 
Ov/ay>i 


>• 


vi 

<ux>i 

<vx>i 

<uy)i 

l<vy>ij 


C  5  ) 


and,  locally  (with  respect  to  each  element)  the  order  of  the 
4  nodes  of  the  e-th  element  be  ordered  according  to  Fig,  2  , 
the  components  of  the  nodal  displacement  parameters  of  this 
element  expressed  in  a  row  vector  form   reads 

{  ue  }  =  lul«   vl»   Uxl»  vxl>  uyi"  vvl«  u2«  v2>  ux2'  vr2> 

uy2»  vy2»  U3,  vg,  Ux3#  vx3*  UY3»  vy3#   U4„   V4, 

"x4»  vx4,  uy4,  vy4j  (  6  1 

According  to  this  arrangement,  the  24  interpolating  functions 

corresponding   to   these  nodal  parameters  are  constructed   as 
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F3(x) 
F2(x) 


p-1 


Fi(y)    F3(y)    F2(y)      F4(y) 


»•( 


=L 


)Mb 


F4(x) 


R  (x)=  JC2L-  +  ^CL-x)  (L-2x) 
L       L2      L 


»  (x).  X    +  x(L-x)   (2x-L) 
L         L2      L 


Ma(x)=;t(L-x)2       )=  x2g-x) 


L  1/ 

Fig.  2  The  local  ordering  system  and  the  shape  function 
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follows . 

»1  =  Fi<X>»Fi<Y> 

N2  =  Ni 

N5  *    Fi<X)»F3<Y> 

N6  ■  N5 

Ng  =  F4<X)»Fi<X) 

N10  =  Ng 

N13  -  F2(X)»F2<Y' 

Nx4  =  Ni3 

N17  =  F2<X)«F4<Y> 

Nia  »  N17 

N21  *  F3<X)»F2<Y> 

N22  =  ^21 


<  7  ) 


N3  =  F3<X)»Fi<Y.) 

N4  ■  N3 
N7  =    F2<X>*Fi<Y> 

NS  ■  N7 
Nil  "  F2(X)»F3(Y> 

N12  -  N11 
N15  =  F4(X)»F2<Y) 

Nl6  ■  Nl5 
Nig  '    Fi<X)»F2<Y) 

N20  "  N19 
N23  ■  Fi<X)»F4<Y) 
N24  "  N23 

These  interpolating  functions  are  used  in  the  3-dimensionol 
plot  program  <  see  appendix  )  and  plotted  in  Fig.  3-14  for 
a  typical  rectangular  shaped  element. 

The  equation  which  describes  the  displacement  of  the  element 
is 

'Ni  0  N3  0  N5  •  •  ■  •   N21  0  N23   0  ' 

(  ue  >  <  8  > 
„  0   N2  0  N4  0  •  •  •  ■    0  N22  0   N24  , 

For  simplicity.   the  displacement  components  of  the   e-th 

element  can  be  rewritten  as 


EN]  (  u«  } 


<  9  > 


The   definition   of   C   N  ]  and  t  ue  )  can   be   redefined   by 
comparing  Eq . (  5  )  and  Eq.(  7  ). 
STRAIN  COMPONENTS 
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Fig. 3  First  order  Hermitian  shape  function  using  nodal  w  of  node  1 
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Fig.  4  First  order  Hermitian  shape  function  using  nodal  w  of  node  2 
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Fig.  5  First  order  Hermitian  shape  function  using  nodal  w  of  node  3 
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Fig.  6  First  order  Hermitian  shape  function  using  nodal  w  of  node  4 
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Fig.  7  First  order  Hermitian  shape  function  using  dw/dx  of  node  1 
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Fig.  8  First  order  Hermitian  shape  function  using  dw/dx  of  node  2 
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Fig.  9  First  order  Hermitian  shape  function  using  dw/dx  of  node  3 
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Fig.  10  First  order  Hermitian  shape  function  using  dw/dx  of  node  4 
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Fig.  11  First  order  Hermitian  shape  function  using  dw/dy  of  node  1 
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Fig.  12  First  order  Hermitian  shape  function  using  dw/dy  of  node  2 
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Fig.  13  First  order  Hennitian  shape  function  using  dw/dy  of  node  3 
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Fig.  14  First  order  Hermitian  shape  function  using  dw/dy  of  node  '4 
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For  the  stress  analysis  of  shear  wall.  it  is  considered 
as  a  plane  stress  problem.  The  strain  components  are 
calculated  by 

{  s  }  =  [  B  1  (  ue  }  <  10  i 

where 

'a 


f3*  1 


7  <  11  );  C  B  1  = 


TxY 


dx 
0 

a 
I  ay 


ay 
a 

ax 


The  detailed  expression  of  C  B  ]  reads 

3N21 


[  B  3  = 


(   3NJ        3N3 

0  0  • 

3x        8x 

3N2  3N4 

0    0  

3y  ay 

3N2  9Ni  8N4  3N3 


3N2? 


0 


0 


ax 

o 


ax 

3N22      3N24 
0   


3y        3y 
3N22  3N21  &N24  3N23 


8y   3x   3y 


ax 


ay 


ax 


ay 


ax  > 


<  13  1 


In  case  of  linear  elastic  isotropic  plane  stress  problem,  the 
stress-strain   relations   are  given  by  Hooke's  law  . 


{  S    ) 


[  d  :  i 


(  14  ) 


where 


C  D  1  = 


1 
0 


o 
1   -  n 


<  15  ) 


(  S  } 


r«x  ] 

<3y 
TXy 


<  16  ) 
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The  strain  energy  of  the  element  can  be  expressed  as  followa 


r  (  x  .  y  )  t  dx  dy 


(  17  ) 


where  the  T*  is  the  strain  energy  density  function 


r  =  —   1_  «  J  C  D  3  < 
2 


<  18  1 


and  t  is  the  function  of  thickness  of  a  wall. 

By  the  substitution  of  Eq .  (  10  >,  one  has  the  strain  energy 

1 


U 


L  "eJ 


t  [  B  3T  [  D  3  C  B  3  dA 


{  ue 


)  (  19  ) 


Note   that  by  the  use  of  the  definition  of  stiffness   matrix, 

the   nodal  force  (  Fe  )  =  C  K  3  (  ue  ) .   Therefore   the   work 

done  by  the  nodal  force 

1 

W   =  .  ue  |   C  K  3  {  ue  )  <  20  ) 

2 

The  expression  of  stiffness  can  be  yielded  by  comparison  of  U 
and  W. 


C  K  3 


C  B  3  ■   t  D  3  [  B  3  dA 


<  21  > 


CALCULATION  OF  GENERALIZED  NODAL  FORCES 

In    terms    of   nodal   forces   and   nodal    displacement 
components,  the  work  done  by  the  e-th  element  is 


W  =  |   ue  _|   [  Fe  ) 

2 


(  22  ■) 


The  work  done  due  to  the  external  loading  intensity  of  fx<s> 
and  fyCs),  which  are  the  x  and  y  components,  respectively,  at 
s  over  the  i-th  loaded  side  si  can  be  calculated  by 
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Wi 


1 
2  J 


.  u 


fx(s) 


fy<») 

The   work   done   by   the   external  load   applied   on   the   Na 
boundaries  of  the  e-th  element  is  the  summation  of  W^ . 


Vs 


Ne   1 

S  

i  =  l   2 


l_v 


t  f  )  ds 


<  24  1 


By  the  substitution  of  Eq.<  8  )  into  Eq . (  24  ),  one  gets 


we 


Ne   1 

2  

i  =  l   2 


C  NT  ]  [  f  )  ds 


(  25  ) 


By   the  use  of  the  principle  of  virtual  work,   one   can   show 
that  the  i-th  component  of  the  C  Fe  )  vectors  is 
e      ^e 


8ui 


Thus 


•      Ne 

I  Fj  ]  =  Z 
i  =  l 


[  N  ]T   t  fs  3  dsi 


<  26  .) 


<  27  1 


Si 


The  vector  of  generalized  nodal  force  components  is  obtained. 
NUMERICAL  INTEGRATION  AND  THE  USE  OF  GAUSSIAN  QUADRATURE 

There  are  several  methods  available  for  the  numerical 
evaluation  of  definite  integrals.  The  use  of  Gaussian 
quadrature  has  been  proved  to  be  most  useful  for  finite 
element  work  because  the  integration  of  a  n-point  Gaussian 
quadrature  can  provide  exact  integration  for  <2n-l)-th 
polynomials.  In  two-dimensional (  rectangular  )  regions,  the 
result  can  be  obtained  by  using  Gaussian  quadrature  with 
respect  to  x  and  y  direction  iteratively. 
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n    n 
I  -    Z   S  «j  »j  f  (  xi  ,  yj  )  a-b 
i»l  3=1 


<  2a  ) 


where  Wj.  and  Wj  ara  tha  weighting  factors  at  the  i-th  and  J- 
th  sample  points  and  Xi  and  yi  ara  tha  x  and  y  values  of  the 
i-th  and  j-th  sample  points. (  30  )  Tha  meaning  of  a  and  b  ara 


In  tha  computer  programs 
developed,  9  Gaussian  points 
over  tha  rectangular  domain 
wara  used.   Thay   ara  shown  in 
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STIFFNESS  OF  THE  REINFORCED  STEEL  BAR 

y  Tha   contribution   of   tha 

..t     Cxb-yb) 

rainforcad    staal    shown   in 
Fig. 16  to  tha  stiff nass  matrix 

of  an  element  is  considered  as 

—  X 

follows.  Sines   tha  reinforced 

Fi3-i6  A  steel  bar  of  an  element 

steal   bar   is  embedded  in  the 

concrete,   it  is   suitable   to  assume   that  the  displacement 

of  tha  rainforcad  staal  bar  at  any  point  is  characterized  by 

tha  aaae  displacement  functions  of  the  element  which  contains 

tha  staal  bar  of  interest,  i.e.. 
Nu 


l>  N, 


t  u«  ) 


<  29  ) 


C  B  J  (  u«  } 


(  30  ) 
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The   strain   of  interest  for  a  steel  bar   is   only   the 
strain  in  the  axis  of  the  bar.  Therefore 


I  <*e  U 


cos20  sin29   — sin20 


5} 

VxY 


L'sj 


t 


C  31  ) 


i.e., 

<  «9  >s  =   j_  T5J   [  B  J  {  ue  )  <  32  ) 

To   avoid  addition  of  concrete  again,   which  is   replaced   by 
steel,  the  Young's  modulus  of  steel  E5  is  replaced  by 

E*  =  Es  -  Ec  (  33  j 

where  Ec  is  the  Young's  modulus  of  concrete.  Therefore 
m 
Ss  =  (  Es  -Ec  >  (  &q   )s  (  34  ) 

The   additional   strain  energy  of  the  reinforced   steel 
bar  of  area  As  will  be 


Ui 


As 

2 


|_ua  J  CB]T  (Ts>  E"  LT,  j   [B:  Cue)  dl 


(  35  ) 


Thus  the  additional  stiffness  due  to  the  s-th  steel  bar   will 
be 


K. 


C  B  3T  C  Ts  )  E*  L  T5  J  [  B  3  dl   <  36  ) 


By   adding   the  summation  of  the  contribution  of  all  bars   to 

the   stiffness   of  concrete  of  the  e-th  structural  element 

the  combined  stiffness  matrix  becomes 
e       e  e 

C  K  :=  [  Kconcrete  }  *  C  K3t  3  (  37  ) 

where   C   Kst   3   denotes  stiffness   of   steel   bars  and   the 

expression  is 


C  Kst  ]  =  2  CK53=  S  As 

3=1        S=l 


I  B  IT(  Ts  )E*  1_T5  j[  B  )dl   C  38  ) 
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CHAPTER  4 
NUMERICAL  EXPERIMENT 

The  dimension  of  model  shear  walls  analysed  are  the 
same  as  that  employed  in  a  test  unit  experimented  by  Paulay< 
27  > .  The  overall  dimension  of  the  model  shear  walls  is  shown 
in  Fig.  17. 

The  model  shear  wall  is  divided  into  52  rectangular 
elements  with  71  nodal  points  as  shown  Fig.  18.  The  3  types 
of  shear  walls  with  different  array  of  reinforcement  as  shown 
in  Fig.  19  and  20  are  analysed.  Properties  of  each  shear 
walls  are  presented  in  Table.  1. 


WALL 
TYPE 

SECTION 

fc 
MPa 

REINFORCEMENT 

HORIZONTAL 

VERTICAL 

DIAGONAL 

TYPE  1 
TYPE  2 
TYPE  3 

rectangular 
rectangular 
rectangular 

27.2 
27.2 
27.2 

D14-120 
D12-120 

D12-140 
D12-140 
4-D12  for 
each  edgej 

D12-200 

The  reinforcement  used  in  type  1  is  exactly  the  same  as  that  used 
by  Paulay  (27) 

MODULUS  OF  ELASTICITY 

For  the  purpose  of  analysis  of  the  shearwalls,  three 
different  kinds  of  modulus  of  elasticity  of  concrete  were 
used  as  shown  in  Stress-strain  curve  of  Fig.  21.  The  first 
value  used  is  the  modulus  of  elasticity  which  is  calculated 
according  to  the  formula  (8-5-1)  given  in  the  ACI  code  (34>. 
For   clarity,   first  value  will  be  called  as  initial  modulus 
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of    elasticity    <    35    )  .    For    f  c    ■    3944    psi 


5.64*10  O.0O22 


Eco    =    57000    J  f c 

=  3579700  psi (  24700  MPa ) 

The  maximum  strain  «£,  at  point 

2  shown  in  Fig. 24,  of  this  concrete 

is  estimated  as  follows. 

2fc 

«2  ■  =  2.203  *  10-3 

E 


Fig. 21    Stress-Strain    Curve     The         stress-strain      curve      can       be 

of  concrete 

interpolated  by  the  use  of  parabolic  curve. 


S(  «  )  =  Ai«e  *  A2*«2 


<  39  ) 


A2 


S(  •  )  -  Ai< 


C  40  ) 


Differentiation  of  Eq.C  39  )  with  respect  to  s  yields 
d6(«) 


=  Ai  *  2  A2  « 


<  41  ) 


Initial  Modulus  of  elasticty  can  be  found  by  substituting  «=0 
into  Eq.  C  41  ) . 

Ai  *  2A2*0  =  3579700 
Therefore, 

Ai  -  3579700  psi 
The   coefficient  A2  can  be  obtained  by  substituting  the 
value  of  «2  and  Ai  into  Eq .  <  40  ) . 
3944-3579700».0022 


A2  = 


(0.0022)2 
812260000 
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Therefore  the  function  of  stresa-strain  curve  is 

S<  e  )  =  3579700*e  -812260000«e2  (  42  ) 

When   the   stress  S<  e  >  •  0.45  fc„   the  strain  can  be 

obtained  by  substituting  0.45fc  =  1774.8  psi  into  Eq .  <  42  >, 

Thus, 

1774.8  =  3579700«i  -  812260000^2 

e  can  be  obtained  by  solving  the  above  equation.  Therefore. 

«1  =  5.6935  *  10~4 

Secant    modulus   of   elasticity   Eci   can   be   obtained    by 

substituting  the  above  «i  into  the  following  equation. 

0.45  fc 


Eel  - 


1774.8 


=  3117000  psi  (  21500  MPa  > 


5.6935»10"4 

Another  secant  modulus  of  elasticty  corresponding  to  fc  -      fi- 
can  be  obtained  by  the  following. 

Ec2  "  


3944 


=  1793000  psi  <  12400  MPa  ) 


0.0022 

Modulus   of  elasticity  of  reinforcement  used  was  29«106 
psi  as  described  in  <8-5-2)  of  ACI  code. 
RESULTS  OF  NUMERICAL  EXPERIMENT 

Three  model  structures  with  different  reinforcement  were 
studied  using  three  different  moduli  of  elasticity  of 
concrete  for  each  model.  Lateral  load  of  200  kN  was  applied 
on  the  upper  left  top  of  the  model  as  shown  in  Fig.   17.   The 
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flow   chart   of   computer  program  as  shown   in   Fig.   33   was 
followed  and  the  programs  explained  in  Chapter  5  were  used. 

Nodal  dispacements  of  each  node  were  found  and  two 
dimensional  plotting  of  dispacement  shape  as  shown  in  Fig.  23 
was  drawn  by  the  computer  program . 

By  the  use  of  the  known  displacements,  the  trajectories 
of  stresses  as  shown  in  Fig.  22  can  be  found.  The  steel 
reinforcements  inside  the  shear  walls  can  be  placed  by 
following  the  trajectories  of  the  stresses.  For  simplicity. 
the  diagonal  reinforcement  was  used  in  wall  type  3. 

By  the  use  of  the  known  displacements  of  each  node, 
strain  and  stress  of  the  points  of  interest  were  calculated. 
The  maximum  strains  of  model  structure  were  found  and 
presented  in  Table.  2.  Also  the  three  dimensional  displays  of 
normal  stresses,  shearing  stresses  and  stresses  of 
reinforcement  are  shown  in  Fig.  24  -  32.  For  simplicity,  the 
results  of  shear  wall  type  1  (  Fig.  24-30)  and  shear  wall 
type  3  <Fig.  31-32)  using  the  intial  modulus  Ec0  are 
presented . 

The  comparison  of  allowable  shear  force  between  the 
estimated  value  and  the  experiment  data  done  by  Paulay  (  27  ) 
are  presented  in  Table.  3.  The  estimated  allowable  shear 
force  is  2.4  times  higher  than  that  of  the  experimental  data. 
The  estimated  cracking  load  is  smaller  than  that  of  the 
experiment  data. 

A  comparison  of  the  ultimate  shear  force  between  the 
estimated   and  the  experiment  data  is  presented  in  Table.   4. 
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The   estimated   ultimate  shear  force  i3  2  times   higher   than 
that  of  the  experiment  data. 

If  the  initail  modulus  of  elasticty  at  zero  stress 
level  is  used  in  the  finite  element  analysis,  the  estimated 
value  of  shear  force  that  could  start  crack  of  concrete  is 
calculated  as  follows. 

<<sc/10)»Pgppiied 


pcrack  - 


strain 
(0.0022/10)»200 


=  389  kN 


1.130213»10-4 
When  the  highest  compressive  stress  of  shear  wall  reaches 
0.45fc,  the  stress  of  rest  parts  of  shear  wall  is  less  than 
0.45fc.  Therefore,  Young's  modulus  of  whole  structure  is 
between  the  value  calculated  according  to  initial  modulus  at 
zero  stress  level  and  that  calculated  at  the  secant  modulus 
of  0.45fc  .  This  suggests  the  average  value  of  allowable  P 
predicted  by  these  two  models  should  be  more  accurate  in 
comparison  with  the  experimental  data.  According  to  this,  the 
average  allowable  P  is  equal  to 
389  *  342** 


365.5  kN 


2 
The   error  is  less  than  1  *  difference  between  the  value   of 
0.45  PexPerimented  and  this  estimation. 

»«  This  value  P=342  kN  is  reported  on  PP.52  which  is  an 
estimated  P  value  that  could  start  the  hair  crack  of  concrete 
in  the  wall. 
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CHAPTER  5 
THE  COMPUTER  PROGRAMS  DEVELOPED 
ELEMENT  STIFFNESS  MATRIX  :MAIN  PROGRAM 

This  program  calculates  the  stiffness  matrices  and  the 
load  vectors  of  various  element.  Input  data  of  this  program 
include  modulus  of  elasticity  and  Poisson's  ratio  of 
concrete*  coordinates  and  node  numbers,  locations  of  steel 
bars,  steel  area.  Young's  modulus  of  steel,  Gaussian  weights 
and  locations,  and  external  loads.  Matrix  [  D  ]  is  calculated. 
Element  stiffness  matrix  and  load  vector  are  calculated  and 
stored  in  STDATA  data  file  and  FCDATA  file  respectively. 
ASSEMBLING  ELEMENT  STIFFNESS  MATRIX :ASSEM  PROGRAM 

This  program  gets  the  <24»24>  stiffness  matrix  of  each 
element  from  a  data  file  named  "STDATA"  and  assembes  element 
by  element  to  form  a  lower  triangular  global  stiffness  matrix 
of  the  structure.  The  data  assembled  by  this  program  are 
stored  row  by  row  in  a  data  file  "STDATAN".  This  is  done  by 
the  use  of  the  global  order  of  the  4  nodes  of  each  element 
assigned  to  the  nodes  of  a  fixed  local  ordering  system  shown 
in  Fig.  34. 
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Fig.  34  Local  ordering  system 
for  a  typical  element 


S= 


Fig.  35  Partitioned  matrix 
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For  example,  let  G3.,  G2.  G3  and  C4  be  the  order  number 
of  the  nodes  lj,»l,  12  =  2,  13  =  3  and  14  =  4  as  shown  in  Fig.  35, 
then  the  location  of  each  submatrix  of  the  stiffness  matrix 
of  the  structural  element  can  be  associated  to  the  required 
location  in  the  global  stiffness  matrix  explained  as  follows. 

In  this  partitioned  stiffness  matrix,  each  submatrix  has 
(6*6)  elements.  The  first  subscript  refers  to  the  local  order 
number  of  the  node  to  which  the  six  nodal  force  components 
are  of  interest,  and  the  second  subscript  refers  to  the 
related  nodes  that  its  six  displacements  are  to  be  used  as 
its  operand.  Note  that  the  row  number  of  the  global  stiffness 
matrix  is  related  to  the  corresponding  nodal  force  component, 
and  the  column  number  is  related  to  the  corresponding  nodal 
displacement  parameter  when  they  are  arranged  according  to 
the  order  of  the  nodes  of  the  structure.  For  simplicity  of 
explanation,  only  the  procedure  to  add  S13  submatrix  into  the 
global  stiffness  matrix  t  K  ]  will  be  explained. 


«. 
«* 
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. 

Fig. 36  Global  stiffness  matrix 


Since  the  global  node  number  of  li  and  I3  are  G]_  and  G3, 
therefore 
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R6  =  Gi»6 

Ci  =  <G3-1)»6 

C3  =  G3»6 
The  submatrix  S13  should  be  superposed  into  the  global  matrix 
[  K  ]  at  the  location  of  a  submatrix  of  K  <  Rl,  Ci  )  and 
KCRg.Cg),  which  are  located  on  the  upper  left  corner  and  on 
the  lower  left  corner  respectively  as  shown  by  the  shaded 
block  in  Fig.  36.  The  superposition  of  all  other  submatrix 
can  be  done  similarly.  Note  that  only  elements  belonging  to 
the  lower  triangular  part  of  the  global  stiffness  matrix  will 
be  stored  in  the  data  file. 

This  program  also  assembles  the  element  force  vector 
into  the  global  force  vector.  This  program  gets  the  (24*1) 
force  vector  of  each  element  from  a  data  file  named  "FCDATA" 
and  assembles  element  by  element  to  form  a  global  force 
vector.  The  data  are  stored  in  a  data  file  "FCDATAN". 

By  the  use  of  this  program,  a  complete  system  of 
simultaneous  equation  can  be  expressed  as 

C  K  :  t  u  )  =  t  F  }  <  43  ) 

REDUCED  STIFFNESS  MATRIX :REDUC  PROGRAM 

This  program  reduces  the  dimensionality  of  stiffness 
matrices  to  form  a  reduced  matrix  equation.  Because  the 
structure  is  a  symmetric  shear  wall,  nodal  displacements  of 
all  the  nodes  to  the  right  of  the  structure  can  be  expressed 
in  terms  of  the  corresponding  components  of  the  nodal 
parameter  of  the  corresponding  nodes  on  the  left  half  of   the 
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structure . 

Let  C  C  ]  be  a  connection  matrix  which  is  constructed 
according  to  the  condition  of  asymmetry  of  the  displacement 
field  over  the  corresponding  nodes.  For  example.  the 
conditions  of  u  and  v  at  the  points  <  xi.  yi  )  and  <  -xi,  yi 
)  are  u  <  Xi.  yi  >  ■  u  (  -Xi,  yi  ),  and  v  <  xj..  yi  )  =  -(- 
*i>  Yi5-  If  the  global  order  numbers  for  these  two  nodes  are 
p  and  q,  the  C(p-l) *6*1]  th  row  of  C  matrix  has  value  of  zero 
at  each  column  except  for  the  value  of  1  at  column  (q-l)»6*l. 
On  the  C<p-1>«6*2]  th  row.  only  the  C<q-1)»6*2]  th  column  has 
value  of  1  and  the  value  of  rest  elements  of  this  low  is 
zero. 

After  the  connection  matrix  C  has  been  constructed,  the 
complete  vector  of  the  global  nodal  displacements  (parameters) 
(  u  J  can  be  expressed  in  terms  of  the  reduced  vector  of 
unknown  nodal  displacements  t  ur  }  as  follows. 

[  U  }  •  [  C  1   (  Ur  ]  (  44  ) 

<n»l)     (n»r>   (r»l) 

Substituting  this  equation  into  Eq .  <  43  .)  ,  one  gets 

C  K  ]   C  C  3  {  ur  }  =  (  F  )  <  45  ) 

(n»n)  (n«r)  <r«l)      (n«l) 

Premultiplication   of  Eq. <  45  )  by  the  transposition  of  C  C  ] 

yields 

[  C  3    :  K  ]  t  C  D  {  ur  )  -  C  C  ]    {  F  )        <  46  ) 
<r»n>   (n»n)  (n»r)  (r»l)     <r*n)   <n*l> 

Let   C  Kr  ]  =  [  C  ]    C  K  ]   C  C  ]  <  47  ) 

<r»r ) 

and      {  Fr  )  =  C  C  ]    (  F  }  <  48  > 

Eq.<  46  )  is  reduced  into  the  following  form 
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t  Kr  3  (  ur  )  =  {  Fr  }  <  49  ) 

However,   the   C   C   ]  and  C  C  ]   matrices  did   not   actually 

constructed  in  the  program  "REDUC".   The  concept  was  used   to 

accomplish  the  process  of  CC:  CK3  CC]  and  [C:  (F}  . 

Note   that  particular  attention  was  given  to  take   care 

the   case  that  only  the  lower  triangular  part  of   matrix   was 

saved.   This   introduced  complicated  manipulation  of   element 

location  in  the  development  of  program  "REDUC" .  The  important 

part  of  the  process  are  reported  as  follows. 

S-th    t-Ch 
column  column 


s-th  row-, 
t-th  row* 
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* 

• 

fs 
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k 

. Fig.  37  Lower  triangular  matrix 
Considering  the  cose  that  the  s-th  row  of  C  C  ]T  matrix 


has  the  following  elements 
Cs 


J 


10. 0, ■ • -0.1.0 .a.O. ■ ■ 

s-th  column    t-th  column 


then   the  s-th  row  of  C  K*  ]  is  the  summation  of 
a*  (_Kt  j   .  i.e.  , 

L.Ks_j   "  |_  Ksj   *  a  »  L_Kt_l 


t  >  s 


where 


t  K*  ]  =  C  C  ]T  C  K  ] 


<  50  ) 


and 


<  51  ) 


<  52  ) 


L"sj 


s-th  row   of  K  matrix 


|_  Kt  _J  =  t-th  row   of  K  matrix 
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Note  that  only  the  shaded  part  of  K  matrix  shown  in 
Fig.  37  is  saved  in  the  memory,  and  the  required  element 
<Kij>  for  J    >  i  needs  to  be  calculated  by 

KiJ  "  Kji  <  53  ) 

In  the  program,  K  matrix  is  saved  row-wisely  from  first 
column  to  the  main  diagonal  row  by  row  in  a  column  matrix, 
i.e.,   (  A  >  ,  then  the  sequence  number  of  Kj±    in  {  A  }  is 


(  54  ) 


2 

Therefore 

Kij  =  A  <  m  >.  (  55  ) 

By   using  this  method  of  locating  desired  elements,   one   can 

find  K3  as  follows: 

For  elements  Ksj  with  j  <_  s, 

•        s»(s-l)  t«(t-l) 

Ks3  "  *  < *  J>  *  a  •  A  < ♦  j )    <  56  > 

2  2 

# 
For  elements  Ksj  with  j  >  s  and  J  <.  t, 

»        J"<J-1>  t»Ct-l> 

Kaj  =  A  ( *  s)  ♦  a  »  A  ( *  ])    <  57  ) 

2  2 

* 

For  Kaj  with  J  >  t, 

3*<J-U  j»(j-1) 

Ksj  =  A  < *  s)  *  a  »  A  < *  j)    (  58  ) 

2  2 

*  * 

The  superposition  of  C  Ks  ]  and  a»  C  Kt  3  to  complete  the  s-th 

column  of 

I  C  ]    [  K  1   [  C  ]  =  [  KR  ]  <  59  ) 

can  be  done  likewise.  For  simplicity,  let 

C  R  3  =  C  Kr  ]  (  60  > 
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then  the  a-th  column  of  R  can  be  expressed 

*  * 

I  8S  I  =  I  K5  )  *  a  <  [  Kt  1  <  61  ) 

Note  that  in  the  expression  of  [  C  ]  C  K  ]  =  C  K  ] 
the  elements  in  the  related  columns  have  the  following 
calculations . 

For  t  >  i  and  i  *    s 

»         t«<t-l> 
Kit  ■  A  <  *  i  )  C  62  ) 

2 

»        t*<t-l)                 t«(t-l> 
K5t  =  A  ( *  s)  *  a  «  A  < *  t)    <  63  ) 


For  i  >  t 


»         i»(i-l) 

Kit  "  A  <  *  s  ) 


For  i  >  s 

•         i«<i-l) 

Kis  ■  A  (  *  s  )  (  65  ) 

2 
»        s«<»-l>  t*(t-l> 

KSs  =  A  <  *  s)  *  a  »  A  < ■<■  s>     <  66  ) 

2  2 

Substituting  these  relationships  into  Eq . <  61  ),  one 
can  find  the  s-th  column  of  the  reduced  stiffness  matrix.  The 
superposition  of  CC3  (F)  can  be  done  in  a  straightforward 
way  which  needs  no  further  explanation. 

Because  the  model  structure  as  shown  in  Fig.  18  has  71 
nodes  and  each  node  has  6  nodal  parameters,  the  size  of 
assembled  global  stiffness  matrix  is  (426*426)  and  the  total 
number  of  elements  of  stiffness  matrix  is  181,476.  The  Zenith 
100   microcomputer   cannot  handle  this  amount    of   data   for 
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lack  of  memory  spaces  and  the  data  used  should  be  reduced.  By 
using  the  advantage  of  the  property  of  symmetry  in  the 
stiffness  matrix,  lower  (  or  upper  >  triangular  part 
including  main  diagonal  elements  can  be  stored.  In  this  case, 
the  total  elements  of  triangular  part  are 
6n»  <6n*l) 


2 
6«71«(6«71+1> 


2 
■  90.951 
where   n   denotes   total   node   numbers.   Total   elements   of 
stiffness  matrix  by  the  use  of  "REDUC"  program  are 
6«39»<6«39i-l) 


27,495. 


2 
Percentage   of  data  space  required  between  the   "REDUC" 
program  and  the  usual  program  not  using  reduction  is 
27,495 


*100  =  30.23  ?s 


90,951 

This   program  requires  only  30.23  t  memory  space.   A   lot   of 

computer   memory  spaces  and  C.P.U.   time  can  be  saved.   Force 

vector  are  also  reduced  by  this  program. 

BOUNDARY  CONDITIONS  OF  ZERO  DISPLACEMENTS IMODIF  PROGRAM 

After   assembling   the   element   stiffness   matrix   and 

element   load  vector  of  each  structural  element,   the   matrix 

equation  of  the  entire  structure  can  be  written  as 

I  K  ]  (  u  }  =  (  F  J  (  67  ) 

m»m    m*l      m*l 

Without  enforcement  of  the  boundary  conditions,  this  equation 
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cannot  be  solved  for  t  u  }  since  the  matrix  [  K  ]  will  be 
singular  and  its  inverse  does  not  exist.  The  physical 
significance  of  this  is  that  the  structure  is  free  to  undergo 
unlimited  rigid  body  motion  unless  some  support  constraints 
are  imposed  to  keep  the  structure  in  equilibrium  under  the 
loads.  Hence  some  boundary  conditions  have  to  be  applied 
before  solving  {  u  ) .  The  rows  and  columns  corresponding  to 
the  zero  nodal  parameters  have  to  be  eliminated  in  order  to 
make  the  stiffness  matrix  non-singular.  The  lower  triangular 
part  of  the  stiffness  matrix  is  stored  as  a  single 
dimensional  array.  Except  for  the  main  diagonal  elements 
which  are  unit,  the  rest  of  the  eliminating  part  are  made  to 
be  zero.  Data  created  in  this  program  is  the  same  as  one 
created  in  "ASSEH"  program  except  for  eliminating  rows  and 
columns . 
BOUNDARY  CONDITIONS  OF  ZERO  STRESSES rENBOUN  PROGRAM 

In  addition  to  the  procedure  of  "MODIF"  program,  the 
zero  stress  components  of  all  the  boundary  nodes  are 
considered  in  this  program. 

In  case  of  6*  =  0  on  a  certain  boundary  node,  then 

•sx  .  n    <=y  =  o  (  68  ) 

where  p.   denotes  the  Poisson's  ratio.   Then, 

«x  •  -  p   ey  <  69  ) 

Coefficients  of  the  row  and  column  related  to  nodal  parameter 
8u/3y  are  superimposed  by  the  coefficients,  which  are 
multiplied  by  -l//i.  of  the  row  and  column  related  to  nodal 
parameter   3u/8x.   The   coefficients   related   to   9u/9x   are 
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eliminated  . 

In  case  of  6y  =  O  on  a  certain  boundary  node,  then 

«r  '  ^%  =  0  (  70  ) 

or 

ey    ■    -A  «x  <    71) 

Superposition    of    coefficients    are    same    as    the    case    of    &yt    =0. 
In    case    of    yxy    =    0    on    a    certain    boundary    node, 

3u  8v 

*     =0  (    72    > 

8y  8x 

or 

3u  9v 

9y  c>x 

The   coefficient   of   the   row  and  column   related   to   nodal 
parameter  9u/9y  are  superimposed  by  the   coefficients,   which 
are  multiplied  by  -1,   of  the  row  and  column  related  to  Su f 8x 
and  the  coefficients  related  to  &u/ 9x  are  eliminated. 
SOLUTION  VECTOR: SKYLINE  PROGRAM 

This  program  finds  the  solution  vectors.  In  the  global 
stiffness  matrix,  a  number  of  zero  terms  must  still  be  stored 
and  manipulated  during  the  solution  process.  The  manipulation 
of  zero  terms  can  be  reduced  by  placing  further  limitations 
on  the  range  over  which  the  Gaussian  forward  elimination 
algorithms  are  applied.  To  impose  these  limitations,  it  is 
convenient  to  change  the  sequence  of  the  application  of  the 
algorithms.  Skylines  are  found  first  and  all  zeros  above  the 
skyline  of  the  stiffness  matrix  are  not  considered  xn 
calculation   process.   The  elements  of  stiffness   matrix   are 
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stored  in  a  one-dimensional  array,  starting  with  the  diagonal 
term  of  each  column.  The  Gaussian  algorithms  can  be  modified 
so  that  elimination  can  be  carried  out  directly  on  the  one- 
dimensional  stiffness  array.  This  program  was  developed  in 
the  Advanced  Structural  Analysis  II  course  (  32  )  and  was 
modified  slightly  in  this  study. 
REARRANGED  PARAMETERS IREPARAM  PROGRAM 

This  program  inserts  proper  values  of  nodal  parameters  to 
the  solution  vectors  found  by  the  use  of  skyline  technique. 
Zero  values  of  displacement  parameters  or  these  dependent 
nodal  parameters  were  eliminated  in  the  previous  program. 
After  restored  in  order,  they  can  be  used  for  finding  strain 
and  stresses  at  every  desired  point  of  the  structure. 
REPRESENTATION  OF  COMPLETE  DISPACEMENT  VECTOR :COMDISP  PROGRAM 

In  REDUC  program,  a  half  of  the  structure  including  nodes 

of   center   line  was   considered.   In   REPARAM   program,   the 

solution   of   a  half  of  the  structure  was   found.   Therefore. 

solution   of  another  half  of  the  structure  must  be   added   in 

actual    solution    vector.    Relationships    between    nodal 

i 
parameters  which  are  active  nodal  displacement  and  eliminated 

nodal   displacement   vectors.    By   this   relationships,   the 

solution   vectors   which  were  eliminated   during   calculating 

process  must  be  added. 

3-D  REPRESENTATION  OF  STRESS : STRESS  GRAPH  PROGRAM 

Nodal  displacement  vectors  were  found  by  COMDISP  program. 

This   program   calculates   the  stresses  and   strains   at   the 

points   of  interest  on  the  shear  walls  by  using  the   solution 
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vectors.   Three   dimensional  graphic  display  of  each   stress, 
such  as  normal  stress  S^,  Gy ,   and  shearing  stress  Txy ,   are 
plotted  by  this  program. 
2-D  REPRESENTATION  OF  DISPLACEMENT:  DISPLACE  PROGRAM 

This   program   will   show   the   two-dimensional    graphic 
display   of  the  displacement  of  shear  wall  by  using   solution 
vectors  calculated  in  COMDISP  program. 
3-D  REPRESENTATION  OF  STRAIN:  STRAIN  PROGRAM 

This   program   is  used  to  calculate  the   strains   at   the 
points   of   interest.   Three  dimensional  graphic  displays   of 
each  strain  are  plotted  . 
TRAJECTORIES  OF  STRESSES: TRAJECTORY  PROGRAM 

This  program  finds  the  trajectories  of  the  stresses  and 
gives  the  two-dimensional  graphic  display.  After  finding 
normal  stresses  <^x,  6y)  and  shearing  stress  Txy  ,  the 
magnitude  of  the  maximum  and  minimum  normal  stresses.  which 
are  principal  stresses,  is  found. 


Fig.  38  Mohr's  circle  for  principal  stresses 
Maximum  and  minimum  normal  stresses  are 

Sk  *  Sy     /  /«k  -  «y 


Si  =  OA  =  OC  *  CD  = 


*  Txy   <  74  ) 
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—   «x  ♦  6y     //«x-6yi 

S2  =  OB  =  OC  -  CD  =  -   /  1  *  i~y       <  75  > 


The   directions  of  principal  stresses  can  also  be   calculated 

from  the  Fig.  38  .   The  angle  DCA  is  the  double  angle  between 

the  stress  Si  and  the  x  axis  and  since  20  is  measured  from   D 

to  A  in  the  clockwise  direction,   the  direction  of  G\    must  be 

as  indicated  in  Fig.  38  . 

DE 
tan  20   =  --_---  (  76  i 

CE 

Regarding  the  sign  of  the  angle  0,  it  must  be  taken  negative 

in   this   case   since  it  is  measured  from  the  x  axis   in   the 

clockwise  direction.  Hence 

DE         2tx/ 

tan  2*  ■  -  -3-  =  -  (  77  ) 

CE       S„  -  Sy 

If   we  take  a  point  at  the  neutral  surface,   then  6*   becomes 

zero.   An   element  at  this  point  will  be  in  the  condition   of 

pure  shear.   The  directions  of  the  principal  stresses  will  be 

at  45*  to  the  x  and  y  axes.   The  trajectories  of  the  stresses 

are   two  systems  of  orthogonal  curves  whose  tangents  at   each 

point  are  in  the  directions  of  the  principal  stresses  at  this 

point.   All  these  curves  intersect  the  neutral  surface  at  45* 

and   have   horizontal  or  vertical  tangents   at   points   where 

shear  stress  Txy  is  zero. 
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CHAPTER  6 
CONCLUSIONS 

From  this  numerical  experimental  research,  the 
following   conclusions  can  be  reached: 

1.  Finite  element  method  provides  point-wise  value  of 
stresses  and  strains  in  the  shear  walls. 

2.  After  the  reduction  of  the  dimensionality  of  the 
matrix  equation,  the  Zenith  100  computer  can  handle 
relatively  large  problem  with  high  efficiency  and  accuracy. 
For  the  use  of  the  micro-computer,  the  reduction  of  matrix 
dimension  is  necessary  if  finite  element  method  is  used. 

3.  By  correctly  finding  the  principal  stresses  and 
principal  direction,  failure  of  shear  walls  can  be  reduced 
and  the  optimal  design  of  reinforcement  in  the  shear  walls 
can  be  obtained. 

4.  Linear  elastic  theory  is  not  adequate  for  predicting 
the  structural  response  of  the  structure  like  shear  wall.  It 
is  necessary  to  develop  the  non-linear  finite  element 
program  to  understand  the  non-linear  behavior  of  shear  walls. 

5.  The  stress  concentration  near  the  intersection  of 
the  wall  and  footing  can  not  be  accurately  calculated  without 
further  development. 

6.  The  experience  learned  through  the  development  has 
provided   the   capability  to  modify  and  improve  the   software 
available. 

7.  The  micro-computer  is  too  slow  for  carrying  out  the 
project  of  this  size. 


69 


8.  Buckling  analysis  of  shear  walls  is  not  included  in 
this  study.  For  thin  shear  walls,  buckling  analysis  should  be 
included. 

SUGGESTIONS  FOR  FURTHER  RESEARCH 

In  order  to  get  better  results  related  to  the  stress 
analysis  of  shear  walls,  the  following  suggestions  for 
further  research  are  recommended. 

1 .  The  number  of  elements  should  be  increased 
especially  near  the  points  of  singularities. 

2.  Need  the  development  of  nonlinear  finite  element 
program  for  the  effective  solution  of  shear  wall  problems. 

3.  Making  the  finite  element  method  as  a  part  of 
optimal  design  of  shear  wall. 

4.  Apply  the  modified  compression-field  theory  for 
reinforced  concrete  elements  subjected  to  shear <  33  )  for  the 
analysis  of  the  shear  walls  and  compare  the  results  to  the 
experimental  data. 

5.  The  additional  stress  due  to  the  interaction  of 
shearing  forces,  axial  forces,  and  bending  moments  to  be 
included  in  the  investigation. 
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10  REM 

20  REM  MAIN  PROGRAM 

30  REM  »***« mm **»*«« « 

40  REM  PROGRAM  TO  FIND   STIFFNESS  MATRIX  OF  THE   ELEMENTS   BY 

USING  FINITE  ELEMENT  METHOD 
50  DIM  NDC71.4) ,NN<4) ,X<71> ,Y<71> ,DC<3.3) ,B<3,24) ,DB<3.24> . 

SMC24.24) ,W#<3) .H#<3) .FX<4) ,FY<4) .PX<4) ,PY<4) , BDB (24.24) . 

SX(4> ,SY<4) ,SXX<4> ,SYX<4) .SXY<4) 
60  DIM  SYY<4) ,T<4) ,N<2.24) .TX<4) ,TY<4) .FFXC2) ,FFY<2> .FXX<2) . 

FYX(2).FXY(2) ,FYY(2) .GF (24) ,FF<24) . FC(24) ,XB <52. 14) . 

YBC52.14) .XEC52.14) .BAR<52) .THC4.9) 
70  DIM  DST(3)  ,DS<3.3)  .SSMC24.24)  ,CSM<24,24)  ,D<3,3>  ,YE(52,14) 
80  EC=2468410!:  REM   modulus   of  elasticity   of   cone.--  unit 

<N/sq.cm  ) 
90   ES=2E*07  :REM  modulus  of  elasticity  of   steel   bar--  unit 

(N/sq.cm  ) 
100  ESEC=ES-EC 

110  M=.2:R£M  ooisson's  ratio 

120  'AHB=2.01 :REM  area  of  horizontal  bar--  unit(sq.cm) 
130  AVB=l.l:REM  area  of  diagonal  bar--  unit(sa.cm) 
140  TH1=60 
150  TH2=10 
160  TH3=55 

170  REM  x,y  coordinate  of  each  node 
180  NM=71:REM  total  node  numbers 
190  NE=52:REM  no  of  total  elements 
200  FOR  1=1  TO  NM 

210  READ  X<I>  ' 

220  NEXT  I 

230  DATA      0.      0,   32.5,   32.5.   40. 0 ,  40. 0.  40 . 0 .  40 .0 , 

40.0.  40.0.  40.0 

240  DATA   50.0.  50.0.  50.0,  50.0.  50.0,  50.0.  50.0 

250  DATA   60.0.  60.0.  60.0.  60.0.  60.0.  60.0.  60.0 

260  DATA  120.0.120.0.120.0.120.0.120.0,120.0.120.0 

270  DATA  190.0.190.0.190.0.190.0.190.0.190.0.190.0 

280  DATA  260.0.260.0.260.0.260.0.260.0,260.0.260.0 

290  DATA  320.0.320.0.320.0.320.0.320.0.320.0.320.0 

300  DATA  330.0.330.0,330.0,330.0,330.0,330.0,330.0 

310  DATA  340.0.340.0.340.0.340.0.340.0.340.0.340.0 

320  DATA  347.5.347.5.380.0,380.0 

330  REM 

340  FOR  1=1  TO  NM 

350  READ  Y<I> 

360  NEXT  I 

370  DATA    0.0.  43.5,193.5.233.5 

380  DATA    0.0.  43.5.  53.5.  63.5.118.5.193.5.233.5 

390  DATA    0.0.  43.5.  53.5.  63.5.118.5,193.5.233.5 

400  DATA    0.0.  43.5.  53.5.  63.5.118.5.193.5.233.5 

410  DATA    0.0.  43.5.  53.5.  63.5.118.5,193.5.233.5 

420  DATA    0.0.  43.5.  53.5.  63.5.118.5.193.5.233.5 

430  DATA    0.0.  43.5.  53.5.  63.5.118.5.193.5.233.5 

440  DATA    0.0,  43.5.  53.5.  63.5.118.5.193.5.233.5 
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450  DATA    0.0.  43.5.  53.5.  63.5.118.5.193.5,233.5 

460  DATA    0.0.  43.5,  53.5,  63.5,113.5,193.5,233.5 

470  DATA  193.5.233.5,   0.0.  43.5 

480  REM 

490  REM  NODE  NO.  CORRESPONDING  TO  EACH  ELEMENT 

500  FOR  EL  =  1  TO  NE 

510  FOR  EN  =  1  TO  4 

520  READ  ND(EL.EN) 

530  NEXT  ENINEXT  EL 

540  DATA   1.5.6.2,5.12.13,6.12,19,20,13,19,26,27,20,26,33.34, 
27 

550  DATA  33.40.41.34.40.47.48,41,47,54.55.48,54,61,62,55,61.7 

0,71,62 
560  DATA  6,13.14,7,13,20,21,14,20,27,28,21,27,34,35.28,34.41, 

42.35 

570  DATA  41.48.49.42.48.55.56.49.55.62.63,56.7.14,15.8,14,21, 
22,15 

580  DATA  21.28,29,22,28,35,36,29,35,42,43,36.42,49,50,43,49.5 
6,57,50 

590  DATA  56.63,64,57,8,15,16,9,15.22,23,16,22,29,30,23,29,36, 
37,30 

600  DATA  36.43,44,37,43,50,51,44,50,57,58.51,57,64.65,58,9,16 

,17,10 

610  DATA  16,23,24,17.23,30,31,24.30,37,38,31,37,44,45,38.44,5 

1,52,45 
620  DATA  51,58,59,52,58,65,66,59,3,10,11,4,10,17,18,11,17,24, 

25,18 

630  DATA  24,31,32,25,31,38,39,32,38,45,46,39,45,52,53,46,52,5 

9,60,53 
640  DATA  59,66,67,60,66,68,69,67 
650  REM  bar  coordinate 
660  FOR  1=1  TO  NE 
670  READ  BAR<I) 
680  NEXT  I 
690  DATA  0,6.6,10,10,10.10,6,6,0,3,3,6,7,7,6.3.3,3,3,6.7,7,6, 

3,3.7,8,12.12.12,12 
700  DATA  8.7,10,10,12,14,14,12,10,10.0.6,6,8,10.10,8,6,6,0 
710  REM 

720  FOR  1=1  TO  NE 
730  BN=BAR(I> 
740  IF  BN  =  0  THEN  780 
750  FOR  J=l  TO  BN 
760  READ  XBd.J) 
770  NEXT  J 
780  NEXT  I 
790  DATA  43.43.43.43,43.43,50.50.50.50.57.57.60.60.70.90,110, 

60,60,70,90,110,120,120,130,150 
800  DATA  170,120,120.130,150,170.190,190,210.230,250.190,190, 

210 , 230 , 250 , 260 . 260 
310  DATA  270,290,310.260.260.270.290.310,320.320.320,320,323, 

323 , 330 , 330 , 330 , 330 , 337 . 337 , 43 . 43 , 43 . 50 , 57 , 57 
820  DATA  70,90,110.60,80,100.130,150,170,120.140.160,180,190, 

210,230,250,200 
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830  DATA  220.240.270.290.310.260.280.300.320.323.323.330,337. 
337 ,  43 .  43 ,  43 .  50 .  37  .  57  .  60 .  80  .  100 .  70 .  90 , 1 10  .  120 

840  DATA  140.160.180.130.150.170,200.220.240,190,210.230.250. 
260 , 280 , 300 , 270 . 290 

850  DATA  310.320.323.323,330,337,337,43,43,43,43.43.43,43.50. 

50 .  50 ,  50 ,  50 .  55 ,  57 ,  57 ,  60 ,  60 ,  60 .  70 ,  90  ,  1 10 ,  60 .  60 
860  DATA  60,75,95.115,120,120,120,130,150,170,120,120  120  135 

,155,175,190,190 

870  DATA  190,210.230.250,190,190,195,215,235,255,260,260,260, 
270 , 290 , 310 . 260 . 260 

880  DATA  260.275,295,315,320,320,320,320,320,320,323,323,330, 

330 ,  330  ,  330 ,  330 ,  337 ,  337 ,43,43,43,45 
890  DATA  43,43,43.43,43,43,50.50,50,50.50.50,50,50,57,57.60.6 

0,60,65,85,105,60,60,60,60,80 

900  DATA  100,120,120,120,125,145,165.185,120,120,120.120,140, 
160,180, 190,190, 190 

910  DATA  190,205,225,245,190,190,190.190,200,220,240,260,260, 
260 , 265 , 285 . 305 , 260 

920  DATA  260,260,260,280.300,320,320,320,325,320.320.320.320, 
323 , 323 , 330 , 330 , 330 , 330 , 330 

930  DATA  330.330,330,337,337,43,43,43.43,43,43,50,50.50,50.57 

, 57 , 60 , 60 , 80 , 100 , 60 , 63 , 83 , 103 , 120 . 120 
940  DATA  140,160,180,120,123,143.163,183,190,190.200.220,240, 

190,190,203,223,243 

950  DATA  260,260,280.300.260,263,283,303,320,320,320,323,323, 

323 , 330 , 330 , 330 , 330 
960  DATA  337,337 
970  FOR  1=1  TO  NE 
980  BN=BAR<I> 
990  IF  BN=0  THEN  1030 
1000  FOR  J  =  l  TO  BN 

1010  READ  YBCI.J) : YB< I , J) = YB< I , J) ♦ . 5 
1020  NEXT  J 
1030  NEXT  I 
1040  DATA  36,16.10,30,2.5,2.5,23,3.23,43,2.5.2.5,33,13,3.3,3 

13.33,43,43.43.33,13,3,3,3,13,33,43,43 
1050  DATA  43.23.3,3,3,3,23,43,43.43,43,33.13,3,3.3.13,33,43.4 

3,43,33,13.13,33,2.5.2.5,23 
1060  DATA  3.23,43,2.5,2.5,50,43,43,43,43,43,43.43,43,53.53,53 

.43,43,43,53,53,53,53,43,43.43,43,53.53 
1070  DATA  53.43,43.43,53.53.53.53,43.43.43.43.43,56.53.53.63. 

53,53,53,53.53,63,63,63,53,53,53,53,63,63 
1080  DATA  63,53,53.53,63,63,63.63,53,53,53,63.63.63.53.53.53. 

63,53.53,96.76,70,90,110,63,63,103.83 
1090  DATA  63,83,103,118,63,63.113,93,73,63,63,63,73,93,113,11 

8,118,118,113,93,73,63,63 
1100  DATA  63,73,93,113,118.118,118,103,83,63,63,63,63,83,100, 

118,118,118,118,113 
1110  DATA  93,73.63,63.63,73,93.113,118,118,118.113,93,73,73,9 

3,113,63,63,103,83,63,83 
1120  DATA  103.63,63.176,156,136,118,130,150,170,190,118.118,1 

83.163.143,123,123,143.163,183,118.118 
1130  DATA  173.153.133,118,118,118,133,153.173,193,193,193.173 
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.153,133.118,118 

1140  DATA  118,118,133.153.173,193,193,193.193.183.163.143,123 
,118,118,118,123 

1150  DATA  143,163.133,193,193,193.173,153,133,118,118.118.133 
,153,173,193,193 

1160  DATA  193,173,153,133,118,133,153,173,193.118,118,183,163 

,143,123,123,143,163,183,118,118 
1170  DATA  216,196,210,230.193,193,223,203,203,223,193,193,213 

,193,193,193.213,230,230,230,213 
1180  DATA  193,193,193,193,213,230.230,230,230,223,203,193,193 

,193,203,223,230 
1190  DATA  230,230,213,193,193,193,213,230,230,230,213,193,213 

,230,193,193,223,203,203 
1200  DATA  223,193,193 
1210  REM 

1220  FOR  1=1  TO  NE 
1230  BN=BAR<I) 
1240  IF  BN  =  0  THEN  1280 
1250  FOR  J=l  TO  BN 
1260  READ  XE(I.J> 
1270  NEXT  J 
1280  NEXT  1 

1290  DATA  50,50.50,50,43,43,60,60,60.60,57,57,70,90,110,120,12 

0,70,90,110,120,120,130,150 
1300  DATA  170,190,190,130,150,170,190,190,210.230,250,260,260 

210,230,250,260 

1310  DATA  260,270,290,310,320,320.270,290,310,320.320,330,330 
330 , 330 , 323 , 323 , 337 , 337 

1320  DATA  337,337,337,337,50,43,43,60,43,43,80,100,120,70,90.1 

10,140,160.180,130,150.170,190,200,220 
1330  DATA  240.260,210,230,250,280,300,320,270.290.310,330,323, 

323,337,337,337,50,43,43,60,57,57,70,90,110.80 
1340  DATA  100,120,130,150.170,190,140,160,180,210,230,250  200 

220 , 240 , 260 , 270 , 290 
1350  DATA  310,280,300.320,330,323,323,337,337,337,50,50,50,50, 

50,43,43,60,60,60,60,60,60,57,57,65,85,105,120 
1360  DATA  120,120,70,90,110,120,120,120,125,145,165,185,190  1 

90,130,150,170,190 

1370  DATA  190,190,205,225,245,260,260,260,210,230.250,260,260 
,260,265,285,305 

1380  DATA  320,320.320.270,290.310.320,320,320,325,330,330.330 
, 330 , 330 , 323 , 323 , 337 , 337 

1390  DATA  337,337.337,337,337,50,50,50,50,50,50,50,50,43.43,6 

0,60,60,60,55.60,60,60,57,57,80,100,120 
1400  DATA  120,120,120,75,95,115,120,120,120,140,160,180,190,1 

90,190,190,135,155 

1410  DATA  175,190,190,190,190,200,220,240,260,260,260,260,195 
,215,235,255,260 

1420  DATA  260,260,280.300.320.320.320,320.275.295,315,320,320 
,320,330,330,330 

1430  DATA  330.330.330,330.330,323.323.337,337,337,337,335,337 

, 337 , 337 , 337 , 337 , 50 , 50 . 50 , 50 , 43 , 43 , 57 , 60 
1440  DATA  60.60,57,57,77,97.117,120,80,100,120.120,137,157,17 
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7.190.190,140,160.180.190 
1450  DATA  190.197.217,237.257.260.200.220.240.260.260.277.297 
,317,320,280.300 

1460  DATA  320,320.330.330.330,330.323,323,337.337,337.337,337 

,337 
1470  FOR  1=1  TO  NE 
1480  BN=BAR<I> 
1490  IF  BN  =  0  THEN  1530 
1500  FOR  J  =  l  TO  BN 

1510  READ  YE<I.J>:YE<I.J)=YE<I.J)*.5 
1520  NEXT  J 
1530  NEXT  1 
1540  DATA  43.23,3.23,43.43.33.13.13.33,43,43.43.43.43.33.13,3 

,3,3,13.33,43,43,43,43,23,3.3.3 
1550  DATA  3,23,43,43,43,33.13.3,3,3.13,33,43.43,43,33,13.3.3, 

3,13,33,43.23,3,23,43,43 
1560  DATA  30,10,16,36.43,43.43.53,53,53.53,53,53,53,53.43.43. 

43,53,53,53,43,43,43,43.53,53,53,53 
1570  DATA  43,43,43,53,53.53,43,43,43,43,53,53,50,53,53.63.63. 

63,53,63.63,63,63,63,53,53,53,63,63,63,63 
1580  DATA  53,53.53,63,63,63,53,53.53,53,63.63,63.53,53,53,63, 

63,63,56,63.63,103,83.63,83.103.118.118 
1590  DATA  113.93,73,73,93,113,118.118,118,118,118,113,93,73,6 

3,63,63,73,93,113,118,118 
1600  DATA  118,118,103.83,63.63.63.63.83.100.118.118.118,113,9 

3,73,63,63,63,73,93 
1610  DATA  113.118,118,118,113.93,73.63,63,63,73,93.113,118,10 

3,83,63,83,103,118,118,110 
1620  DATA  90,70,76,96,118,118,183,163,143.123.123,143.163.183 

,193,193,193,173.153.133,118,133 
1630  DATA  153,173.193.193,193,193,193,173,153.133,118.118.118, 

133,153,173,193,193,193 
1640  DATA  183,163,143,123,118,118,118,123,143,163,183.193.193, 

193.193,173,153 
1650  DATA  133,118,118,118,118,133,153,173,193,193,193,173.153, 

133,118,118,118 
1660  DATA  133,153,173,183.163,143,123,123,143,163,183,193,193, 

190,170,150,130,118,136 
1670  DATA  156,176,193,193,223,203,203,223,230.230,230,213,193, 

213,230,230,230,230,230.213,193,193,193 
1680  DATA  213,230,230,230.223,203,193,193,193.203,223,230,230, 

230,230,213,193 
1690  DATA  193.193,193.213,230,230,230,213.193,193.193,213.223. 

203 . 203 , 223 , 230 , 230 , 230 
1700  DATA  210,196,216,230,230 
1710  FOR  1=1  TO  NE 
1720  BN=BAR<I> 
1730  FOR  J  =  l  TO  BN 

1740  PRINT  I.J.XB<I.J>.YB<I.J> ,XE<I.J) .YECI.J) 
1750  NEXT  J 
1760  PRINT 
1770  NEXT  I 
1780  FOR  1=1  TO  3 
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1790  READ  W#lI).H#<I> 

1800  NEXT  X 

1810  REM 

1820  DATA  -.774596669241483,  .555555555555556 

1830  DATA   .0  ,  .888888888888889 

1840  REM 

1850  DATA   .774596669241483,  .555555555555556 

I860  D1=EC/<1-M»M> 

1870  DC<1.1)=D1:DC<1,2)=D1»M 

1880  DC<2,1)=DC<1.2) : DC (2, 2) =D1 : DC (3,3) = (1-M) »Dl/2 

1890  OPEN  "R".#2."B:FCDATA',.4 

1900  FIELD  #2.4  AS  Cs 

1910  FOR  1=1  TO  (NE»24) 

1920  B  =  0 

1930  LSET  CS=MKSS<B> 

1940  PUT  #2,1 

1950  NEXT  I 

1960  CLOSE  #2 

1970  OPEN  "R",#1,"B:STDATA".4 

1980  FIELD  #1,4  AS  AS 

1990  FOR  1=1  TO  (NE»576) 

2000  A=0 

2010  LSET  AS  =  MKSsS<A) 

2020  PUT  #1,1 

2030  NEXT  I 

2040  CLOSE  #1 

2050  FOR  EL  =  1  TO  NE 

2060  FOR  11=1  TO  24 

2070  FOR  JJ=1  TO  24 

2080  SM ( 1 1 , J J ) =0 : CSM  < 1 1 , J J ) =0 : SSM  < I I . J J ) =0 

2090  NEXT  JJ 

2100  NEXT  II 

2110  ST=0 

2120  FOR  K=l  TO  4 

2130  NNCK)=0 

2140  NEXT  K 

2150  PRINT  "ELEMENT  NO":EL 

2160  FOR  K  =  l  TO  4 

2170  NN<K)=ND(EL.K> :PRINT  NN<K): 

2180  NEXT  K 

2190  FOR  1=1  TO  3 

2200  X#=0:HX#=0 

2210  X#=W#<I) :HX#=H#<I> 

2220  X=CX(NN<1)) tX(NN<2)))/2f(X(NN<2) )-X(NN<D) )»X#/2 

2230  V=X-X<NN<1>> : L=X ( NN < 2) > -X ( NN < 1 ) > 

2240  GOSUB  5000 

2250  FX ( 1 ) =F1 : FX < 2 ) =F2 : FX ( 3 ) =F3 : FX < 4 ) =F4 

2260  PX(1)=F5:PX<2)=F6:PX<3)=F7:PX(4)=F8 

2270  FOR  J=l  TO  3 

2280  Y#=0:HY#=0 

2290  Y#=W#(J) :HY#=H#<J) :HXY#=HX#»HY# 

2300  Y=<Y<NN(l)>*Y<NN<4>)>/2*<Y<NN<4>>-Y<NNCl)))*Y#/2 

2310  V=Y-Y<NN(1)):L=Y(NN<4))-Y(NN<1)) 
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2320  GOSUB  5000 

2330  FY<1)=F1:FY(2>=F2:FYC3>=F3:FY(4)=F4 

2340  PY<1)=F5:PY(2>=F6:PY<3>=F7:PY<4>=F8 

2350  GOSUB  5100 

2360  NEXT  J, I 

2370  A=<XCNN(2))-X<NN<l))>/2:B=<Y<NN<4>>-Y<NNCl>)>/2 

2380  AB=A»B 

2390  IF  EL<=10  THEN  AB=AB»TH1 

2400  IF  EL>=43  THEN  AB=AB«TH3 

2410  IF  EL>10  AND  EL<43  THEN  AB  =  AB»TH2 

2420  FOR  11=1  TO  24 

2430  FOR  JJ  =  1  TO  24 

2440  CSM(II,JJ)=SM<II,JJ)»AB 

2450  SMCII,JJ)=0 

2460  NEXT  JJ.II 

2470  REM 

2480  PRINT 

2490  ST=1 

2500  PRINT  "    STEEL  PART" 

2510  PRINT 

2520  BN=BAR(EL> 

2530  IF  BN=0  THEN  2940 

2540  FOR  EL1=1  TO  BN 

2550  PRINT  "      ELEMENT"; EL:"  BAR  N0";EL1 

2560  XB1=XB<EL.EL1> 

2570  XE1=XE<EL,EL1) 

2580  YB1=YB<EL.EL1> 

2590  YE1=YE<EL.EL1> 

2600  LENG=SQR< (XE1-XB1) ~2* < YE1 -YB1) ~25 

2610  HLEN=LENG/2 

2620  'IF  <YE1-YB1)=0  THEN  BT  =  2  ELSE  BT  =  1 

2630  SN=<YE1-YB1)/LENG 

2640  CS=<XE1-XB1)/LENG 

2650  DST<1)=CS~2 

2660  DST<2)=SN~2 

2670  DST<3)=SN»CS 

2680  FOR  11=1  TO  3 

2690  FOR  JJ=1  TO  3 

2700  DS(II.JJ)=DST<II)»DST<JJ) 

2710  NEXT  JJ.II 

2720  FOR  1=1  TO  3 

2730  X#=0:HX#=0 

2740  X#=W#<I) :HX#=H#<I> :HXY#*HX#»ESEC 

2750  X=<XEl*XBl)/2-KXEl-XBl)»X#/2 

2760  V=X-X(NN<1)):L=X<NN<2))-X(NN<1)) 

2770  GOSUB  5000 

2780  FX<1>=F1:FX<2)=F2:FX<3>=F3:FX<4)=F4 

2790  PX<1)=F5:PX<2)=F6:PX<3)=F7:PX<4)=F8 

2800  Y=<YEl+YBl)/2*<YEl-YBl)»X#/2 

2810  V=Y-Y<NN<1)):L=Y<NN<4))-Y<NN<1)> 

2820  GOSUB  5000 

2830  FY<1)=F1:FY<2>=F2:FY<3>=F3:FY<4>=F4 

2840  PY<1>=F5:PY<2>=F6:PY<3>=F7:PYC4)=F8 
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2850  GOSUB  5100 

2860  NEXT  I 

2870  A  =  AVB 

2880  AB=A»HLEN 

2890  FOR  11=1  TO  24 

2900  FOR  JJ  =  1  TO  24 

2910  SSM  < 1 1 , J J ) =SM  < 1 1 , JJ  > ♦ AB -SSM  < I I , JJ ) 

2920  SM<II,JJ)=0 

2930  NEXT  J J, II 

2940  NEXT  ELI 

2950  FOR  11  =  1  TO  24 

2960  FOR  JJ*1  TO  24 

2970  SSM  <  1 1 ,  J  J  >  =CSM  <  1 1 ,  J  J  >  *SSM  < II , JJ ) 

2980  NEXT  J J, II 

2990  OPEN  "R",#1,"B:STDATA",4 

3000  FIELD  #1,4  AS  AS 

3010  FOR  11  =  1  TO  24 

3020  FOR  JJ  =  1  TO  24 

3030  IJ=<II-1)«24+JJ-KEL-1)»24»24 

3040  LSET  AS  =  MKSi5<SSM(II,JJ>> 

3050  PUT  #1,1 J 

3060  NEXT  J J, II 

3070  CLOSE  #1 

3080  'PRINT  : 'PRINT 

3090  'PRINT  "ELEMENT  NO.=";EL 

3100  'PRINT 

3110  'FOR  11=1  TO  24 

3120  'FOR  JJ=1  TO  24 

3130  'PRINT  ■■SM<";II;JJ;")  =  ";SM(II,JJ);TAB(30);"SM<";JJ;II;") 

=";SM<JJ,II> 
3140  'NEXT  JJ 
3150  'PRINT 
3160  'NEXT  II 

3170  IF  EL=>11  AND  EL=<42  THEN  4980 
3180  FOR  1  =  1  TO  4 

3190  SX<I)=0:SY<I)=0:SXX(I)=0:SYX(I)=0:SXY<I)=0:SYY<I)=0:T<I) 

=0:TX<I)=O:TY<I)=0 
3200  NEXT  I 
3210  FOR  11=1  TO  24 
3220  FC<II>=0 
3230  NEXT  II 
3240  FOR  1=1  TO  4 
3250  READ  SX < I) ,SY ( I ) ,SXX < I ) ,SYX < I ) ,SXY < I ) ,SYY < I) ,T (I ) ,TX (I ) , 

TYCI) 
3260  NEXT  I 

3270  DATA  0,1940.44,0,10.2129,0,0,     O,  -8.3100,0 
3280  DATA  0,1531.93,0,10.2129,0,0,  297.42,-6.561,0 
3290  DATA  0,0.0,0,0,0,0,0,0 
3300  DATA  0,0,0,0.0,0,0,0,0 

3310  DATA  0,1531.93.0,10.2129,0,0,  297.42,-6.561,0 
3320  DATA  0,1429.80,0,10.2129,0,0,  360.84,-6.123,0 
3330  DATA  0,0,0,0,0,0,0,0,0 
3340  DATA  0,0,0,0,0,0,0,0,0 
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3350 

DATA 

3360 

DATA 

3370 

DATA 

3380 

DATA 

3390 

DATA 

3400 

DATA 

3410 

DATA 

3420 

DATA 

3430 

DATA 

3440 

DATA 

3450 

DATA 

3460 

DATA 

3470 

DATA 

3480 

DATA 

3490 

DATA 

3500 

DATA 

3510 

DATA 

3520 

DATA 

3530 

DATA 

3540 

DATA 

3550 

DATA 

0.1429.80,0,10.2129,0,0,  360.84,-6.123,0 

0,1327.67,0,10.2129,0,0,  419.89.-5.686,0 

0,0,0,0,0,0,0,0,0 

0,0.0,0,0,0,0,0,0 

0,1327.67,0,10.2129,0,0,  419.89,-5.686,0 

0,  714.90,0,10.2129,0,0,  682.32,-3.062,0 

0,0,0,0,0,0,0,0,0 

0,0,0,0,0,0,0,0,0 

0,  714.90,0,10.2129,0,0,  682.32,-3.062,0 

0,     0   ,0,10.2129,0,0,  789.47,      0,0 

0,0,0,0,0,0,0,0,0 

0,0,0,0,0,0,0,0,0 

0,     0   ,0,10.2129,0,0,  789.47,      0,0 

0,-714.90,0,10.2129,0,0,  682.32,  3.062,0 

0,0,0,0,0,0,0,0,0 

0,0,0,0,0,0,0,0,0 

0,-714.90,0,10.2129,0,0,  682.32,  3.062,0 

0,-1327.67,0,10.2129,0,0,  419.89,  5.686,0 

0,0,0,0.0,0,0,0,0 

0,0,0.0,0,0,0,0,0 

0,-1327.67.0,10.2129,0,0,  419.89,  5.686,0 
3560  DATA  0,-1429.8,0,10.2129,0,0,  360.84,  6.123,0 
3570  DATA  0,0,0,0,0,0,0,0,0 
3580  DATA  0,0,0,0,0,0,0,0,0 

3590  DATA  0,-1429.8,0,10.2129,0,0,  360.84,  6.123,0 
3600  DATA  0,-1531.93,0,10.2129,0,0,  297.42,  6.561,0 
3610  DATA  0,0,0,0,0,0,0,0,0 
3620  DATA  0,0,0,0,0,0,0.0,0 

3630  DATA  0,-1531.93,0,10.2129,0,0,  297.42.  6.561,0 
3640  DATA  0,-1940.44,0,10.2129,0,0,       0,  8.31  ,0 
3650  DATA  0,0,0,0,0,0,0,0,0 
3660  DATA  0,0,0,0,0,0,0,0,0 
3670  DATA  0,0,0,0,0,0,0,0,0 
3680  DATA  0,0,0,0,0,0,0,0,0 
3690  DATA  0,0,0,0,0,0,  88.54,  11.518,0 
3700  DATA  0,0,0,0,0,0,      O,  12.094,0 
3710  DATA  0,0,0,0,0,0,0,0,0 
3720  DATA  0,0,0,0,0,0,0,0,0 
3730  DATA  0,0,0,0,0,0,199.88,  10.750,0 
3740  DATA  0,0,0,0,0,0,  88.54,  11.518,0 
3750  DATA  0,0,0,0,0,0,0,0,0 
3760  DATA  0,0,0,0,0,0,0,0,0 
3770  DATA  0,0,0,0,0,0,303.54,   9.982,0 
3780  DATA  0,0,0,0,0,0,199.88,  10.750,0 
3790  DATA  0,0,0,0.0,0,0,0,0 
3800  DATA  0,0,0,0,0,0,0.0,0 
3810  DATA  0.0.0,0,0,0,764.26,   5.375,0 
3820  DATA  0,0,0,0,0,0,303.54,   9.982,0 
3830  DATA  0,0.0,0,0.0,0.0,0 
3840  DATA  0,0,0,0.0,0,0,0,0 
3850  DATA  0,0,0,0,0,0,952.38,       0,0 
3860  DATA  0,0.0,0,0,0,764.26,   5.375,0 
3870  DATA  0,0,0,0,0,0,0,0,0 
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3880 
3890 
3900 
3910 
3920 
3930 
3940 
3950 
3960 
3970 
3980 
3990 
4000 
4010 
4020 
4030 
4040 
4050 
4060 
4070 
4080 
4090 
4100 
4110 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4380 
4390 
4400 


DATA 
DATA 


DATA 
DATA 
DATA 
DATA 


.375,0 
0,0 


-9.982,0 
-5.375,0 


-10.750,0 
-9.982,0 


0,-12.094,0 
88.54,-11.518,0 


DATA  0,0.0,0,0,0,0,0,0 
DATA  0,0,0.0,0,0,764.26, 
DATA  0,0,0.0,0,0,952.38, 
DATA  0,0,0,0,0,0,0,0,0 
DATA  0,0,0,0,0,0,0,0,0 

0,0,0,0,0,0,303.54, 

0,0,0,0,0,0.764.26, 
DATA  0,0,0,0,0,0,0,0,0 
DATA  0,0,0,0.0,0,0,0,0 

0,0,0,0,0,0,199.88, 

0,0,0,0,0,0,303.54, 

0,0,0,0,0,0,0,0,0 

0,0.0,0,0,0,0.0,0 
DATA  0,0,0.0,0,0,  88.54,-11.518, 
DATA  0,0,0,0,0,0,199.88,-10.750, 
DATA  0,0,0,0,0,0,0,0,0 
DATA  0,0,0,0,0,0,0,0,0 
DATA  0,0,0,0,0,0, 
DATA  0,0.0,0,0,0, 
FOR  SIDE=1  TO  4 
FOR  1=1  TO  24 
GF<I)=0:FF(I)»0 
NEXT  I 

I1=SIDE:I2=SIDE*1 
IF  I2>4  THEN  12=1 
L1=X(NN<I2))-X<NN(I1>) 
L2=Y<NN<I2) >-Y<NN<Il) ) 
LN=SQR(L1»L1*L2«L2> 
UX=L2/LN:UY=-L1/LN 
FOR  M=l  TO  2 

IF  M=l  THEN  Z  =  I1  ELSE  Z=I2 
FFX<M)=UX»SX<Z>*UY»T(Z> 
FFY<M)=UY«SY<Z)*UX»T(Z) 
FXX(M)=UX«SXX(Z)*UY»TX(Z) 
FYXCM)=UY*SYX<Z)*UX«TX<Z) 
FXY(M)=UX«SXY<Z)«-UY»TY<Z> 
FYY<M)=UY»SYY<Z>-UX»TY<Z> 
NEXT  M 

IF  SIDE=2  OR  SIDE=4  THEN  4340 
Y#  =  -l 

IF  SIDE=3  THEN  Y#=l 

Y=<Y<NN<l))+Y(NN<4)))/2*(Y(NN<4))-Y(NN(l)))»Y#/2 
V=Y-Y<NN<1)> :L=Y<NN(4))-Y<NN<1)> 
GOSUB  5000 

FY(l)=Fl:FY<2)=F2:FY(3)=F3:FY<4)=F4 
GOTO  4400 
X#  =  l 

IF  SIDE=4  THEN  X#=-l 

X=(X<NN(l))»X(NN(2)>>/2f (X < NN (2) ) -X (NN ( 1 ) ) ) «X#/2 
V=X-X<NN<1)> :L=X<NN<2))-X<NN<1>) 
GOSUB  5000 

FX<l)=Fl:FX(2)=F2:FX<3)=F3:FX<4)=F4 
FOR  GQ=1  TO  3 
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4410  IF  SIDE=2  OR  SIDE=4  THEN  4490 
4420  X#  =  0 
4430  X#=W#<GQ> 

4440  X=<X<NN(l))*X(NN<2)))/2-<X<NN(2))-X(NN(l>)).X#/2 
4450  V=X-X<NN<1>)  :L  =  X  (NN  (2)  )  -X  <NN<1>  ) 
4460  GOSUB  5000 

4470  FX<1)=F1:FX<2)=F2:FX<3>=F3:FX<4>=F4 
4480  GOTO  4550 
4490  Y#  =  0 
4500  Y#=W#(GQ) 

4510  Y=(Y<NN(l))*Y(NN<4)))/'2*<Y<NN<4))-Y<NN(l))).Y#/2 
4520  V=Y-Y<NN<1>):L=Y<NN<4)>-Y<NN<1>) 
4530  GOSUB  5000 

4540  FYC1)=F1:FY<2)=F2:FY<3>=F3:FY<4)=F4 

4550  FOR  11=1  TO  4 

4560  C=<II-1>»6 

4570  IF  11=1  OR  11=4  THEN  K=l  ELSE  K=2 

4580  IF  II<  =  2  THEN  KK  =  1  ELSE  KK  =  2 

4590  N(1,1*C)=FX<K>*FY(KK> 

4600  N(1.3+C)=FX<K+2)»FY(KK) 

4610  N(l,5tC)=FX<K)«FY(KK*2) 

4620  N<2.2*C)=N<1,H-C> 

4630  N(2,4*C)=N<1,3*C> 

4640  N<2.6*C)=N<1,5+C) 

4650  NEXT  II 

4660  PX=0:PY=0 

4670  FOR  M=l  TO  2 

4680  IF  M=l  THEN  PT=<I1-1>.6  ELSE  PT=<I2-1>»6 

4690  PX  =  FFX<M)»N<l,PTi-l)*FXX<M)«N(l,PTt3)*FXY<M).N<l.PT*5)»PX 

4700  PY  =  FFY<M>.N<2.PT  +  2>*FYX<M).N<2,PT*4>*FYY<M>.N<2,PT<-6>*PY 

4710  NEXT  N 

4720  FOR  1=1  TO  24 

4730  GF<I)=N<1.I).PX*N(2.I)«PY 

4740  NEXT  I 

4750  FOR  1=1  TO  24 

4760  FF<I)=FF<I)*GF<I).H#<GQ> 

4770  NEXT  I 

4780  NEXT  GQ 

4790  FOR  1  =  1  TO  24 

4800  IF   SIDB-1  OR  SIDE=3   THEN  L= (X <NN <2> ) -X <NN < 1 > > > /2   ELSE 

L=(Y(NN<4))-Y<NN<l))>/2 
4810  FCCI)=FF(I>«L*FC(I> 
4820  NEXT  I 
4830  NEXT  SIDE 
4840  PRINT  : PRINT 

4850  PRINT  "FORCE  VECTOR  OF  ELEMENT  NO.=":EL 
4860  PRINT  : PRINT 
4870  FOR  1  =  1  TO  24 
4880  PRINT  "FC(":I:")=":FC<IJ 
4890  NEXT  I 

4900  OPEN  •■R",#2."B:FCDATA".4 
4910  FIELD  #2.4  AS  CS 
4920  FOR  11=1  TO  24 
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4930  PT=<EL-1)«24*II 

4940  LSET  CS  =  MKS!3  (  FC  i  II)  ) 

4950  PUT  #2.PT 

4960  NEXT  II 

4970  CLOSE  #2 

4980  NEXT  EL 

4990  END 

5000  LV=L-V 

5010  F1  =  LV/L*V»LV«<L-2»V)/L'S3 

5020  F2  =  V/L  +  V*LV«(2»V-L)/'L~3 

5030  F3=V»LV"2/L"2 

5040  F4=-LV«V~2/L"2 

5050  F5  = (-6»V'L»6«V"2)/L"3 

5060  F6=<6»V«L-6»V"2)/L"-3 

5070  F7=<L~2-4*L»V*3»V~2)/L~2 

5080  F8=<-2»L«V*3«V'~2)/L"2 

5090  RETURN 

5100  FOR  11=1  TO  4 

5110  C=(II-1>»6 

5120  IF  11=1  OR  11=4  THEN  K=l  ELSE  K=2 

5130  IF  II  <  =  2  THEN  KK  =  1  ELSE  KK  =  2 

5140  B(1,1»C)=PX<K>»FY<KK> : B (2 , 2»C) =FX <K> »PY <KK> 

5150  B<1.3*C)=PX<K+2)»FY<KK) :B (2.4-C) =FX(K*2) »PY (KK) 

5160  B(1,5*C)=PX<K)»FY(KK»2> :B (2 ,6*C> =FX <K) »PY (KK*2) 

5170  B<3.1*C)=B<2.2»C):B<3.2tC)=B<l.l*C) 

5180  B<3.3*C)=B(2,4*C) :B<3.4*C) =B < 1 , 3*C) 

5190  B<3.5*C)=B<2.6*C> : B<3 ,6*C> =B C 1 . 5+C) 

5200  NEXT  II 

5210  FOR  11=1  TO  3 

5220  FOR  JJ  =  1  TO  3 

5230  IF  ST=0  THEN  5240  ELSE  5260 

5240  D<II.JJ)=DC(II,JJ> 

5250  GOTO  5270 

5260  D<II,JJ)=DS<II,JJ) 

5270  NEXT  JJ 

5280  NEXT  II 

5290  FOR  11=1  TO  3 

5300  FOR  JJ=1  TO  24 

5310  DB<II.JJ)=0 

5320  FOR  K=l  TO  3 

5330  DB  < 1 1 . J J ) =DB  < 1 1 , J J  >  *D  < 1 1 , K >  »B  <K , J J ) 

5340  NEXT  K.JJ.II 

5350  FOR  11=1  TO  24 

5360  FOR  JJ=1  TO  24 

5370  BDB<II»JJ)=0 

5380  FOR  K=l  TO  3 

5390  BDB(II,JJ)=BDB<II.JJ)-B<K,II)*DB<K,JJ) 

5400  NEXT  K.J J. II 

5410  FOR  11=1  TO  24 

5420  FOR  JJ=1  TO  24 

5430  SM<II,JJ)=SM<I1,JJ>-BDB(II,JJ)*HXY# 

5440  NEXT  JJ.II 

5450  RETURN 
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10  REM  PROGRAM  NAME  :  ASSEM.  PROGRAM 

20  REM  ............................. 

30  REM   Program  of  assembling  the  element   stiffness   matrix 

stored  in  the  random 
40  REM  access  file  "STDATA"  and  of  converting  to  the   global 

stiffness   matrix   which   will  be  stored  in  a   new   random 

access  file  "STDATAN" 

so  b=o:s=o:z=o:zz=o:zi=o:zzi=o:bb=o:5S=o:bk=o:bs=q 

70  REM 

80  DIM  ND<4.4) 

90  FOR  1=1  TO  4 

100  FOR  J  =  l  TO  4 

110  READ  ND(I.J) 

120  NEXT  J.I 

130  DATA   1,  4,  5.  2.   4,  7.  8.  5.    2.  5,  6.  3.    5.  a.  9.  6 

140  REM 

150  REM 

160  OPEN  "R".#3,,,B:STDATAN",4 

170  FIELD  #3,4  AS  BS 

180  FOR  1=1  TO  1500 

190  B  =  0 

200  LSET  B«MMKS«KB) 

210  PUT  #3,1 

220  NEXT  I 

230  CLOSE  #3 

240  B=0 

250  OPEN  "R",#1,"B:STDATA",4 

260  FIELD  #1,4  AS  AS 

270  OPEN  ■•R",#2."B:FCDATA".4 

280  FIELD  #2,4  AS  CS 

290  OPEN  "R",#3."B:STDATAN".4 

300  FIELD  #3,4  AS  BS 

310  OPEN  "R",#4."B:FCDATAN",4 

320  FIELD  #4,4  AS  DS 


330 

FOR  1=1 

TO 

4 

340 

FOR  J=l 

TO 

4 

350 

REM 

360 

REM 

370 

REM 

380 

NJ=ND(I, 

,J>  ; 

IPRINT  NJ 

390  NR=(NJ-1)«6 

400  FOR  K=l  TO  6 

410   NK=NR*K:REM  NK  is  the  row  no.   in  the   global   stiffness 

matrix 
420  ELM=<NK-1) «NK/2:REM  Total  nos  of  lower  triangular  element 

before  node  no.  considered. 
440   5R=(J-1) »6*K:REM   SR   is   the  row   no.   in   the   element 

stiffness  matrix 
450  FOR  L=l  TO  4 
460  NL=ND(I,L) 
470  CN=(NL-1) »6 
480  FOR  M=l  TO  6 
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490  KC  =  CN<-M:REM  KC  is  the  column  no.   in  the  global  stiffness 

matrix 
500  SC=(L-1) «6»« 
510  IF  KONK  THEN  700 

520  REM  Z  is  the  position  of  the  element  "STDATA" 
522  REM  global  stiffness  matrix 
530  Z=<5R-1)»24*SC+(I-1)  »24»24 
540   REM  ZZ  is  the  corresponding  position  of  the   element   in 

"STDATAN" 
550  Z2=ELM*KC 
560  GET  #1,Z 
570  S=CVSCAS) 
530  GET  #3.ZZ 
590  B  =  CVS(BSO  :BK  =  B 
600  B=B*S:BS=B 
610  LSET  BS=MKSS(B> 
620  PRINT"Z=":Z:TAB<10) j"ZZ=":ZZ! :TAB<21) : "S=" :S:TAB<38) :"B=" 

;BK:TAB<56) ;"BS=";BS 
630  REM 
640   REM   Z   is  the  position  of  the  element   in   the   element 

stiffness   matrix  stored  530  REM  in  a  1-dimensional  array. 

and  zz  is  its  corresponding  position  in  the 
660  REM 
670  PUT  #3.ZZ 

680  b=o:s=o:z=o:zz=o:bs=o:bk=o 

690  'PRINT  "B  = " : B : "S= " : S : "KC= " ; KC : "NK= " : NK : "ZZ= " : ZZ : "Z= " ; Z 

700  NEXT  M 

710  NEXT  L 

720  Z1=<I-1)»24*<J-1)»6*K 

730  ZZ1=NK 

740  GET  #2.Z1 

750  SS=CVS(CS) 

760  GET  #4,ZZ1 

770  BB=CVS<DS) 

780  BB=BB+SS 

790  LSET  DS=MKSS(BB> 

S00  PUT  #4.ZZ1 

810  BB=0:ss=0:Zl=o:zzi=0 

820  NEXT  K 

830  NEXT  J 

840  NEXT  I 

850  CLOSE  #1 

860  CLOSE  #2 

870  CLOSE  #3 

880  CLOSE  #4 

890  END 
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10  REM  PROGRAM  NAME:   REDUC.  PROGRAM 

20  REM  «••».*.»»••..,..»«..„„„„„,,,„ 

30  REM   Program  of  assembling  the  element   stiffness   matrix 

stored  in  the  random  access  file  "STDATA"  and  of  converting 
40  REM  to  the  reduced  global  stiffness  matrix  which  will   be" 

stored  in  a  new  random  access  file  "STDATAN" 
60  B=0:s=o:Z=o:ZZ=o:zi=0:ZZ1=0:bb=0:SS=0:bk=0:BS=o 
70  REM 

80  DIM  ND(54.4> ,C0MP(426.2) 
90  TN=71:REM  TN=T0TAL  NODE  NUMBER 
100  N=TN»6 
110  TEL=52:REM  TEL=T0TAL  ELEMENT  NO. 

112  LN=39:REM  LAST  NODE  NUMBER  OF  THE  REDUCED  STRUCTURE 

113  L1=LN«6:L2=L1*1 
120  FOR  1=1  TO  TEL 
130  FOR  J  =  l  TO  4 
140  READ  ND(I.J) 
150  NEXT  J.I 

160  DATA   1.  5.  6.  2.   5.12.13.  6.   12.19.20.13.   19.26.27.20. 

26.33.34,27 
170  DATA  33.40.41.34,  40,47,48.41,   47,54,55.48.   54.61,62.55. 

61,70.71.62 
180  DATA   6.13.14.  7.  13.20.21.14.   20.27.28.21.   27.34.35,28. 

34,41,42,35 
190  DATA  41,48,49,42,  48.55,56,49,   55,62,63,56,    7,14,15.  8, 

14,21,22.15 
200  DATA  21,28.29.22,  28.35.36,29,   35,42,43,36,   42,49.50,43. 

49,56,57,50 
210  DATA  56,63,64,57,   8,15,16,  9,   15,22,23,16,   22,29,30,23. 

29,36,37,30 
220  DATA  36,43.44,37,  43,50,51,44,   50,57,58,51.   57.64.65,58, 

9,16,17,10 
230  DATA  16,23,24.17,  23,30,31,24,   30,37,38,31,   37,44,45.33. 

44.51,52,45 
240  DATA  51,58,59.52,  58,65.66,59.    3,10.11.  4,   10,17,16.11. 

17,24,25,18 
250  DATA  24,31,32.25.  31.38.39,32,   38,45,46,39,   45,52,53,46, 

52,59,60,53 
260  DATA  59.66,67,60.   66.68.69.67 
270  IJ  =  0 

280  FOR  1=1  TO  LI 
290  C0MP<I.1)=I 
300  NEXT  I 
310  FOR  I=L2  TO  N 
320  READ  COMP(I.l) 
330  NEXT  I 

340  DATA  151,152,153,154,155,156,157,158,159,160,161,162.163, 

164,165,166,167,168 
342  DATA  169,170,171,172,173,174.175.176.177.178.179.180,181, 

182,183,184,185,186 
344  DATA  187.188,189.190,191,192.109,110.111,112.113.114,115, 

116,117,118.119,120 
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346  DATA  121,122.123,124,125.126,127.128,129,130,131,132.133. 

134.135.136.137.138 
348  DATA  139.140.141.142.143.144.145.146.147.148.149.150 

403  DATA  67.68.69,70,71,72,73,74.75,76.77,78.79,80,81.82.83.8 
4.85.86.37,88.89.90 

404  DATA  91.92.93.94.95.96.97.98.99.100,101,102,103,104.105.1 
06,107,108 

410  DATA  25,26,27.28,29.30.31,32,33.34.35,36,37,38,39,40,41.4 
2.43.44.45.46,47 

420  DATA  48,49,50,51,52,53.54.55,56.57.58,59.60.61.62.63.64.6 

5,66 
430  DATA  13,14,15.16,17,18,19,20,21.22,23.24 
440  DATA  1,2,3.4.5.6,7.8,9,10,11,12 

450  FOR  1  =  1  TO  LI 

460  C0MP<I,2)=1 

470  NEXT  I 

480  FOR  I=L2  TO  N 

490  READ  C0MP(I.2> 

500  NEXT  I 

510  DATA  1,-1,-1,1,1,-1.   1,-1.-1.1,1,-1,  1,-1,-1,1,1,-1.  1,- 

1,-1.1,1,-1 
520  DATA  1,-1,-1.1.1,-1,   1,-1,-1.1,1,-1,  1,-1,-1,1,1,-1,  1,- 

1,-1,1,1,-1 
530  DATA  1,-1,-1,1,1,-1,   1,-1,-1,1,1,-1,  1,-1,-1,1,1,-1,  1,- 

1,-1.1.1,-1 
540  DATA  1,-1,-1,1,1,-1,   1,-1,-1,1,1,-1.  1.-1,-1.1,1,-1,  1,- 

1.-1.1,1.-1 
550  DATA  1,-1,-1,1,1,-1,   1,-1,-1,1,1,-1,  1,-1,-1,1,1,-1,  1.- 

1,-1,1,1,-1 
560  DATA  1.-1.-1,1,1,-1.   1,-1,-1.1,1,-1,  1,-1,-1.1,1,-1.  1,- 

1.-1.1,1.-1 
570  DATA  1,-1.-1.1.1.-1.   1.-1.-1,1,1,-1,  1,-1,-1,1,1,-1,  1,- 

1.-1.1,1,-1 
580  DATA  1.-1,-1.1,1,-1,   1,-1,-1,1.1,-1.  1.-1,-1.1,1,-1,  1,- 

1,-1.1,1,-1 
590  OPEN  "R".#3,"B:STDATAN",4 
600  FIELD  #3,4  AS  BS 
610  FOR  1=1  TO  <Ll»L2/2) 
620  B=0 

630  LSET  BS=MKSS<B) 
640  PUT  #3,1 
650  NEXT  I 
660  CLOSE  #3 

670  OPEN  "R",#4,"B:FCDATAN",4 
680  FIELD  #4,4  AS  DS 
690  FOR  1=1  TO  LI 
700  B=0 

710  LSET  Df3  =  MKSS(B) 
720  PUT  #4.1 
730  NEXT  I 
740  CLOSE  #4 

750  OPEN  "R" .#1 ,"B:STDATA",4 
760  FIELD  #1,4  AS  AS 
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770  OPEN  "R" ,#2, "B : FCDATA" . 4 

780  FIELD  #2,4  AS  CS 

790  OPEN  ,'R".#3."B:STDATAN",4 

800  FIELD  #3,4  AS  BS 

310  OPEN  "R",#4,"B:FCDATAN",4 

820  FIELD  #4,4  AS  DS 

830  FOR  ELMT=1  TO  TEL 

835  PRINT  "ELEMENT  NO":ELMT 

840  FOR  RCRN=1  TO  4 

850  REM 

860  REM 

870  REM 

880  NP=NDCELMT,RCRN) :PRINT  NP 

890  NR=(NP-1>»6 

900  FOR  K=l  TO  6 

910  NNK=NR*K:REM  NK  is  the  row  no.   in  the  global   stiffness 

matrix 
920  NK=COMP(NNK.l) 
930  FC=C0MP<NNK,2> 
940  ELM=<NK-1)«NK/2:REM  Total  nos  of  lower  triangular  element 

before  node  no.  considered. 
960   SR=(RCRN-1)»6*K:REM   SR  is  the  row  no.   in   the   element 

stiffness  matrix 
970  FOR  L=l  TO  4 
980  NL=ND<ELMT,L) 
990  CN=(NL-1)»6 
1000  FOR  M=l  TO  6 
1010  KKC  =  CN*M:REM  KC  is  the  column  no  in  the  global  stiffness 

matrix 
1020  KC=C0MPCKKC,1>  :FFC  =  C0MP  CKKC2) 
1030  SC=<L-1) »6+M 
1040  IF  KONK  THEN  1230 

1050  REM  2  is  the  position  of  the  element  "STDATA" 
1060  Z=<SR-l)*24-SO<ELMT-l)»24»24 
1070   REM  ZZ  is  the  corresDonding  position  of  the  element   in 

"STDATAN" 
1080  GET  #1,Z 
1090  S  =  CVS<As3) 
1100  ZZ=ELM*KC 
1110  GET  #3,ZZ 
1120  B  =  CVS(BS!>  :BK»B 
1130  B=B*S»FC«FFC 
1140  LSET  BS=MKSS<B> 
1150  'PRINT"Z=":Z:TAB(9> ;"ZZ=" ;ZZ;TAB(19) ; "S»" ;S; TAB<35> ;"B=" 

;BK;TAB<53)  ;  "BS  =  "  ;  B;  TAB  < 67)  ;  "FC""  ;  FC;TAB<73>  ;  "FFC="  ;  FFC 
1160  REM 
1170   REM   Z  is  the  position  of  the  element   in   the   element 

stiffness  matrix  stored  530  REM  in  a  1-dlmensional  array 

,  and  zz  is  its  corresponding  position  in  the 
1180  REM  global  stiffness  matrix 
1190  REM 
1200  PUT  #3,ZZ 
1210  B=o:S=o:z=o:zz=0:bk=o 
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1220    'PRINT    ••B  =  ";B;"S  =  ";5;"KC  =  ";KC;"NK  =  ";NK;"22=":ZZ;"2=";Z 

1230  NEXT  M 

1240  NEXT  L 

1250  Z1=<ELMT-1>»24+<RCRN-1)»6*K 

1260  ZZ1=NK 

1270  GET  #2,Z1 

1280    SS  =  CVS<Cs3) 

1290  GET  #4,ZZ1 

1300  BB=CVS<DS) 

1310  BB  =  BB«-SS»FC 

1320  LSET  DS=MKSS<BBJ 

1330  PUT  #4,ZZ1 

1340  BB=0:SS=0:Z1=0:ZZ1=0 

1350  NEXT  K 

1360  NEXT  RCRN 

1370  NEXT  ELMT 

1380  CLOSE  #1 

1390  CLOSE  #2 

1400  CLOSE  #3 

1410  CLOSE  #4 

1420  END 
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10  REM  PROGRAM  NAME  :  MODIF  PROGRAM 

20  REM  ............................ 

30  REM  In  this  program  boundary  condition  is  considered.   In 

case  of  fixed  end, 
40  REM  set  rows  and  columns  of  the  global  stiffness  matrix  tc 

be  zero  except 
SO  REM  making  diagonal  elements  be  zero. 
60  REM 

70  DIM  ZER0C6) 
80  FOR  1=1  TO  6 
90  READ  ZERO (I) 
100  NEXT  I 
110  REM 

120  DATA  1.2,19,20,37,38 
130  FOR  1=1  TO  38 
140  P=0 
150  REM 
160  REM 

170  FOR  J=l  TO  6 
180  IF  P=l  THEN  200 
190  IF  I=2ER0(J)  THEN  P=l 
200  NEXT  J 
210  IF  P=0  THEN  330 
220  PRINT  I.P 
230  FOR  K=l  TO  I 
240  S=0 

250  IF  I=K  THEN  S=l 
260  XK-((X-1>»I/2)*K 
270  OPEN  "R",#3,"B:STDATAN",4 
280  FIELD  #3,4  AS  AS 
290  LSET  AS=MKSS(S) 
300  PUT  #3,  IK 
310  CLOSE  #3 
320  NEXT  K 
330  NEXT  I 
340  FOR  1=1  TO  38 
350  P=0 
360  REM 

370  REM  The  following   makes  the  column  elements  be  zero 
380  REM 

390  FOR  J  =  l  TO  6 
400  IF  P=l  THEN  420 
410  IF  I=ZER0<J)  THEN  P=l 
420  NEXT  J 
430  IF  P=0  THEN  540 
440  PRINT  I.P 
450  FOR  K=I  TO  53 
460  S=0 

470  IK=K»<K+l)/2+I 
480  OPEN  "R",#3,"B:STDATAN",4 
490  FIELD  #3,4  AS  AS 
500  LSET  AS=MKSS(S> 
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510  PUT  #3, IK 

520  CLOSE  #3 

530  NEXT  K 

540  NEXT  I 

550  FOR  1=1  TO  6 

560  OPEN  "R",#4,"B:FCDATAN".4 

570  FIELD  #4,4  AS  BsS 

580  S«0 

590  K=ZER0<I> 

600  PRINT   K 

610  LSET  BS=MKSS<S> 

620  PUT  #4,K 

630  CLOSE  #4 

640  NEXT  I 

650  END 
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10  REM 

20  REM  PROGRAM  NAME:  ENBOUN 

30  REM  *»**»•*»•***««»*»**•* 

40  REM   This  program  considers  matrix  superposition   of   the 

known  nodai  parameters 
50  REM  and  deletes  the  rows  and  columns  which  have  zero  nodal 

parameters,  and  finally  gives  the  reduced  lower  triangular 

matrix . 

70  DIM  ND(50> ,ACTND(50,3> .ZER0(50.9> ,CBN(3> .WKK300) ,WK2(300> 

. WK3(300> .WK4C300) 
80    DIM    ZER0K50) 
90    INPUT'POISSON'S    RATI0":M 
100    INPUT'TOTAL    NODE    N0":TN 
110    TL=TN»6 

120  INPUT"H0W  MANY  NODES  ARE  CONSIDERED" ; HN 
130  FOR  1  =  1  TO  HN 
140  PRINT  I 

150  INPUT"N0DE  NUMBER" ; ND( I > 
160  PT=(ND(I>-1) »6 

170  INPUT"        SIGMA  X  =  0  (yes type  1   else   type  0)  "  ;  Z 

ISO  IF  Z  =  l  THEN  ISO  ELSE  200 

190  ACTND<I.1)=PT*6:ZER0<I.1)=PT*3:G'0T0  210 

200  ACTND(I.1>=0 

210  INPUT"        TAU   XY  =  0  (yes type  1   else   type  0)  " : Z 

220  IF  Z  =  l  THEN  230  ELSE  240 

230  ACTND(I,2)=PT*4:ZER0(I,2)=PT*5:G0T0  250 

240  ACTND (I. 2) =0 

250  INPUT"        SIGMA  Y=0  (yea type  1   else   type  0)  ":Z 

260  IF  Z=l  THEN  270  ELSE  280 

270  ACTND(I,3)=PT*3:ZER0(I,3)=PT*6:G0T0  290 

280  ACTND(I,3)=0 

290   INPUT"         U  =  0        (yes type  1   else    type   0) 

":Z 

300  IF  Z  =  l  THEN  ZERO ( I ,4) =PT*1  ELSE  ZER0(I.4>=0 

310   INPUT"        V=0        (yes type  1   else     type   0) 

";z 

320  IF  Z=l  THEN  ZER0( I . 5) =  PT*2  ELSE  ZERG(1.5>=0 

330   INPUT"         U.y  =  0      (yes type  1   else    type   0) 

":Z 
340  IF  Z  =  l  THEN  ZER0( I .6) =  PT*5  ELSE  ZER0(I.6)=0 

350   INPUT"        V.x  =  0      (yes type  1   else     type   0) 

"  S  Z 
360  IF  Z=l  THEN  ZERO ( 1 . 7) =  PT*4  ELSE  ZER0(I,7>=0 
370   INPUT"         SIGMA   X  =  0  AND  SIGMA  Y  =  0   (yes type  1 

else    type  0)  ":Z 
380  IF  Z  =  l  THEN  390  ELSE  400 

390  ZER0(I,8)=PT+3:ZER0(I.9)=PT+6:G0T0  410 
400  ZER0(I.8)=0:ZER0(I„9>=0 
410  NEXT  I 
420  FOR  1=1  TO  HN 
430  LPRINT  "ACTND(":I;"1)=";ACTND(I,1) ;TAB(20) : "ACTND ( " : I : "2) 

= " ; ACTND (1.2): TAB ( 40 ) : "ACTND ( " ; I ; "3 ) = " : ACTND (1.3) 
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440  NEXT  I 

450  CBN  ( 1 )  =  -  M  :  CBN  (  2 )  =  - 1 :  CBN  <  3  >  =  -  M 

460  OPEN  "R".#4."3:FCDATAN",4 

470  FIELD  #4,4  AS  Ds 

480  OPEN  "R-.tfa.-BrSTDATAN"^ 

490  FIELD  #3,4  AS  AS 

500  FOR  CASE=1  TO  3 

510  FOR  1=1  TO  HN 

520  PRINT  I 

530  IK=0:IK1=0 

540  IF  ACTNDCI.CA5E>>0  THEN  550  ELSE  1080 

550  XN=ACTND(I,CASE) 

560  IF  CASE=1  THEN  XNl=XN-3 

570  IF  CASE=2  THEN  XN1=XN»1 

580  IF  CASE=3  THEN  XN1=XN*3 

590  IK=(XN-l)»XN/2 

600  FOR  J=l  TO  XN 

610  PT=IK*J 

620  GET  #3,PT 

630  WKKJ>=CVS<AS)  :  WK3  <J>  =WK1  (J) 

640  NEXT  J 

650  IKl=CXNl-l)«XNl/2 

660  FOR  J=l  TO  XN1 

670  PT1=IK1*J 

680  GET  #3.PT1 

690  WK2(J)=CVS(AS> :WK4CJ>=WK2cJ) 

700  NEXT  J 

710  FOR  K  =  XN  TO  TL-1 

720  PT=K«<K*1)/2*XN 

730  GET  #3,PT 

740  WKKK*1)=CVSCAS)  :  WK3  <K»1 )  =  WK1  <K*1 ) 

750  NEXT  K 

760  FOR  K  =  XN1  TO  TL-1 

770  PT1=K»(K-1)/2*XN1 

780  GET  #3.PT1 

790  WK2(K*1)=CVS(AS) : WK4 <K* 1 ) =WK2 (K-l ) 

800  NEXT  K 

810  FOR  J=l  TO  TL 

820  WK1<J)=WKKJ)*WK2(J)«CBN(CASE) 

830  NEXT  J 

840  WK4<XNJ=WKKXN1) 

850  WK3(XN)=WKKXN) 

860  FOR  J=l  TO  TL 

870  WK3(J)=WK3(J) *WK4 < J ) »CBN ( CASE) 

880  NEXT  J 

890  IK=(XN-1) »XN/2 

900  FOR  K=l  TO  XN-i 

910  PT=IK-K 

920  LSET  AS=MKSS(WK1(K) ) 

930  PUT  #3,PT 

940  NEXT  K 

950  FOR  K=XN  TO  TL 

960  PT=<K-1)«K/2*XN 
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970  LSET  A«»MKSS(WK3(K) ) 

980  PUT  #3.PT 

990  NEXT  K 

1000  GET  #4.XN 

1010  3B=CVSCDS) 

1020  GET  #4,XN1 

1030  BB1=CVS(DS) 

1040  BB=BB+BB1»CBN(CASE) 

1050  LSET  DS=MKSS(BB) 

1060  PUT  #4,XN 

1070  BB=0:BB1=0 

1080  NEXT  I 

1090  NEXT  CASE 

1100  CLOSE  #3 

1110  CLOSE  #4 

1120  PT=0 

1130  FOR  1=1  TO  HN 

1140  FOR  J  =  l  TO  9 

1150  IF  ZERO<I,J>>0  THEN  1160  ELSE  1180 

1160  PT=PT»1 

1170  ZEROKPT)=ZEROCI,J) 

1180  NEXT  J 

1190  NEXT  I 

1200  FOR  1=1  TO  PT 

1210  FOR  J  =  l  TO  PT-1 

1220  IF  ZER0KJ>>ZER01(J»1)  THEN  1230  ELSE  1260 

1230  TEMP  =  ZEROKJ) 

1240  ZER0KJ>=ZER01<J*1) 

1250  ZER0KJ»1)=TEMP 

1260  NEXT  3 

1270  NEXT  I 

1280  FOR  1=1  TO  PT 

1290  LPRINT  "ZER01<"»Ij")«"»ZBR01(I3 

1300  NEXT  I 

1310  REROW  =  TL-PT:REEL  =  REROW»  <RER0W<-1)  /2 

1320  OPEN"R" , #2 , "B : RSTDATA" , 4 

1330  FIELD  #2.4  AS  BS 

1340  FOR  1=1  TO  REEL 

1350  B=0 

1360    LSET    BS  =  MKSS8<B> 

1370  PUT  #2,1 

1380  NEXT  I 

1390  CLOSE  #2 

1400  0PEN"R" , #1 , "B : RFCDATA" . 4 

1410  FIELD  #1,4  AS  CS 

1420  FOR  1=1  TO  REROW 

1430  B=0 

1440  LSET  CSS  =  MKSS<B) 

1450  PUT  #1,1 

1460  NEXT  I 

1470  CLOSE  #1 

1480  ZZ=0:ZZ1=0 

1490  OPEN"R" , #4 . "B : FCDATAN" , 4 


98 


1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 


FIELD  #4.4  AS  DS 
0PEN"R" . #3 . "B : STDATAN" . 4 

AS  AS 

"B:RSTDATA".4 

AS  BS 

••B:RFCDATA".4 

AS  CS 


FIELD  #3.4 

0PEN"R",#2 

FIELD  #2.4 

0PEN"R",#1 

FIELD  #1.4 

FOR  1=1  TO  TL 

PRINT  I 

FOR  IJ=1  TO  PT 

RN=ZER01CIJ) 

IF  I=RN  THEN  1800 

NEXT  IJ 

FOR  J=l  TO  I 

FOR  IJ-1  TO  PT 

NC  =  ZEROKIJ) 

IF  J=NC  THEN  1740 

NEXT  IJ 

ZZ=ZZ*1 

IJK»I«<I-1)/2*J 

GET  #3.  UK 

ST=CVS<AS) 

LSET  BS=MKSS<ST) 

PUT  #2.ZZ 

NEXT  J 

ZZ1=ZZ1»1 

GET  #4,1 

FC=CVS(DS) 

LSET  Cs=MKSS<FC> 

PUT  #1.ZZ1 

NEXT  I 

CLOSE  #1 

CLOSE  #2 

CLOSE  #3 

CLOSE  #4 

END 


99 


10  REM 

20  REM       PROGRAM  NAME:  SKYLINE 

30  REM       ..................... 

40  REM   THIS   PROGRAM   FIND  THE  SOLUTION   VECTORS   BY   USING 

SKYLINE  TECHNIQUE. 
50  REM 

60  INPUT"TOTAL  R0W5":N 
70  NP=N-l:N0=N+2 
80  DIM  MK242)  ,PR<242>  ,X<240> 
90  OPEN  ••R",#3."B:RSTDATA",4 
100  FOR  1=1  TO  N 
110  IK=I»<I-l)/2 
120  FOR  J  =  l  TO  I 
130  FIELD  #3.4  AS  AS 
140  IJ=IK+J 
150  GET  #3.IJ 
160  B=CVS(AS) 
170  IF  B=0  THEN  210 

180  IF  BXO  AND  MI(I>=0  THEN  MI(I)=J 
190  PRINT  "11I(";I:")=":«I(I) 
200  GOTO  220 
210  NEXT  J 
220  NEXT  I 
230  CLOSE  #3 

240  OPEN  "R".#4,"B:RFCDATA",4 
250  FOR  1=1  TO  N 
260  FIELD  #4,4  AS  BS 
270  GET  #4,1 
280  A=CVS<BS) 
290  LPRINT  "FC("iI;")="sA 
300  NEXT  I 
310  FOR  1=1  TO  N 
320  FIELD  #4,4  AS  BS 
330  GET  #4,1 
340  A=CVS(BS) 
350  IF  A=0  THEN  370 

360  MI(NP)=I:PRINT  "Mil t" :NP: "}■" :MI (NP) : GOTO  380 
370  NEXT  I 
380  CLOSE  #4 
390  M=l 

400  FOR  1=1  TO  N 
410  M1=I-MI(I)*1 
420  IF  M<M1  THEN  M=M1 
430  NEXT  I 
440  PR<1)=0 
450  FOR  1=1  TO  N 
460  PR<I-1)=PR<I)*I-MI(I)-H 
470  PRINT  "PR<":I*l:',)  =  "iPR<I  +  l) 
480  NEXT  I 

490  PR(NQ)=PR(NP)*N-MI(NP> *1 
500  NT=PR(NQ) 
510  DIM  S(14000> 
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IK=";IK 


:IJ 


,#4,,•B:RFCDATA,' 
4  AS  BS 


520  LPRINT  "NT  =  ,,:NT 

530  OPEN  "R".#3,"B:RSTDATA 

540  FIELD  #3,4  AS  AS 

550  IJ=0 

560  FOR  1  =  1  TO  N 

570  PRINT  "I=":l 

580  J1  =  MKI) 

590  K=PR<I) 

600  IK=I*<I-1)/2:PRINT 

610  FOR  J=J1  TO  I 

620  IJ=IK-J:PRINT  "IJ= 

630  K=K»1 

640  GET  #3, I J 

650  S(K)=CVS(AS) 

660  PRINT  ,,SK<,":K:">=":S(K) 

670  NEXT  J 

680  NEXT  I 

690  CLOSE  #3 

700  OPEN  "R" 

710  FIELD  #4, 

720  KK=PR<NP) 

730  K=0 

740  II  =  MKNP> 

750  FOR  J=II  TO  N 

760  GET  #4. J 

770  K  =  K*1 

780  S(KK*K)=CVS<BS) 

790  PRINT  "SKC":KK-K;")=":S(KK*K) 

800  NEXT  J 

810  CLOSE  #4 

820  REM  The  following  is  the  elimination 

830  N1=N-1 

840  FOR  1=2  TO  N 

850  I2=I*M-l:IF  I2>N  THEN  I2=N 

860  FOR  J=I  TO  12 

870  R1  =  MKI):IF  MKJ>>MI(I)  THEN  R1  =  MI(J) 

880  I1=I-1:IF  R1>I1  THEN  940 

890  FOR  R=R1  TO  II 

900  IJ=PR<J)*I-MI(J)*l:RI=PR(I)*R-MI(I)*l:RJ=PR(J)*R-MI<J)*l 

910  RR=PR(R)-R-MI(R)tl 

920  S<IJ)=S(IJ)-S<RI) »S<RJ)/S(RR) 

930  NEXT  R 

940  NEXT  J 

950  R1=MKNP>:IF  MICI)>MI(NP)  THEN  R1=MI(I) 

960  IF  R1>I-1  THEN  1030 

970  FOR  R=R1  TO  II 

980  RI=PR<I)-R-MI<I)-H:RR=PR<R)*R-MI(R)*l:IJ=PR(NP)*I-MI<NP) 

♦  1 
990  RJ  =  PR(NP)*R-MKNP)*1 
1000  S<IJ)=S(IJ)-S<RI)»S<RJ)/S<RRJ 
1010  PRINT  S<IJ) 
1020  NEXT  R 
1030  NEXT  I 
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1040  FOR  1=1  TO  MICNP)-1:X<I)=0:NEXT  I 

1050  FOR  I=MI(NP)  TO  N   :RI=PR (NP) -I-MI (NP> *1 :X<I> =S<RI> : NEXT 

I 
1060  RR=PR(N) ♦N-MI<N)*1 
1070  X(N)=X<N)/S(RR) 
1080  FOR  11=1  TO  Nl 
1090  NR=N*1-II:XR=XCNR) 

1100  J1=MKNR) :J2=NR-l:lR=PR<NR)-MICNR)*l  • 
1110  FOR  I=J1  TO  J2 
1120  IJ=IR*I 
1130  X<I)=X<I)-S<IJ)«XR 
1140  NEXT  I 

1150  RR  =  PR<J2)"-J2-MI(J2)-1 
1160  X(J2)=XCJ2)/S(RR) 
1170  NEXT  II 

1180  LPRINT  "SOLUTION  VECTOR" 
1190  FOR  1=1  TO  N 

1200  LPRINT  "X(":I;"> ":XCI> 

1210  NEXT  I 

1220  OPEN  "R".#5,"B:DISP",4 

1230  FIELD  #5,  4  AS  DS 

1240  FOR  1=1  TO  N 

1250  LSET  DS  =  MKSsp<X<I>  ) 

1260  PUT  #5.1 

1270  NEXT  I 

1280  CLOSE  #5 

1290  END 
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10  REM 

20  REM  PROGRAM  NAME:  REPARAM 

30  REM  ..................... 

40  REM   This  program  adds  all  known  nodal  parameters  to   the 
reduced  nodal  parameters. 

50  DIM  ACTND(50> .ZER0K1O0) ,U<300> ,UU<300) 

60  M=.17 

70  INPUT  "TOTAL  NODE  N0.":NP 

SO  N=NP»6 

90  INPUT"H0w  MANY  NODAL  PARAMETERS  MUST  BE  MODIFIED" : MP 

100  FOR  1=1  TO  MP 

110  PRINT  I 

120  INPUT  ACTND(I) 

130  NEXT  I 

140  INPUT  "NO.  OF  GIVEN  PARAMETERS" :NGP 

ISO  FOR  1=1  TO  NGP 

160  PRINT   I 

170  INPUT  ZEROKI) 

160  NEXT  I 

190  NN=N-NGP 

200  OPEN  "R",#5."B:DISP".4 

210  FIELD  #5.4  AS  ES 

220  FOR  1=1  TO  NN 

230  GET  #5.1 

240  U(I)=CVS<ES) 
250  NEXT  I 

260  CLOSE  #5 

270  FOR  1=1  TO  NGP 

280  ZER0=ZER01CI) 

290  FOR  J=I  TO  N-l 

300  IF  J=ZERO  THEN  320  ELSE  310 

310  NEXT  J 

320  FOR  JJ=J  TO  N 

330  UU<JJ)=U<JJ) 

340  NEXT  JJ 

350  U<ZER0)=0 

360  FOR  K=J*1  TO  N 

370  U<K)=UU(K-1) 

380  NEXT  K 

390  NEXT  I 

400  FOR  1=1  TO  N 

410  PRINT  I.U(I) 

420  NEXT  I 

430  FOR  1=1  TO  MP 

440  Z=0 

450  Z=ACTND<I) 

460  IF  INT<Z/6)=Z/6  THEN  490 

470  IF  INT< (Z»2>/6>=<Z*2)/6  THEN  500 

480  IF  INT< (Z-3)/6)=(Z»3)/6  THEN  510 

490  U<Z-3)=U<Z)«<-M> :G0T0  520 

500  ucz*i)=-U(Z) :goto  520 

510  U<Z*3)=U(Z)«<-M) :G0T0  520 


103 


520  NEXT  I 

530  OPEN  "R".#6."B:RDISP".4 

540  FIELD  #6,4  AS  FS 

550  FOR  1=1  TO  N 

560  LSET  FS=MKSS(UCI>) 

570  PUT  #6,1 

580  NEXT  I 

590  CLOSE  #6 

600  LPRINT  "SOLUTION  VECTORS" 

610  LPRINT  "•••••»•••••«••»« ":LPRINT 

620  Z*0 

630  FOR  1=1  TO  N-5  STEP  6 

640  Z=2*l 

650  LPRINT  "U<":Z:"> 

660  LPRINT  "V<";2:") 

670  LPRINT  "DU/DX/(";Z:") 

680  LPRINT  "DV/DX/<":Z:") 

690  LPRINT  "DU/DY/<":Z:"> 

700  LPRINT  "DV/DY/<"jZ;") 

710  NEXT  I 


:U<I) 

:U<I*1) 

:U<I*2) 

:U(I*3) 

:U(I»4) 

:U<I*5) 
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PROGRAM  NAME 


COMDISP 


10  REM 

20  REM  ...................... 

30  REM   This  program  adds  all  known  nodal  parameters  to   the 

reduced  nodal  parameters. 
40  DIM  SYM<50,2),CBN(6>,U<426> 
50  INPUT  "TOTAL  NODE  NO . " ; NP 
60  N=NP»6 

70  INPUT"H0W  MANY  NODES  ARE  SYMMETRIC" ;NP5Y 
80  FOR  1  =  1  TO  NPSY 
90  PRINT  I 

100  INPUT  "ADDING  NODE  NO. " ;SYM (I , 1) 
110  INPUT  "ORIGINAL  NODE  NO. " ;SYM < I ,2) 
120  NEXT  I 
130  NN=N-NPSY 

140  OPEN  "R",#5,"B:RDISP",4 
150  FIELD  #5,4  AS  E9 
160  FOR  1=1  TO  NN 
170  GET  #5,1 
180  U(I)=CVS(ES) 
190  NEXT  I 
200  CLOSE  #5 

210  CBN(1>=1:CBN(2)=-1:CBN<3)=-1:CBN(4>=1:CBN<5>=1:CBN<6>=-1 
220  FOR  1=1  TO  NPSY 
230  ND=CSYM(I,1>-1>«6 
240  NDl=<SYM(I,2)-l)-6 
250  FOR  J=l  TO  6 
260  PT=ND*J 
270  PT1=ND1*J 
280  U(PT)=UCPT1)*CBN(J) 
290  NEXT  J 
300  NEXT  I 

310  OPEN  "R".#6,"B:RRDISP",4 
320  FIELD  #6,4  AS  FS 
330  FOR  1=1  TO  N 
340  LSET  FS=MKSS<U<I>> 
350  PUT  #6,1 
360  NEXT  I 
370  CLOSE  #6 
380  LPRINT  "SOLUTION  VECTORS" 


;U<I) 

;U<I-1> 

;U<I*2) 

;U<I*3> 

:U(I+4> 

:U<I*5) 


390  LPRINT 
400  Z=0 

»":LPRINT 

410  FOR  1=1 

TO  N-5  STEP  6 

420  2=2*1 

430  LPRINT 

440  LPRINT 
450  LPRINT 
460  LPRINT 
470  LPRINT 

(l 

'nu  /nY  //«•••» 

■DV/DX/C" 
"DU/DY/C" 

2 
Z 

")-- 

»• . 

480  LPRINT 

•DV/DY/<" 

Z 

")  -- 

... 

490  NEXT  I 

500  END 
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10  REM 

20  REP!  PROGRAM  NAME:  STRESS 

30  REM  *»»»»»»#«»»»*#»»»»»» 

40  REM   PROGRAM  TO   FIND  THE  UNIT  STRAIN  OF  CONCRETE   AND   OF 
STEEL 

50  DIM  ND(71,4)  ,NN(4>  ,X(71)  ,Y(7i:>  ,B<3,24:>  ,FX(4>  ,FY(4>  ,PX(4)  , 
PY(4) 

60  DIB  XB(52.14) , YB (52, 14) ,XE (52, 14) . YE (52. 14) . BAR (52) . U (426 > 

,UU(24) 
70  DIM  DST(3) ,EP(3) ,XST(20) ,XND(20) ,SIG(3) ,D(3,3) 
80  REM  x,y  coordinate  of  each  node 
90  NM=71:REM  total  node  numbers 
100  NE=52:REM  no  of  total  elements 
110  PI=3. 1415927* 
120  CSl=C0S(PI/4) 
130  SNl=SIN(PI/4) 
140  EC=1236364! :M=.2:FCP=2720 
150  ES=2E+07:ESEC=ES-EC 
160  D1=EC/(1-M»M) 
170  D(1,1)=D1:D(1,2)=D1«M 
180  D(2,1)=D<1,2) :D<2,2)=D1 
190  D(3,3)=(l-M)»Dl/2 
200  FOR  1=1  TO  NM 
210  READ  X(I) 
220  NEXT  I 
230  DATA      0.      0.   32.5,   32.5,   40.0,  40.0,  40.0,  40.0, 

40.0,  40.0,  40.0 
240  DATA  50.0,  50.0,  50.0,  50.0,  50.0,  50.0,  50.0 
250  DATA  60.0,  60.0,  60.0,  60.0,  60.0,  60.0,  60.0 
260  DATA  120.0,120.0,120.0,120.0,120.0,120.0,120.0 
270  DATA  190.0,190.0,190.0,190.0,190.0,190.0,190.0 
280  DATA  260.0,260.0,260.0,260.0,260.0,260.0,260.0 
290  DATA  320.0,320.0,320.0,320.0,320.0,320.0,320.0 
300  DATA  330.0,330.0,330.0,330.0,330.0,330.0,330.0 
310  DATA  340.0,340.0,340.0,340.0,340.0,340.0,340.0 
320  DATA  347.5,347.5,380.0,380.0 
330  REM 

340  FOR  1=1  TO  NM 
350  READ  YCI) 

43.5,193.5,233.5 

43.5,53.5,63.5,118.5,193.5,233.5 
43.5,53.5,63.5,118.5,193.5,233.5 
43.5,53.5,63.5,118.5,193.5, 233 . 5 
43.5,53.5,63.5,118.5,193.5,233.5 
43. 5, 53. 5, 63. 5, 118. 5, 193. 5, 233 . 5 
43.5,53.5,63.5,118.5,193.5,233.5 
43.5,53.5,63.5,118.5,193.5,233.5 
43.5,53.5,63.5,118.5,193.5, 233 . 5 
43 . 5 , 53 . 5 , 63 . 5 , 1 18 . 5 , 193 . 5 , 233 . 5 
233.5,   0.0,  43.5 


360 

NEXT 

I 

370 

DATA 

0 

.0, 

380 

DATA 

0 

.0, 

390 

DATA 

0 

.0, 

400 

DATA 

0 

.0, 

410 

DATA 

0 

.0, 

420 

DATA 

0 

.0, 

430 

DATA 

0. 

.0, 

440 

DATA 

0 

,0, 

450 

DATA 

0 

,0. 

460 

DATA 

0, 

.0, 

470 

DATA 

193, 

,5. 

480 

REM 

106 


490  REM  NODE  NO.  CORRESPONDING  TO  EACH  ELEMENT 

500  FOR  EL»1  TO  NE 

510  FOR  EN  =  1  TO  4 

520  READ  NDCEL.EN) 

530  NEXT  EN: NEXT  EL 

540   DATA    1.5.6,2.  5,12.13.6,  12,19,20,13,  19,26,27.20. 

26,33.34,27 
550  DATA  33,40,41,34,  40.47,48,41,  47,54.55,48,  54,61,62,55. 

61,70,71,62 
560  DATA  6,13.14,7,  13,20,21,14,  20,27,28,21,  27,34,35.28, 

34,41,42,35 
570  DATA  41,48,49,42,  48,55.56,49,  55,62,63,56,  7.14,15,  8. 

14,21,22,15 
580  DATA  21,28.29,22.  28,35,36,29,  35,42,43,36,  42,49,50,43. 

49,56,57,50 
590  DATA  56,63,64,57,  8,15,16,9,  15,22,23,16,  22,29,30,23, 

29,36,37,30 
600  DATA  36,43,44,37,  43,50,51,44,  50,57,58,51,  57,64,65,58, 

9,16,17,10 
610  DATA  16,23,24,17,  23,30,31,24,  30,37,38,31,  37,44,45,38, 

44,51,52,45 
620  DATA  51,58,59,52,  58,65,66,59,  3,10,11,4,  10,17,18,11, 

17,24,25,18 
630  DATA  24.31,32,25,  31,38.39,32,  38,45,46.39,  45,52,53,46, 

52,59,60,53 
640  DATA  59,66,67,60,  66,68,69,67 
650  REM  bar  coordinate 
660  FOR  1=1  TO  NE 
670  READ  BAR(I) 
680  NEXT  I 
690  DATA  0,6,6,10,10,10,10,6,6,0,3.3.6.7,7,6.3,3,3,3,6,7,7,6, 

3,3,7,8,12,12,12,12 
700  DATA  8,7.10.10,12,14,14,12,10,10,0,6,6,8,10,10,8,6,6,0 
710  FOR  1=1  TO  NE 
720  BN=BAR(I) 
730  IF  BN=0  THEN  770 
740  FOR  J=l  TO  BN 
750  READ  XBCI,J) 
760  NEXT  J 
770  NEXT  I 
780  DATA  43,43,43,43,43,43,50,50,50,50,57,57,60,60,70,90,110, 

60,60,70,90,110,120,120,130,150 
790  DATA  170,120,120,130,150,170,190,190,210,230,250,190,190, 

210 , 230 , 250 , 260 , 260 
800  DATA  270,290,310,260,260,270,290,310,320,320,320,320.323, 

323 , 330 , 330 , 330 , 330 , 337 , 337 , 43 , 43 , 43 , 50 , 57 , 57 
810  DATA  70,90,110,60,80,100,130,150,170,120,140,160,180,190, 

210,230,250,200 
820  DATA  220,240,270,290,310,260,280,300,320,323,323,330,337, 

337,43,43,43,50,57,57,60,80, 100,70,90, 110, 120 
830  DATA  140,160,180,130,150,170,200,220,240,190,210,230,250, 

260 , 280 , 300 , 270 , 290 
840  DATA  310,320,323,323,330,337,337,43,43,43,43,43,43.43,50, 
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50  ,  50 ,  50 ,  50 ,  55 ,  57 ,  57 ,  60 ,  60 ,  60 ,  70 ,  90 ,  1 10 ,  60  ,  60 
850  DATA  60,75,95,115,120,120.120,130,150,170,120,120,120,135 

,155,175,190,190 
860  DATA  190,210,230,250,190,190,195,215,235,255,260.260,260, 

270,290,310,260,260 
870  DATA  260,275,295,315,320,320,320.320,320,320.323.323.330. 

330 , 330 , 330 , 330 , 337 ,337,43,43,43,45 
880  DATA  43,43,43,43,43,43,50,50,50,50,50,50,50,50,57,57,60,6 

0 ,  60 ,  65 ,  85  ,  105 ,  60 ,  60 ,  60 ,  60 ,  80 
890  DATA  100,120,120,120,125,145,165,185,120,120,120,120,140, 

160,180,190,190,190 
900  DATA  190,205,225.245,190,190,190.190,200.220,240,260.260. 

260 ,  265 ,  285 ,  305 ,  260 
910  DATA  260,260,260,280,300,320,320.320,325,320.320,320,320, 

323 , 323 , 330 , 330 , 330 , 330 , 330 
920  DATA  330,330,330,337,337,43,43,43.43,43,43,50,50,50,50,57 

,  57 ,  60 ,  60 ,  80 ,  100 ,  60 ,  63 ,  83 ,  103 , 1 20 ,  120 
930  DATA  140,160,180,120,123.143,163,183,190,190,200,220,240. 

190,190,203,223,243 
940  DATA  260,260,280,300,260,263,283,303,320,320,320,323,323, 

323 , 330 , 330 , 330 , 330 
950  DATA  337,337 
960  FOR  1=1  TO  NE 
970  BN=BAR<I) 
980  IF  BN  =  0  THEN  1020 
990  FOR  J  =  l  TO  BN 

1000  READ  YBCI, J) : YBCI, J)=YB<1, J)*,5 
1010  NEXT  J 
1020  NEXT  I 
1030  DATA  36,16,10,30,2.5.2.5,23,3,23,43,2.5.2.5.33,13,3.3.3. 

13,33,43,43,43,33,13,3,3,3,13,33,43,43 
1040  DATA  43,23,3,3,3,3,23,43,43,43,43,33,13,3,3,3,13.33,43,4 

3,43,33,13,13,33,2.5,2.5,23 
1050  DATA  3,23,43,2.5,2.5,50.43,43,43,43,43,43,43,43,53,53,53 

,43,43,43,53,53,53,53,43,43,43,43,53,53 
1060  DATA  53,43,43,43,53,53.53,53,43,43,43,43,43,56,53,53,63, 

53,53,53,53,53,63,63,63,53,53,53,53,63,63 
1070  DATA  63,53,53,53,63,63,63,63,53,53,53,63,63.63,53.53,53, 

63 , 53 , 53 , 96 , 76 , 70 , 90 , 1 10 , 63 , 63 , 103 , 83 
1080  DATA  63,83,103,118,63,63,113,93,73,63,63,63,73.93,113,11 

8,118,118,113,93,73,63,63 
1090  DATA  63,73,93,113,118,118,118,103,83,63,63,63,63,83,100. 

118,118,118,118,113 
1100  DATA  93,73,63,63,63,73,93,113.118,119,118,113,93,73,73.9 

3,113,63,63,103,83,63.83 
1110  DATA  103,63,63,176,156,136.118,130,150,170,190,118,118,1 

83,163,143,123,123,143,163,183,118,118 
1120  DATA  173,153,133,118,118,118,133,153,173,193,193,193,173 

,153,133,118,118 
1130  DATA  118, 118, 133. 153, 173. 193, 193, 193, 193, 183, 163, 143. 123 

,118,118,118,123 
1140  DATA  143,163,183,193,193,193,173,153,133,118,118,118,133 
,153,173,193,193 
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1150  DATA  193,173,153,133,113.133,153.173,193,118,118,133.163 

,143,123,123,143,163,183,118, 118 
1160  DATA  216,196,210,230,193,193,223,203,203,223,193,193,213 

, 193,193,193,213,230,230,230,213 
1170  DATA  193,193,193,193,213,230,230,230,230,223,203,193,193 

,193,203,223,230 
1180  DATA  230.230,213,193.193,193,213,230,230,230,213.193,213 

,230,193, 193,223,203,203 
1190  DATA  223,193,193 
1200  FOR  1=1  TO  NE 
1210  BN=BAR(I> 
1220  IF  BN  =  0  THEN  1260 
1230  FOR  J=l  TO  BN 
1240  READ  XECI.J) 
1250  NEXT  J 
1260  NEXT  I 
1270  DATA  50,50,50,50,43,43,60,60,60,60.57,57,70,90,110,120.1 

20,70,90,110,120,120,130,150 
1280  DATA  170,190,190,130,150,170,190,190,210,230,250,260,260 

,210,230,250,260 
1290  DATA  260,270.290.310.320.320,270,290.310.320.320.330,330 

, 330 , 330 , 323 , 323 , 337 ,337 
1300  DATA  337,337,337,337,50,43,43,60,43,43,80,100,120,70,90, 

110,140,160,180,130,150,170,190,200,220 
1310  DATA  240,260,210,230,250,280,300,320,270,290,310,330,323 

,323,337, 337 , 337 , 50 , 43 , 43 , 60 , 57 , 57 , 70 , 90 , 1 10 , 80 
1320  DATA  100,120,130,150,170,190,140,160,180,210,230,250,200 

, 220 , 240 , 260 , 270 , 290 
1330  DATA  310,280,300,320,330,323,323,337,337,337,50,50,50,50 

, 50 , 43 , 43 , 60 , 60 , 60 , 60 , 60 , 60 , 57 , 57 , 65 , 85 , 105 , 120 
1340  DATA  120,120,70,90,110,120,120,120,125,145,165,185,190,1 

90,130,150,170,190 
1350  DATA  190,190,205,225,245,260,260,260,210,230,250,260,260 

,260,265,285,305 
1360  DATA  320 , 320 , 320 , 270 , 290 , 310 , 320 , 320 , 320 , 325 . 330 , 330 . 330 

,330,330,323,323,337,337 
1370  DATA  337,337,337,337,337,50,50,50,50,50,50,50,50,43,43,6 

O , 60 , 60 , 60 , 55 , 60 , 60 , 60 , 57 , 57 , 80 , 100 , 120 
1380  DATA  120,120,120,75,95,115,120,120,120,140,160,180,190,1 

90,190,190,135,155 
1390  DATA  175,190,190,190,190,200,220,240,260,260,260.260.195 

.215,235,255,260 
1400  DATA  260,260.280.300.320,320,320,320,275,295.315.320.320 

,320,330,330,330 
1410  DATA  330,330,330.330,330,323,323.337,337.337,337,335.337 

, 337 , 337 , 337 , 337 , 50 , 50 , 50 , 50 , 43 , 43 , 57 . 60 
1420  DATA  60,60,57,57,77,97,117,120,80,100,120,120,137.157,17 

7,190,190,140,160,180,190 
1430  DATA  190,197,217,237,257,260,200,220,240,260,260,277,297 

,317,320,280,300 
1440  DATA  320,320,330,330,330,330,323,323,337,337,337,337,337 

,337 
1450  FOR  1=1  TO  NE 
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1460  BN=BAR(I> 

1470  IF  BN=0  THEN  1510 

1480  FOR  J*l  TO  BN 

1490  READ  YE<I, J) :YE(I,J>=YE(I, J)*.5 

1500  NEXT  J 

1510  NEXT  I 

1520  DATA  43, 23. 3.23, 43, 43, 33,13. 13. 33, 43, 43.43, 43. 43, 33, 13. 3 

,3,3,13,33,43,43,43,43,23,3,3,3 
1530  DATA  3,23,43,43,43,33,13,3.3,3,13,33,43,43,43.33,13,3,3, 

3,13,33,43,23,3,23,43,43 
1540  DATA  30,10,16,36,43,43.43,53,53,53,53,53,53,53,53,43,43, 

43,53,53,53,43,43,43,43,53,53,53,53 
1550  DATA  43,43,43,53,53,53,43,43,43,43,53,53,50,53,53,63,63, 

63 , 53 , 63 , 63 , 63 , 63 , 63 , 53 , 53 , 53 , 63 , 63 , 63 , 63 
1560  DATA  53,53,53,63,63,63,53,53,53,53,63,63,63,53,53,53,63, 

63,63,56,63,63, 103,83,63,83, 103, 118, 118 
1570  DATA  113,93,73,73,93.113.118,118,118,118,118,113,93,73,6 

3,63,63,73,93,113,118,118 
1580  DATA  118,118,103,83,63,63,63,63,83,100,113,118.118,113,9 

3,73,63,63,63,73,93 
1590  DATA  113,118.118.118.113,93,73,63,63,63,73,93,113.118.10 

3,83,63,83,103,118,118,110 
1600  DATA  90,70,76.96,118,118,183,163,143,123,123,143,163,183 

,193,193,193,173,153,133,118,133 
1610  DATA  153,173,193,193,193,193,193,173,153,133,118,118,118 

,133,153,173,193,193, 193 
1620  DATA  183,163,143,123,118,118,118,123,143,163,183,193,193 

,193,193,173,153 
1630  DATA  133,118,118,118,118,133,153,173,193,193,193,173,153 

,133,118,118,118 
1640  DATA  133,153,173,183.163,143.123.123.143.163,183.193.193 

,190,170,150, 130,118,136 
1650  DATA  156,176,193,193,223.203,203,223.230,230,230,213,193 

,213, 230 , 230 , 230 , 230 , 230 ,213,193,193,193 
1660  DATA  213,230,230,230,223,203,193,193,193,203,223,230,230 

,230,230,213,193 
1670  DATA  193,193.193.213,230,230,230,213,193,193,193.213.223 

,  203 ,  203 ,  223 ,  230 ,  230 ,  230 
1680  DATA  210,196,216,230,230 
1690  FOR  1=1  TO  NE 
1700  BN=BAR<I) 
1710  FOR  J=l  TO  BN 

1720  LPRINT  I,J,XB<I,J> ,YBCI,J> ,XE(I,J) ,YECI,J> 
1730  NEXT  J 
1740  LPRINT 
1750  NEXT  I 

1760  OPEN  "R",#5,"B:RRDISP",4 
1770  FIELD  #5,4  AS  DS 
1780  FOR  1=1  TO  6»NW 
1790  GET  #5,1 
1800  U(I)=CVS(DS) 
1810  NEXT  I 
1820  CLOSE  #5 
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1830  DX=5:DY=5 

1840  FOR  CASE=1  TO  8 

1850  MAX=-1000:MIN=1000 

I860  FOR  EL=27  TO  42 

1870  FOR  K=l  TO  4 

1880  NN(K)=0 

1890  NEXT  K 

1900  FOR  K  =  l  TO  4 

1910  NN(K)=ND<EL,K) : 'PRINT  NN(K) 

1920  NEXT  K 

1930  T  =  0 

1940  FOR  11=1  TO  4 

1950  FOR  JJ  =  1  TO  6 

1960  T  =  T  +  1 

1970  IJ=<NN<II)-1)»6*JJ 

1980  UUCT)=U(IJ) 

1990  NEXT  JJ.II 

2000  XN=<X(NN(2))-X<NNC1)))/5*1 

2010  YN=(YCNNC4))-Y(NN(D)  )/5*l 

2020  FOR  1=1  TO  YN 

2030  XSTCI)=X(NN(1) ) 

2040  XNDCI)=X<NN<2)> 

2050  NEXT  I 

2060  IF  EL>27  THEN  2080 

2070  CLS 

2080  FOR  1=1  TO  4 

2090  J=I*1:IF  J=5  THEN  J=l 

2100  XP=5*X<NN(I))*Y(NN(I))»CS1 

2110  YP=10*Y(NN(I) )»SN1 

2120  XQ=5*X(NN(J))+Y<NN(J)).CS1 

2130  YQ=10*Y(NN(J))«SN1 

2140  LINE  (XP»1. 3-55. 180-YP« . 65) - (XQ»1 . 3-55. 180-YQ*. 65) 

2150  NEXT  I 

2160  IF  CASE>=7  THEN  2460 

2170  FOR  1=1  TO  4 

2180  X1  =  X(NN(D)  :Y1  =  Y<NNCD) 

2190  12=1*1 

2200  IF  12  =  5  THEN  12  =  1 

2210  X2=X<NN<I2>> : Y2=Y (NN (12) ) 

2220  DPX=(X2-X1)/(XN-1) :DPY= C Y2-Y1 ) / ( YN-1) 

2230  IF  DPX  =  0  THEN  ED  =  YN  ELSE  ED  =  XN 

2240  FOR  J=l  TO  ED 

2250  X=X1*(J-1)»DPX 

2260  Y=Y1*(J-1)»DPY 

2270  GOSUB  3200 

2280  NEXT  J 

2290  NEXT  I 

2300  FOR  IK=1  TO  YN 

2310  XO=XST(IK) 

2320  XEND=XND<IK) 

2330  DDX=(XEND-X0)/(XN-1) 

2340  Y=Y(NN(1) )*5»(IK-1) 

2350  FOR  J=l  TO  XN 


Ill 


2350 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 


X  =  X0-*-(J-l.WDDX 

G05UB  3200 

IF  MAX<MGT  THEN  2390  ELSE  2400 

MAX=MGT:XMX=X: YMX  =  Y 

IF  MIN>MGT  THEN  2410  ELSE  2420 

MIN=MGT:XMN=X:YMN=Y 

NEXT  J 

NEXT  IK 

IF  CASE  <7  THEN  2820 

REM  calculate  steel  part 

BN=BAR(EL> 

FOR  EL1  =  1  TO  BN 

XB1=XB(EL,EL1) : Y31= YB (EL, ELI > 

XE1=XE<EL,EL1) : YE1=YE (EL, ELI > 

LENG=SQR( (XE1-XB1>~2+ ( YE1 -YB1 > ~2> 

SN=(YE1-YB1)/LENG 

CS=  (XE1-XBD/LENG 

IF  CASE=7  THEN  2550 

IF  CASE=8  THEN  2560 

IF  SN  =  0  OR  (YE1-YB1X0  THEN  2570  ELSE  2810 

IF  SN>0  OR  (YE1-YB1>>0  THEN  2570  ELSE  2810 

DST(1)=CS~2 

DST(2)=SN*2 

DST(3)=SN«CS 

DX=(XEl-XBl>/20 

DY=(YEl-YBl)/20 

FOR  J=l  TO  21 

X=XB1-DX»(J-1> 

V=X-X(NN(1>> 

L=X(NN<2> )-X(NN(l) ) 

GOSUB  3100 

FXC1)=F1:FX(2)=F2:FX(3)=F3:FX(4}=F4 

PX ( 1 >  =F5 : PX ( 2  >  =F6 : PX ( 3  >  =F7 : PX ( 4  >  =  F8 

Y=YB1*DY»(J-1) 

V=Y-Y(NN(1) ) 

L=Y(NN(4) )-Y(NN(l> ) 

GOSUB  3100 

FY(1)=F1:FY(2>=F2:FY(3)=F3:FY(4>=F4 

P Y ( 1 >  =F5 : P Y ( 2  >  =  F6 : PY ( 3  >  =F7 : P Y ( 4  >  =  F8 

GOSUB  3340 

IF  MAX<ET  THEN  2770  ELSE  2780 

MAX=ET:XMX=X:YMX=Y 

IF  MIN>ET  THEN  2790  ELSE  2800 

MIN=ET:XMN=X: YMN  =  Y 

NEXT  J 

NEXT  ELI 

NEXT  EL 

IF  CASE=1  THEN  2910 

IF  CASE=2  THEN  2930 

IF  CASE=3  THEN  2950 

IF  CASE=4  THEN  2970 

IF  CASE=5  THEN  2990 

IF  CASE=6  THEN  3010 
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2890  IF  CASE=7  THEN  3030 

2900  IF  CASE=8  THEN  3050 

2910  LOCATE  23,5:PRINT  "MAXIMUM  PRINCIPAL  STRESS 

2920  GOTO  3060 

2930  LOCATE  23.5:PRINT  "MINIMUM  PRINCIPAL  STRESS 

2940  GOTO  3060 

23.5:PRINT 


FOR   CRACKINGCOR   CRUSHING)   OF 


,5:PRINT  "NORMAL  STRESS  IN  Y-DIRECTION" 


2950   LOCATE 

CONCRETE" 

2960  GOTO  3060 

2970  LOCATE  23.5:PRINT  "NORMAL  STRESS  IN  X-DIRECTION 

2980  GOTO  3060 

2990  LOCATE  23 

3000  GOTO  3060 

3010   LOCATE  23,5:PRINT  "SHEARING  STRE55  WITH  RE5PECT  TO 

AXES" 
3020  GOTO  3060 

3030   LOCATE   23,5:PRINT  "5TRESS  OF  5TEEL   BAR (horizontal 

greater  than  90  deg>" 
3040  GOTO  3060 

STRESS  OF  STEEL  BAR (vertical  or  1 


X .  Y 


MAX";MAX;"x  coor . " : XMX j 
;XMN;"y  coor . " ; YMN 
THEN  3070 


3050   LOCATE  23,5:PRINT 

than  90  deg>" 
3060  LOCATE  24,9:PRINT 

"MIN";MIN;"x  coor. 
3070  AS=INKEYS:IF  AS  =  "" 
3080  NEXT  CASE 
3090  END 
3100  LV=L-V 

3110  F1=LV/L*V«LV*(L-2»V)/L~3 
3120  F2=V/L*V»LV»(2»V-L)/L~3 
3130  F3=V«LV~2/L~2 
3140  F4  =  -LV»V'S2/L"2 
3150  F5=  (-6»V«L  +  6«V*2)/'L*3 
3160  F6=<6«V«L-6«V-2)/'L*3 
3170  F7=(L~2-4«L»V+3»V-2)/L~2 
3180  F8=(-2«L>V*3«V*2)/LA2 
3190  RETURN 

3200  V1=X-X(NN(1.V)  :L1=X (NN (2) > -X (NN ( 1 > } 
3210  V2=Y-YCNN(1>> : L2=Y (NN (4) ) -Y (NN ( 1 ) ) 
3220  FOR  11=1  TO  2 
3230  V=V1:L=L1 
3240  IF  11=2  THEN  V  =  V2 
3250  IF  11=2  THEN  L=L2 
3260  GOSUB  3100 
3270  IF  11=2  THEN  3310 

3280  FX(1)=F1:FX(2:>=F2:FX(3)=F3:FX(4>=F4 
3290  PX  ( 1  )  =F5  :  PX  ( 2 )  =F6  :  PX  <  3 )  =F7  :  PX  ( 4 )  =F8 
3300  GOTO  3330 

3310  FY<1)=F1:FY<2}=F2:FYC3)=F3:FY(4>=F4 
3320  PY ( 1 ) =F5 : PY (2>  =F6 : PY ( 3) =F7 : PY (4  >  =F8 
3330  NEXT  II 
3340  FOR  11=1  TO  4 
3350  C=(II-1).6 
3360  IF  11  =  1  OR  11=4  THEN  K  =  l  ELSE  K  =  2 


y  coor . 


:  YMX : 
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3370  IF  II  <=2  THEN  KK=1  ELSE  KK  =  2 

3380  B(l,H-C>=PX(K>»FY(KIO  :B  (2.  2  +  C)  =  FX  (K)  »PY<KK> 

3390  B(1,3*C)=PXCK*2)*FYCKK> : B C2 , 4*C> =  FX CK-2) »PY (KK) 

3400  B(1,5*C)=PX(K).FY(KK*2) : B C2, 6*C> =FX <K> »PY (KK-2) 

3410  B(3,1*C)=B(2,2*C) : B (3,2+C) =B< 1 , 1*C) 

3420  B(3,3*C)=B(2,4*C) :B(3,4*C)=B(1,3*C> 

3430  BC3,5*C>=B<2,6<-C>:B<3,6--C)=B<1,5«-C) 

3440  NEXT  II 

3450  FOR  11  =  1  TO  3 

3460  EP(I15=0 

3470  FOR  Jl  =  l  TO  24 

3480  EP(I1)=EP(I1)+B(I1,J1)»UU(J1> 

3490  NEXT  Jl 

3500  NEXT  II 

3510  FOR  11=1  TO  3 

3520  SIG(II)=0 

3530  FOR  KK=1  TO  3 

3540  SIG(II)=SIG(II)+D<II,KK)«EP(KK) 

3550  NEXT  KK 

3560  NEXT  II 

3570  IF  CA5E<7  THEN  3650 

3580  ET=0 

3590  FOR  Jl=l  TO  3 

3600  ET=DST<J1) «EPCJ1)+ET 

3610  NEXT  Jl 

3620  ET=ET»ESEC 

3630  Z=ET:SC=.03 

3640  GOTO  3800 

3650  IF  CASE<4  THEN  3660  ELSE  3780 

3660  HE1=(SIG(1)*SIGC2) )/2 

3670  HE2=(  (SIG(U-SIG(2))/2)  "2 

3680  HE3=(HE2*SIG(3)*2)-.5 

3690  IF  CASE=3  THEN  3730  ELSE  3700 

3700  IF  CASE=1  THEN  3710  ELSE  3720 

3710  MGT=HE1*HE3:Z=MGT:SC=.3:G0T0  3800 

3720  MGT=HE1-HE3:Z=MGT:SC=.3:G0T0  3800 

3730  S1=SIG(1) :S2=SIG<2) :MGT=0 

3740  MGT= .23» ( (S1-S2) -2*S2~2*S1 "23 / <FCP~2/10- .9*FCP* CS1*S2>  > 

3750  WGT=MGT*.54«(S1-.17«S2)»1.1.FCP/(FCP»(FCP*S1+S2)/10».17- 

FCP« (FCP/10-S1-S2) ) 
3760  Z=MGT:SC=50 
3770  GOTO  3800 
3780  MGT=SIG(CASE-3) 
3790  Z=MGT:SC=.3 
3800  XP=5+X*Y*CS1 
3810  YP=10*Y»SN1*SC«2 
3820  IF  J>1  THEN  3850 

3830  LINE  <1.3*XP-55.180-YP*.65>-<1.3*XP-55.180-YP».65> 
3840  GOTO  3860 

3850  LINE  -<1  .3*XP-55,180-YP«.65> 
3860  RETURN 
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10  REM 

20  REM         PROGRAM  NAME  :   DIPLACE 

30  REM         ........................ 

40  REM  This  program  gives  the  graph  of  displacement. 

50  DIM  XC71),Y<71),ND<52,4),NN<4) , U (426) , UU<24) 

60  REM  x,y    coordinate  of  each  node 

70  NM=71:REM  total  node  numbers 

80  NE=52:REM  no  of  total  elements 

90  FOR  1=1  TO  NM 

100  READ  X<I) 

110  NEXT  I 

120  DATA  0.  0,32.5. 

130  DATA 

140  DATA 

150  DATA 


,32.5,40.0.40.0.40.0,40.0.  40.0.  40.0,  40.0 
50.0,  50.0,  50.0,  50.0,  50.0,  50.0,  50. 
60.0,  60.0.  60.0,  60.0,  60.0,  60.0,  60. 
120.0,120.0,120.0,120.0,120.0,120.0,120. 
160  DATA  190.0,190.0,190.0,190.0,190.0,190.0,190. 
170  DATA  260.0,260.0,260.0,260.0,260.0,260.0,260. 
180  DATA  320.0.320.0,320.0,320.0,320.0,320.0,320.0 
190  DATA  330.0,330.0,330.0,330.0,330.0,330.0,330.0 
200  DATA  340.0,340.0,340.0,340.0,340.0,340.0,340.0 
210  DATA  347.5,347.5,380.0,380.0 


220 

REM 

230 

FOR  1=1  TO 

NM 

240 

READ 

Y(I) 

250 

NEXT 

I 

260 

DATA 

0.0, 

43.5, 

193 

.5. 

270 

DATA 

0.0, 

43.5, 

53 

.5, 

280 

DATA 

0.0, 

43.5, 

53 

.5, 

290 

DATA 

0.0, 

43.5, 

53 

.5, 

300 

DATA 

0.0, 

43.5, 

53 

.5. 

310 

DATA 

0.0, 

43.5, 

53 

■  5, 

320 

DATA 

0.0. 

43.5, 

53. 

.5, 

330 

DATA 

0.0, 

43.5, 

53 

,5, 

340 

DATA 

0.0, 

43.5, 

53, 

,5, 

350 

DATA 

0.0, 

43.5, 

53, 

,5, 

360 

DATA 

193.5, 

233.5, 

0, 

.0, 

233.5 
63.5,118.5,193.5,233 
63.5,118.5,193.5,233 
63.5,118.5,193.5,233 
63. 5, 118. 5, 193. 5, 233, 
63.5,118.5,193.5,233, 
63.5,118.5.193.5,233, 
63.5,118.5,193.5,233, 
63.5,118.5,193.5,233, 
63.5,118.5,193.5,233, 
43.5 
370  REM 

380  REM  NODE  NO.  CORRESPONDING  TO  EACH  ELEMENT 
390  FOR  EL=1  TO  NE 
400  FOR  EN=1  TO  4 
410  READ  ND(EL.EN) 
420  NEXT  EN: NEXT  EL 
430  DATA   1,  5.  6.  2. 

26,33,34,27 
440  DATA  33,40,41,34, 

61,70,71,62 
450  DATA   6.13.14,  7. 

34,41,42,35 
460  DATA  41,48,49,42, 

14,21,22,15 
470  DATA  21,28,29,22.  28,35,36.29,  35,42.43.36, 
49,56,57,50 


5.12.13.  6. 
40.47,48,41. 
13,20.21,14, 


12,19,20,13, 
47.54,55,48, 
20,27,28,21, 


48,55,56,49,  55,62,63,56,     7,14 


42,49, 


19,26,27,20, 
54,61,62,55, 
,35,28, 
15,  8, 
50,43, 


27,34, 
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480  DATA  56,63,64,57,   8,15,16,  9,  15,22,23,16,   22,29,30,23, 

29,36,37,30 
490  DATA  36,43,44,37.  43.50,51,44,  50.57.58.51,   57,64,65,58, 

9,16,17,10 
500  DATA  16,23.24,17.  23.30.31,24.  30,37,38,31,   37,44,45,38. 

44,51,52,45 
510  DATA  51,58,59.52.  58.65.66,59.   3.10.11,  4,   10,17,18,11, 

17,24,25,18 
520  DATA  24,31,32.25.  31,38,39.32.  38.45.46.39.   45,52,53.46. 

52,59,60,53 
530  DATA  59,66,67,60,  66,68,69,67 
540  OPEN  "R",#5,"B:RRDISP",4 
550  FIELD  #5,4  AS  DS 
560  FOR  1=1  TO  426 
570  GET  #5,1 
580  UCI)=CVS<DS> 
590  NEXT  I 
600  CLOSE  #5 
610  REM 

620  FOR  EL=1  TO  52 
630  FOR  K  =  l  TO  4 
640  NN(K)=ND(EL,K> 
650  NEXT  K 
660  FOR  11=1  TO  4 
670  FOR  JJ=1  TO  6 
680  T  =  T  +  1 

690  IJ=CNN<II>-1>»6*JJ 
700  UU<T)=U(IJ> 
710  NEXT  JJ.II 
720  IF  EL>1  THEN  740 
730  CLS 

740  FOR  1=1  TO  4 
750  J=I+1:IF  J=5  THEN  J  =  l 
760  XP=20*X<NN<I>> :YP=Y(NN(I>) 
770  XQ=20*X<NN<J)> :YO=Y<NN<J>) 
780  LINE  <XP»1.3,200-YP)-(XQ»1.3,200-YQ) 
790  NEXT  I 
800  FOR  1=1  TO  4 
810  J=I+l:IF  J=5  THEN  J=l 

820  XP=20*X<NN<I))*UO<<I-1)»6*1)«50:YP=Y<NN<I))*UU< (I-l>*6*2> 

•  50 

830  XQ=20*X(NN<J))*UU((J-l)»6*l)»50:Ya=Y«NN(J))*UU((J-l)«6*2) 

•  50 

840  LINE  <XP«1.3.200-YP)-(XQ»1.3.200-YQ) 
850  LOCATE  24, 22: PRINT  "DISPLACEMENT  SHAPE" 
860  AS  =  INKEYS:IF  As3=""  THEN  860 
870  END 
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PROGRAM  NAME:  STRAIN 


10  REM 
20  REM 
30  REM  ».».,.........»...,. 

40   REM   PROGRAM  TO  FIND  THE  UNIT  STRAIN  OF  CONCRETE   AND   OF 

STEEL 

50  DIM  ND<71.4> ,NN<4> ,X<71) ,Y<71) ,B<3.24) ,FX<4) ,FY<4) ,PX(4) , 

PY<4) 
60  DIM  XB<52,14> ,YB(52,14) ,XE(52,14) ,YE<52,14) ,BAR<52> ,U(426) , 

UU<24) 
70  DIM  DST(3) ,EP<3> .X5TC20) ,XND(20> 
80  REM  x,y  coordinate  of  each  node 
90  NM=71:REM  total  node  numbers 
100  NE=52:REM  no  of  total  elements 
110  PI  =  3. 1415927* 
120  CSl=C0S<PI/4> 
130  SNl=SIN<PI/4> 
140  FOR  1=1  TO  NM 
150  READ  X<I) 


40.0 


160 

NEXT 

I 

170 

DATA 

0.0,32.5 

,32.5,40, 

.0,40.0,40 

.0 

.40.0,40.0 

,40.0 

180 

DATA 

50, 

.0,  50, 

.0,  50 

.0, 

.  50, 

.0,  50 

.0, 

.  50, 

.0,  50, 

.0 

190 

DATA 

60 

.0,  60 

.0,  60 

.0, 

,  60 

.0,  60 

.0, 

.  60 

.0,  60 

.0 

200 

DATA 

120 

.0,120, 

.0,120 

.0, 

,120, 

.0,120, 

.0, 

,120 

.0,120, 

.0 

210 

DATA 

190 

.0,190, 

.0,190 

.0, 

,190, 

.0,190, 

.0, 

,190 

.0,190, 

.0 

220 

DATA 

260 

.0,260, 

.0,260 

.0, 

,260, 

,0,260, 

,0, 

,260, 

.0,260, 

.0 

230 

DATA 

320 

.0,320, 

.0,320, 

.0, 

320, 

.0,320, 

,0, 

,320 

.0,320, 

,0 

240 

DATA 

330, 

.0,330, 

,0,330, 

.0. 

330, 

,0,330, 

,0, 

,330, 

.0,330. 

0 

250 

DATA 

340, 

,0,340, 

.0,340, 

.0, 

340, 

.0.340, 

,0, 

340, 

.0,340, 

,0 

260 

DATA 

347, 

.5,347. 

,5,380, 

.0, 

380. 

,0 

270  REM 

280  FOR  1=1  TO  NM 

290  READ  Y<I) 
I 

0.0. 
0.0, 
0.0. 
0.0, 
0.0, 
0.0, 
0.0, 


43. 
43. 
43. 


300  NEXT 
310  DATA 
320  DATA 
330  DATA 
340  DATA 
350  DATA 
360  DATA 
370  DATA 
380  DATA    0.0, 
390  DATA    0.0, 
400  DATA    0.0, 
410  DATA  193.5, 
420  REM 
430  REM  NODE 
440  FOR  EL=1 
450  FOR  EN=1 
460  READ  ND<EL,EN) 
470  NEXT  EN: NEXT  EL 
480  DATA   1,  5,  6,  2 
26,33,34,27 


5,193.5,233.5 

5,53.5,63.5,118.5.193.5,233.5 

5, 53. 5, 63. 5, 118. 5, 193. 5, 233 . 5 

43 .  5 ,  53  . 5 ,  63 .  5 , 1 18  .  5  ,  193 .  5 ,  233  .  5 

43.5,53.5,63.5,118.5,193.5, 233 . 5 

43.5,53.5,63.5,118.5,193.5, 233 . 5 

43 . 5 , 53 . 5 , 63 . 5 , 1 18 . 5 , 193 . 5 , 233 . 5 

43.5,53.5,63.5,118.5,193.5,233.5 

118.5,193.5,233.5 

118.5,193.5,233.5 

5 


43.5,53.5,63.5 

43.5,53.5,63.5 

233.5,   0.0.  43 


NO.  CORRESPONDING  TO  EACH  ELEMENT 
TO  NE 
TO  4 


5,12.13,  6,   12,19,20,13,   19,26,27,20, 
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490  DATA  33.40.41.34,  40.47,48,41.   47,54,55,48,   54,61,62,55 

61.70,71,62 
500  DATA   6,13,14,  7,  13,20,21,14,   20,27,28,21,   27,34,35,28, 

34,41,42,35 
510  DATA  41,48,49,42,  48,55,56,49,   55,62,63,56,    7,14,15,  8 

14,21,22.15 
520  DATA  21,28.29,22,  28,35,36,29,   35,42,43,36,   42,49,50,43, 

49,56,57,50 
530  DATA  56,63,64,57,   8,15,16,  9,   15.22,23,16,   22,29,30,23. 

29,36,37,30 
540  DATA  36,43,44,37,  43,50,51,44,   50,57,58,51,   57,64,65,58, 

9.16,17,10 
550  DATA  16,23,24,17,  23,30,31,24,   30,37,38,31,   37,44,45,38, 

44,51,52,45 
560  DATA  51,58,59,52,  58,65,66,59,    3,10,11,  4,   10,17,18,11, 

17,24,25,18 
570  DATA  24,31,32,25,  31,38,39,32,   38,45,46,39,   45,52,53,46 

52,59,60,53 
580  DATA  59,66,67,60,  66,68,69,67 
590  REM  bar  coordinate 
600  FOR  1=1  TO  NE 
610  READ  BAR<I) 
620  NEXT  I 

630  DATA  0,6,6,10,10,10,10,6,6,0.3,3,6,7,7,6,3,3,3,3,6,7,7,6 
3,3,7,8,12,12,12,12 

640  DATA  8,7,10,10,12,14,14,12,10,10,0,6,6,8,10,10,8,6,6,0 

650  FOR  1=1  TO  ME 

660  BN=BAR<I) 

670  IF  BN=0  THEN  710 

680  FOR  J=l  TO  BN 

690  READ  XB<I,J) 

700  NEXT  J 

710  NEXT  I 

720  DATA  43,43,43,43,43,43,50,50,50,50,57,57,60,60,70,90,110, 
60,60,70,90,110,120,120,130.150 

730  DATA  170,120,120,130,150,170,190,190,210,230,250,190,190, 
210 , 230 , 250 , 260 , 260 

740  DATA  270,290,310,260,260,270,290,310,320,320,320,320,323, 
323,330,330,330,330,337,337,43,43,43,50,57,57 

750  DATA  70,90.110,60,80,100,130,150,170,120,140,160,180,190, 
210,230,250,200 

760  DATA  220,240,270,290,310,260,280,300,320,323,323,330,337, 
337,43,43,43,50,57,57,60,80,100,70,90,110,120 

770  DATA  140,160,180,130,150,170,200,220,240,190,210,230,250, 
260 , 280 , 300 , 270 . 290 

780  DATA  310,320,323,323,330,337,337,43,43,43,43,43,43,43,50, 
50,50,50,50,55,57,57,60,60,60,70,90,110,60,60 

790  DATA  60,75,95,115,120,120,120,130,150,170,120,120,120,135 
,155,175,190,190 

800  DATA  190,210,230,250,190,190,195,215,235,255,260,260,260, 
270 , 290 , 31 0 , 260 ,260 

810  DATA  260,275,295,315,320,320,320.320,320,320,323,323,330, 
330 , 330 , 330 , 330 , 337 , 337 ,43,43,43,45 
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820  DATA  43,43,43,43,43,43,50,50,50,50,50,50,50,50,57,57,60,6 

0 , 60 . 65 , 85 , 105 , 60 , 60 , 60 , 60 , 80 
830  DATA  100,120,120,120,125,145,165,185,120,120,120,120,140, 

160 . 180 , 190 , 190 , 190 
840  DATA  190.205,225,245,190,190,190,190,200,220,240,260,260, 

260 , 265 , 285 , 305 , 260 
850  DATA  260,260,260,280,300,320,320,320,325,320,320,320,320, 

323 , 323 , 330 , 330 , 330 , 330 , 330 
860  DATA  330,330,330,337,337,43,43,43,43,43,43,50,50,50,50.57 

, 57 , 60 , 60 , 80 , 100 , 60 , 63 , 83 , 1 03 , 120 . 120 
870  DATA  140,160,180,120,123,143,163,183,190,190,200,220,240, 

190 , 190 , 203 , 223 , 243 
880  DATA  260,260,280,300,260,263,283,303,320,320,320,323,323, 

323 , 330 , 330 , 330 , 330 
890  DATA  337,337 
900  FOR  1=1  TO  NE 
910  BN=BAR<I> 
920  IF  BN  =  0  THEN  960 
930  FOR  J=l  TO  BN 

940  READ  YB(I,J>:YB<I,J)=YB<I,J)*.5 
950  NEXT  J 
960  NEXT  I 
970  DATA  36,16,10,30,2.5,2.5,23,3,23,43,2.5,2.5,33,13,3,3,3,13 

,33,43.43,43,33.13,3,3,3,13,33,43,43 
980  DATA  43,23,3,3,3,3,23,43,43,43,43,33,13,3,3,3,13,33,43,43, 

43,33,13,13.33,2.5.2.5,23 
990  DATA  3,23,43,2.5,2.5,50,43,43,43,43.43,43,43,43,53,53,53,4 

3 . 43 , 43 , 53 , 53 , 53 , 53 . 43 , 43 , 43 , 43 , 53 , 53 
1000  DATA  53,43,43,43,53,53,53,53,43,43,43,43,43,56,53,53,63,5 

3,53,53,53,53,63,63,63,53,53,53,53,63,63 
1010  DATA  63,53,53,53,63,63,63,63,53,53,53,63,63,63,53,53,53,6 

3,53,53,96,76,70,90,110,63,63,103,83 
1020  DATA  63,83,103,118,63,63,113,93,73,63,63,63,73,93,113,118 

,118,118,113,93,73,63,63 
1030  DATA  63,73,93,113,118,118,118,103,83,63,63,63,63,83,100,1 

18,118,118,118,113 
1040  DATA  93,73,63,63,63,73,93,113,118,118,118,113,93,73,73,93 

, 113 , 63 , 63 , 103 , 83 . 63 , 83 
1050  DATA  103,63,63,176,156,136,118,130,150.170,190.118,118.18 

3.163,143.123,123,143,163.183,118.118 
1060  DATA  173.153,133,118,118.118.133.153,173,193,193,193,173, 

153,133,118,118 
1070  DATA  118,118,133,153,173,193,193,193,193,183.163.143,123, 

118,118,118,123 
1080  DATA  143,163,183,193,193,193.173,153,133,118,118,118,133, 

153,173,193.193 
1090  DATA  193.173,153,133.118.133,153,173,193,118,118,183,163, 

143,123,123.143,163,183,118,118 
1100  DATA  216,196,210,230,193,193,223,203,203,223,193,193,213, 

193,193,193,213,230.230.230.213 
1110  DATA  193,193.193,193,213,230,230,230,230.223.203,193,193, 

193,203,223,230 
1120  DATA  230,230,213,193.193.193,213,230,230,230,213,193.213, 
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230 ,193.193. 223 . 203 . 203 
1130  DATA  223,193.193 
1140  FOR  1=1  TO  NE 
1150  BN=BAR(I) 
1160  IF  BN=0  THEN  1200 
1170  FOR  J=l  TO  BN 
1180  READ  XECI.J) 
1190  NEXT  J 
1200  NEXT  I 
1210  DATA  50,50,50,50,43,43,60,60,60,60,57,57,70,90,110,120, 

120,70,90,110,120,120,130,150 
1220  DATA  170,190,190,130,150,170,190,190,210,230,250,260,26 

0,210.230,250,260 
1230  DATA  260,270,290,310,320,320,270,290,310,320,320,330.33 

0 , 330 , 330 , 323 , 323 , 337 , 337 
1240  DATA  337  ,  337 ,  337  ,  337  ,  50 .  43 ,  43 ,  60 ,  43 ,  43 ,  SO  ,  100 , 1 20 ,  70 ,  90 . 

110,140,160,130,130,150.170,190,200,220 
1250  DATA  240,260,210,230,250,280,300,320.270,290.310.330,323 

,323,337.337,337.50,43,43,60,57,57,70,90,110,80 
1260  DATA  100,120,130,150.170,190,140,160,180,210,230,250,200 

, 220 , 240 , 260 , 270 , 290 
1270  DATA  310,280,300,320,330,323,323,337,337,337,50,50,50,50 

,  50  ,  43, 43 ,  60 .  60 .  60  ,  60 ,  60 ,  60 ,  57 ,  57 ,  65 ,  85 .  105 .  120 
1280  DATA  120,120,70,90,110,120,120,120,125,145,165,185,190,1 

90,130,150,170,190 
1290  DATA  190,190,205,225,245,260,260,260,210,230,250,260,260, 

260,265,285,305 
1300  DATA  320,320,320,270,290,310,320,320,320,325,330,330,330, 

330 , 330 , 323 , 323 , 337 , 337 
1310  DATA  337,337,337,337,337,50,50,50,50,50,50,50,50,43,43,6 

0 ,  60  ,  60 ,  60 ,  55 ,  60  ,  60 ,  60 ,  57 ,  57 ,  80 ,  100 , 1 20 
1320  DATA  120,120,120,75,95,115,120,120,120,140,160,180,190,1 

90,190.190,135,155 
1330  DATA  175,190,190,190,190,200,220,240,260,260,260,260,195 

,215,235,255,260 
1340  DATA  260,260,280,300,320,320,320,320,275,295,315,320,320 

.320.330,330,330 
1350  DATA  330,330,330,330,330,323,323,337,337,337,337,335.337 

,  337 ,  337 ,  337 ,  337 ,  50  ,  50 ,  50  ,  50  ,  43 ,  43 ,  57 ,  60 
1360  DATA  60,60,57,57,77,97,117,120,80,100,120,120,137,157,17 

7 . 190 , 190 , 140 , 160 , 180 , 190 
1370  DATA  190,197,217,237,257,260,200,220.240,260,260,277,297 

,317,320,280,300 
1380  DATA  320.320.330.330.330,330.323,323,337,337,337,337,337 

,337 
1390  FOR  1=1  TO  NE 
1400  BN=BAR(I> 
1410  IF  BN=0  THEN  1450 
1420  FOR  J=l  TO  BN 

1430  READ  YE<I,J)  :  YE  <  I ,  J)  =YE<  I ,  J  )  «•  .5 
1440  NEXT  J 
1450  NEXT  I 
1460  DATA  43,23.3,23,43,43.33,13,13,33,43,43,43,43,43,33,13,3 
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.3.3.13.33.43.43.43.43.23.3.3.3 
1470  DATA  3.23,43.43,43.33,13,3,3,3.13,33.43,43,43,33,13,3,3, 

3.13,33,43,23,3.23,43,43 
1480  DATA  30,10,16,36,43,43,43,53,53,53,53,53,53,53,53,43,43, 

43 . 53 , 53 , 53 , 43 , 43 , 43 , 43 . 53 , 53 . 53 , 53 
1490  DATA  43,43,43,53,53,53,43,43,43,43,53,53,50,53,53,63,63, 

63 ,  53 ,  63 ,  63  ,  63  ,  63 ,  63 ,  53  ,  53 ,  53 ,  63 ,  63 ,  63 ,  63 
1500  DATA  53,53,53,63,63,63,53,53,53,53,63,63,63,53,53,53,63, 

63  ,  63 ,  56 ,  63 ,  63 ,  103  ,  83 ,  63 ,  83 ,  103 , 1 18 , 1 18 
1510  DATA  113,93.73,73,93,113,118,118,118,118,118,113,93,73,6 

3.63,63,73,93,113,118,118 
1520  DATA  118,118,103,83,63,63,63,63,83,100,118,118,118,113,9 

3 , 73 , 63 , 63 , 63 , 73 , 93 
1530  DATA  113,118,118,118,113,93,73,63,63,63,73,93,113,118,10 

3,83,63,83,103,118,118,110 
1540  DATA  90,70,76,96,118,118,183,163,143,123,123,143,163,183 

,193,193,193.173,153,133,118,133 
1550  DATA  153,173,193,193,193,193,193,173,153,133,118.118,118 

,133,153,173,193,193,193 
1560  DATA  183,163,143,123,118,118,118,123,143,163,183,193,193 

,193,193,173.153 
1570  DATA  133,118,118,118,118,133,153,173,193,193,193,173,153 

.133.118.118,118 
1580  DATA  133.153,173,183,163,143,123,123,143,163,183,193,193 

,190,170.150,130,118,136 
1590  DATA  156,176,193,193,223,203,203,223,230,230,230,213,193 

.213,230,230.230,230,230,213,193,193,193 
1600  DATA  213,230,230,230,223,203,193,193,193,203,223,230,230 

,230.230,213,193 
1610  DATA  193,193,193.213,230,230,230,213,193,193,193,213,223 

, 203 , 203 , 223 , 230 , 230 , 230 
1620  DATA  210,196,216,230,230 
1630  FOR  1=1  TO  NE 
1640  BN=BAR(I) 
1650  FOR  J  =  l  TO  BN 

1660  LPRINT  I,J,XB(I,J)  , YB < I , J> , XE < I , J) , YE  CI ,  J> 
1670  NEXT  J 
1680  LPRINT 
1690  NEXT  I 

1700  OPEN  "R",#5,"B:RRDISP",4 
1710  FIELD  #5,4  AS  DS 
1720  FOR  1=1  TO  6»NB 
1730  GET  #5,1 
1740  U<I)=CVS(DSO 
1750  NEXT  I 
1760  CLOSE  #5 
1770  DX=5:DY=5 
1780  FOR  CASE=1  TO  7 
1790  MAX=-1000:MIN=1000 
1800  FOR  EL=27  TO  42 
1810  FOR  K=l  TO  4 
1820  NN<K>=0 
1830  NEXT  K 
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1840  FOR  K=l  TO  4 

1850  NN(K)=ND<EL,K) :'PRINT  NNCK) 

1860  NEXT  K 

1870  T  =  0 

1880  FOR  11=1  TO  4 

1890  FOR  JJ  =  1  TO  6 

1900  T=T*1 

1910  IJ=<NN<II)-1)«6*JJ 

1920  UU<T)=U(IJ) 

1930  NEXT  JJ.II 

1940  XN=CX<NN<2))-X(NN<1))>/5*1 

1950  YN=(Y(NN<4) ) -Y(NN<1> ) )/5*l 

1960  FOR  1=1  TO  YN 

1970  XST(I) =X<NN(1) ) 

1980  XND<I)=X<NN<2) ) 

1990  NEXT  I 

2000  IF  EL>27  THEN  2020 

2010  CLS 

2020  FOR  1=1  TO  4 

2030  J=I*1:IF  J=5  THEN  J=l 

2040  XP=5+X<NN<I) )*Y<NN(I) )*CS1 

2050  YP=10*Y(NN(I) >»SN1 

2060  XQ=5*X(NN<J))+Y<NNCJ) ) »CS1 

2070  YQ=10*Y(NN<J) ) »SN1 

2080  LINE  (XP«1. 3-55. 180-YP» . 65) - (XQ-1 . 3-55. 180-YQ». 65) 

2090  NEXT  I 

2100  IF  CASE>  =  6  THEN  2400 

2110  FOR  1=1  TO  4 

2120  X1=XCNN(I) ) :Y1=Y<NN<I) > 

2130  I2=I*1 

2140  IF  12=5  THEN  12=1 

2150  X2=X(NN<I2) ) :Y2=Y(NN(I2) ) 

2160  DPX=(X2-X1)/<XN-1) :DPY= ( Y2-Y1 ) / < YN-1 ) 

2170  IF  DPX=0  THEN  ED=YN  ELSE  ED=XN 

2180  FOR  J=l  TO  ED 

2190  X=X1*(J-1)«DPX 

2200  Y=Y1+(J-1) »DPY 

2210  GOSUB  3110 

2220  NEXT  J 

2230  NEXT  I 

2240  FOR  IK=1  TO  YN 

2250  XO=XST<IK) 

2260  XEND=XND<IK> 

2270  DDX=<XEND-X0)/(XN-1) 

2280  Y=Y(NN(1))*5*(IK-1) 

2290  FOR  J  =  l  TO  XN 

2300  X=X0*<J-1)»DDX 

2310  GOSUB  3110 

2320  IF  MAX<MGT  THEN  2330  ELSE  2340 

2330  MAX=MGT:XMX=X: YMX=Y 

2340  IF  MIN>MGT  THEN  2350  ELSE  2360 

2350  MIN=MGT:XMN=X:YMN=Y 

2360  NEXT  J 
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2370  NEXT  IK 

23S0  IF  CASE  <6  THEN  2760 

2390  REM  calculate  steel  Dart 

2400  BN=BAR(EL) 

2410  FOR  EL1=1  TO  BN 

2420  XB1=XB(EL,EL1> : YB1=YB (EL, ELI) 

2430  XE1=XE(EL.EL1)  :  YE1  =  YE  <EL,  ELI ) 

2440  LENG=SQR< (XE1-XB1 ) ~2* < YE1 - YB1 > ~2> 

2450  SN=(YE1-YB1)/LENG 

2460  CS=(XE1-XB1)/LENG 

2470  IF  CASE=6  THEN  2490 

2480  IF  CASE=7  THEN  2500 

2490  IF  SN  =  0  OR  (YE1-YB1X0  THEN  2510  ELSE  2750 

2500  IF  SN>0  OR  (YE1-YB1)>0  THEN  2510  ELSE  2750 

2510  DST(1)=CS~2 

2520  DST(2)=SN~2 

2530  DST(3)=SN«CS 

2540  DX=<XEl-XBl)/20 

2550  DY=(YEl-YBl)/20 

2560  FOR  J=l  TO  21 

2570  X=XB1*DX»<J-1) 

2580  V=X-X(NN(1)> 

2590  L=X(NN(2))-X(NN(1) ) 

2600  GOSUB  3010 

2610  FX(1)=F1:FX(2)=F2:FX(3)=F3:FX(4>=F4 

2620  PX ( 1 ) =F5 : PX ( 2) =F6 : PX ( 3  >  =F7 : PX ( 4 ) =F8 

2630  Y=YB1-DY»(J-1> 

2640  V  =  Y-Y(NN(D) 

2650  L  =  Y(NN(4))-Y(NN(D) 

2660  GOSUB  3010 

2670  FY(1)=F1:FY(2)=F2:FY(3)=F3:FY(4)=F4 

2680  PY(1)=F5:PY(2)=F6:PY(3)=F7:PY(4)=F8 

2690  GOSUB  3250 

2700  IF  BAX<ET  THEN  2710  ELSE  2720 

2710  MAX=ET:XMX=X:YMX=Y 

2720  IF  MIN>ET  THEN  2730  ELSE  2740 

2730  MIN=ET:XMN=X:YMN=Y 

2740  NEXT  J 

2750  NEXT  ELI 

2760  NEXT  EL 

2770  IF  CASE=1  THEN  2840 

2780  IF  CASE=2  THEN  2860 

2790  IF  CASE=3  THEN  2880 

2800  IF  CASE=4  THEN  2900 

2810  IF  CASE  =  5  THEN  2920 

2820  IF  CASE=6  THEN  2940 

2630  IF  CASE=7  THEN  2960 

2840  LOCATE  23.5:PRINT  "MAXIMUM  PRINCIPAL  STRAIN  " 

2850  GOTO  2970 

2860  LOCATE  23,5:PRINT  "MINIMUM  PRINCIPAL  STRAIN  " 

2870  GOTO  2970 

2880  LOCATE  23.5:PRINT  "NORMAL  STRAIN  IN  X-DIRECTION' 

2890  GOTO  2970 
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2900 
2910 
2920 

2930 
2940 

2950 
2960 

2970 

2980 
2990 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 


X.Y 


or 


LOCATE  23,5:PRINT  "NORMAL  STRAIN  IN  Y-DIRECTION" 
GOTO  2970 

LOCATE  23.5:PRINT  "SHEARING  STRAIN  WITH  RESPECT  TO 

AXES" 
GOTO  2970 

LOCATE   23.5:PRINT  "STRAIN  OF  STEEL   BAR<horizontal 
greater  than  90  deer)" 
GOTO  2970 

LOCATE  23.5:PRINT  "STRAIN  OF  STEEL  BARCvertical  or  leas 

than  90  deg)" 
LOCATE  24. 9: PRINT  "MAX" ; MAX ; "x  coor . " :XMX ; "y  coor . " ; YMX 
;"MIN":MIN:"x  coor . ":XMN: "y  coor.":YMN 
AS=INKEYS:IF  AS  =  ""  THEN  2980 
NEXT  CASE 
END 

LV=L-V 

F1=LV/L»V»LV» <L-2»V>/L~3 
F2  =  V/L*V»LV«<2»V-L)/L'-3 
F3=V»LV-2/L~2 
F4=-LV«V~2/L~2 
F5=(-6«V»L*6>V~2)/L*3 
F6=(6»V«L-6«V~2)/L~3 
F7=<L~2-4»L»V*3»V~2)/L"2 
F8=(-2*L«V*3»V*2)/L"2 
RETURN 

V1=X-X<NN(1>) :L1=X(NN<2) )-X(NN(l> ) 
V2=Y-Y(NNC1>> :L2=Y(NN<4))-Y(NN(1) ) 
FOR  11=1  TO  2 
V=Vl:L=Ll 
IF  11=2  THEN  V*V2 
IF  11=2  THEN  L=L2 
GOSUB  3010 
IF  11=2  THEN  3220 

FX(l)=Fl:FX(2)=F2:FX<3)=F3:FX(4)=F4 
PX<1)=F5:PX<2)=F6:PX(3)=F7:PX<4)=F8 
GOTO  3240 

FY<l)=Fl:FY(2)=F2:FY<3)=F3:FY<4)=F4 
PY<1)=F5:PY(2)=F6:PY<3)=F7:PY<4)=F8 
NEXT  II 
FOR  11=1  TO  4 
C=(II-1)«6 

IF  11  =  1  OR  11=4  THEN  K=l  ELSE  K  =  2 
IF  II  <=2  THEN  KK=1  ELSE  KK=2 

B<1.1*C)=PX(K)»FY<KK) :B <2 ,2*C) =FX <K) *PY (KK) 
B<1,3*C>=PXCK*2)*FY<KK) : B<2, 4*C> =FX(K*2) *PY <KK) 
B<1.5*C)=PX<K)»FY(KK*2) : B <2,6*C) =FX <K) »PY (KK*2) 
B<3,1*C)=B<2,2*C)  :B<3,2->-C)=B<l,l*C) 
B<3.3*C)=BC2,4*C):B(3,4*C)=B(1,3*C) 
B<3,5+C)=B(2,6*C) : B <3,6*C> =B < 1 , 5-C) 
NEXT  II 
FOR  11=1  TO  3 
EP<I1)=0 
FOR  Jl=l  TO  24 
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3390  EP(I1)=EPCI1)*BCI1,J1>«UU(J1) 

3400  NEXT  Jl 

3410  NEXT  II 

3420  IF  CASE<6  THEN  3490 

3430  ET=0 

3440  FOR  Jl=l  TO  3 

3450  ET=DST(J1)«EP<J1)+ET 

3460  NEXT  Jl 

3470  Z  =  ET 

3480  GOTO  3590 

3490  IF  CASEO  THEN  3500  ELSE  3570 

3500  HE1=(EP(1)*EP<2) )/2 

3510  HE2=< (EP<1)-EP(2) >/2) "2 

3520  HE3=<HE2*EP<3)~2>~.5 

3530  IF  CASE=1  THEN  3540  ELSE  3550 

3540  MGT=HE1*H£3:Z=MGT:G0T0  3590 

3550  MGT=HE1-HE3:Z=MGT 

3560  GOTO  3590 

3570  MGT=EP(CASE-2> 

3580  Z=MGT 

3590  XP=5»X*Y»CS1 

3600  YP=10*Y»SN1*500000!»Z 

3610  IF  J>1  THEN  3640 

3620  LINE  (1.3»XP-55.180-YP».65) - < 1 . 3»XP-55 . 180-YP" .65) 

3630  GOTO  3650 

3640  LINE  -<1.3»XP-55,180-YP».65> 

3650  RETURN 
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1  '  REM 

20  REM  PROGRAM  NAME:   TRsJACTORIES 

30  REM  *»*»»**#*•**»»**■»■•.*.,,»♦*..» 

40  REM  This   program  gjve3   the  graph   of  tra lector ;e3  of 

stresses . 
50  DIM  X<71>  . YC71)  .NDC52.41  ,W#<3)  ,H#<3>  ,NN<4)  ,FX<4>  ,?X<4>  , 

4) , PY<4) ,BC3,24> ,E(3) 
SO  DIM  0(3,3) ,DB(3,24) ,U(426) .UUC24) .UTC426) ,3IG<3) .ELMTC 
70  EC=2468410!  :REM   modulus  of   elasticity   of   cone.--  I 

(N/sq.cm  ) 
80  M=.2:REM  poisson's  ratio 
90  TH1=60 
100  TH2=10 
110  TH3=55 

120  REM  x,y  coordinate  of  each  node 
130  NM=71:REM  total  node  numbers 
140  NE=52:REM  no  of  total  elements 
150  FOR  1=1  TO  NM 
160  READ  X(I) 
170  NEXT  I 
180  DATA      0.       0.   32.5.   32.5.   40 . 0 .  40 . 0 .  40 . 0 .  4C 

40.0,  40.0,  40.0 
130  DATA   50.0,  50.0,  50.0,  50.0,  50.0,  50.0,  50.0 
200  DATA   60.0,  60.0,  60.0,  60.0,  60.0,  60.0,  60.0 
210  DATA  120.0,120.0,120.0,120.0,120.0,120.0,120.0 
220  DATA  190.0,190.0,190.0,190.0,190.0,190.0,190.0 
230  DATA  260.0,260.0,260.0,260.0,260.0,260.0,260.0 
240  DATA  320.0,320.0,320.0,320.0,320.0,320.0,320.0 
250  DATA  330.0,330.0,330.0,330.0,330.0,330.0,330.0 
260  DATA  340.0,340.0,340.0,340.0,340.0,340.0,340.0 
270  DATA  347.5,347.5,380.0,380.0 
280  FOR  1=1  TO  NM 
290  READ  Yd) 


300 

NEXT 

I 

310 

DATA 

0 

.0 

.43 

.5,: 

L93 

.5,233 

.5 

320 

DATA 

0 

.0, 

43 

.5, 

53 

.5, 

63 

,5, 

,118 

,3, 

,  193 

,  5 

.233 

,5 

330 

DATA 

0 

.0, 

43 

.5, 

53 

.5, 

S3 

.5, 

,118 

.5, 

,193 

.5 

■  233 

.5 

340 

DATA 

0 

•  0, 

43 

.5, 

53 

■  5, 

S3 

.5, 

,118 

.5, 

,193 

.5 

,233, 

,5 

350 

DATA 

0 

.0, 

43 

.5. 

S3 

.5. 

S3 

.5, 

118 

.5, 

,193 

.5 

,233, 

,3 

360 

DATA 

0, 

,0, 

4  3 

.5, 

53 

.5, 

63 

.5, 

118 

.5, 

,193, 

,5 

,233, 

.5 

370 

DATA 

0, 

•  0, 

43 

.5, 

53 

.5, 

63 

.5, 

113. 

.5, 

,193 

,5 

,233, 

,  5 

380 

DATA 

0, 

.0, 

43 

.5. 

53 

,5, 

63 

,5, 

118, 

5, 

193, 

.5, 

,233, 

,5 

390 

DATA 

0, 

.0, 

43 

.5, 

53 

.5. 

63, 

,5, 

118, 

.5, 

193, 

,5, 

,233. 

,3 

400 

DATA 

0 

•0, 

43, 

,5. 

33. 

.5, 

63, 

.5, 

118, 

,3, 

193, 

,5, 

,233. 

,3 

410 

DATA 

193, 

.5. 

233, 

.5, 

0, 

0, 

4  3 

,5 

420  REM 

430  REM  NODE  NO.  CORRESPONDING  TC  EACH  ELEMENT 

440  FOR  EL=1  TO  NE 

450  FOR  EN=1  TO  4 

460  READ  ND(EL,EN) 

470  NEXT  EN:NEXT  EL 

430  DATA   1.  5,  6.  2,   5.12.13.  6.  12.19,20,13, 
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26,33,34,27 
490  DATA  33,40.41.34 

61,70,71,62 
500  DATA   6,13.14.  7 

34,41,42,35 
510  DATA  41.4.3,49.42 

14,21,22,15 
520  DATA  21,28,29.22 

49,56,57,50 
530  DATA  56,63.64,57 

29,36,37,30 
540  DATA  36.43,44,37 

9,16,17,10 
550  DATA  16.23.24,17 

44,51,52,45 
560  DATA  51.58,59,52 

17,24,25,18" 
570  DATA  24,31,32,25 

52,59,60,53 
580  DATA  59,66,67,60,  66,68,69 
590  D1  =  ECV<1-M*M) 
600  D<1 ,1) =D1 :D(1 ,2>=D1*M 
610  D(2,1)=DC1 ,2) :D<2,2) =D1 
620  D(3,3>= (1-M) *Dl/2 
630  OPEN  "R",#5, "B:RRDISP" .4 
640  FIELD  #5,4  AS  DS 
650  FOR  1=1  TO  426 
660  GET  #5,1 
670  U(I>=CVS<DS) 
680  NEXT  I 
690  CLOSE  #5 
700  REM 

710  PI=3.1415927# 
720  CS=C0S(PI/4) 
730  3N=SIN(PI/4J 
740  INFUT"NO.  OF  LINES  " ; NLS 
750  INPUT"N0.  OF  LINES  IN  HORIZ 
760   INPUT"HOW   MANY   ELEMENTS 

GRAPHCtype  no.)":HM 
770  FOR  1=1  TO  HM 

780  INPUT"ELEMENT  NO.Ctyoe  no.) 
790  NEXT  I 
800  FOR  EL=1  TO  HM 

810  IF  ELMT<EL.)=30  OR  ELMT(EL)  = 
820  IF  EL=1  THEN  SIGN=-1  ELSE  S 
830  IF  EL=1  THEN  Z  =  2  ELSE  Z=l 
840  IF  EL=1  THEN  SN=1  ELSE  SN= 
850  FOR  K=l  TO  4 
860  NN(K)=ND<ELMT(EL> ,K> 
870  NEXT  K 

880  5P=Y<NN<3> > -YCNNC2) > 
890  DPTH=XCNN(2) >-XCNN(l) > 
900  T=0 
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910  FOR  11=1  TO  4 

920  FOR  JJ=1  TO  6 

930  T»T*1 

940  IJ»(NN<II)-1)»S*JJ 

9S0  UUCT>=U(IJ) 

960  NEXT  JJ.II 

970  IF  EL>1  THEN  990 

980  CLS 

990  FOR  1=1  TO  4 

1000  J=I*1:IF  J=5  THEN  J=l 

1010  XP=20f XCNNCI)  >  :  YP=YCNNCI)  )  »PP 

1020  XQ=20*XCNNC  J>  )  :  YQ  =  Y  CNN  ( J  )  )  »PP 

1030  LINE  CXP-4-520,200-YP»2) - ( XQ-4-520 , 200-Y0-2 > 

1040  NEXT  I 

1050  DLTY=SP/ CNLS-1) 

1060  DLTX=DPTH/100 

1070  DLTTX=DPTH/NLSS 

1080  FOR  JJ  =  1  TO  2 

1090  Y=0:X=0 

1100  FOR  IY=1  TO  NLS 

1110  Y=YCNN<2> )*CIY-1>»DLTY 

1120  FOR  IX=1  TO  100 

1130  YT  =  Y 

1140  X=XCNNC2> )*(IX-1> "DLTX-SIGN 

1150  XT=X 

1160  X=X*DLTX»SIGN 

1170  GOSUB  1500 

1180  IF  Y>Y(NN(3))  OR  Y<YCNNC2))   THEN  1230 

1190  IF  X>XCNNC2))  OR  X<XCNNC1)>   THEN  1230 

1200  IF  YT>Y<NN(3)>  OR  YT<Y(NN(2J)   THEN  1230 

1210  IF  XT>XCNNC2>>  OR  XT<X<NNC1>)   THEN  1230 

1220  NEXT  IX 

1230  NEXT  IY 

1240  NEXT  JJ 

1250  FOR  JJ=1  TO  2 

1260  Y=0:X=0 

1270  FOR  XX=1  TO  NLS5-1 

1230  IF  EL  =  1  AND  JJ  =  1  THEN  Z  =  3 

1290  IF  EL  =  1  AND  JJ  =  2  THEN  2  =  2 

1300  IF  EL=2  AND  JJ=1  THEN  2=4 

1310  IF  EL=2  AND  JJ=2  THEN  Z=l 

1320  XS=XCNNC2> ) *XX*DLTTX »SIGN : Y= Y CNN C2 ) ) 

1330  FOR  IX=1  TO  100 

1340  YT=Y 

1350  X=XS*CIX-1)»DLTX«5IGN 

1360  XT=X 

1370  X=X*DLTX»SIGN 

1380  GOSUB  1500 

1390  IF  Y>YCNNC3))  OR  Y<YCNNC2))   THEN  1440 

1400  IF  X>X(NNC2>>  OR  X<XCNNC1)>   THEN  1440 

1410  IF  YT>Y(NNC3)>  OR  YT<YCNN(2)>   THEN  1440 

1420  IF  XT>X(NNC2>>  OR  XT<XCNNCD)   THEN  1440 

1430  NEXT  IX 


128 


1440  NEXT  XX 

1450  NEXT  JJ 

1460  NEXT  EL 

1470  LOCATE  24, 22  .-PRINT  "TRAJECTORIES  OF  THE  STRESSES 

1480  AS  =  INKEYSi:IF  AS  =  ""  THEN  1480 

1490  END 

1500  V1=X-XCNN<1> > :L1=X(NN<2> ) -XCNN(l) 1 

1510  V2=Y-Y(NN(D) :L2=Y(NN<4> ) -YCNN(l) ) 

1520  FOR  11=1  TO  2 

1530  V=V1:L=L1 

1540  IF  11=2  THEN  V=V2 

1550  IF  11=2  THEN  L=L2 

1560  GOSUB  2070 

1570  IF  11=2  THEN  1S10 

1580  FXU)=F1 :FX<2)=F2:FX<3)=F3:FX(43=F4 

1590  PX  <  1 )  =F5  :  PX  C  2 >  =F6  :  PX  C  3 >  =F7  :  PX  <  4  >  =F8 

1600  GOTO  1630 

1610  FY(1)=F1 :FY(2)=F2:FY(3)=F3:FY<4)=F4 

1620  PY<1)=F5:PY(2)=F6:PY(3)=F7:PY(4)=F8 

1630  NEXT  II 

1640  FOR  11=1  TO  4 

1650  C=<II-1>«6 

1660  IF  11=1  OR  11=4  THEN  K=l  ELSE  K=2 

1670  IF  II<  =  2  THEN  KK  =  1  ELSE  KK  =  2 

1680  B(l ,1-C) =PX(K) »FY(KK) : B ( 2 , 2-C) =FX (K ) *PY < KK ) 

1690  B(1,3*C)=PX(K*2)«FY<KK) :B (2, 4»C> =FX (K*2> »PY (KK> 

1700  B(1,5*C)=PX(K)*FY(KK»23 : B <2 . 6*C) =FX (K) «PY (KK*2) 

1710  B(3, 1*C> =B<2,2*C) 

1720  B(3,2*C)=B(1.1+C) 

1730  B<3,3*C)=B(2,4*C) 

1740  BC3,4*C>  =  B<1 ,3<-C) 

1750  B(3,5»C)=B(2,6i-C) 

1760  B(3,6*C)=B(1,5*C) 

1770  NEXT  II 

1780  FOR  K=l  TO  3 

1790  E(K)=0 

1800  FOR  KK=1  TO  24 

1810  E(K)=E<K)*B(K,KK)«UU(KK) 

1820  NEXT  KK,K 

1830  FOR  K=l  TO  3 

1840  SIG(K)=0 

1850  FOR  KK=1  TO  3 

1860  SIG<K)=SIG(K)*D(K,KK) »E<KK) 

1870  NEXT  KK,K 

1880  PHI=0 

1890  S1=SIG(1) :S2=SIG(2) :TXY=SIGC3) 

1900  SS=S1-S2 

1910  IF  SS=0  THEN  PHI=0:GOTO  1980 

1920  IF  EL=1  AND  JJ=2  THEN  1960 

1930  IF  EL=1  AND  JJ=1  THEN  1970 

1940  IF  EL=2  AND  JJ=2  THEN  1960 

1950  IF  EL=2  AND  JJ=1  THEN  1970 

1960  PHI=ATN<-2»TXY/(S1-S2) )/2:G0T0  1980 
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1970    PHI=<ATNC-2«TXY/(S1-S2) ) *PI)/2 

1980    TN=TAN(PHI> 

1990    Y=YT*TN»DLTX»SN 

2000  IF  Y>Y(NNC3>>  OR  Y<YCNN(2)>   THEN  2060 

2010  IF  X>XCNN(2>>  OR  X<X(NN<D)   THEN  2060 

2020  IF  YT>YCNN(3)  :>  OR  YT<Y(NN<2>>   THEN  2060 

2030  IF  XT>XCNNC2)>  OR  XT<X<NN(1)>   THEN  2060 

2040  XPT  =  XT*20:YPT  =  YT->-PP:YP  =  Y<-PP:XP  =  X*20 

2050  LINE  (XPT»4-520,200-YPT'2>-(XP*4-520,200-YP*2> 

2060  RETURN 

2070  LV=L-V 

2080  F1=LV/L*V«LV« (L-2»V)/L"3 

2090  F2=V/L*V«LV» <2»V-L) /L~3 

2100  F3=V»LV*2/L"2 

2110  F4=-LV»V~2/L~2 

2120  F5=(-6«V»L+6«V-2)/L-3 

2130  F6=(6«V»L-6»V*2)/L-3 

2140  F7= (L-2-4«L»V*3»V~2)/L'2 

2150  F8= C-2»L»V*3»V"2)/L"2 

2160  RETURN 


130 


10  REM 

20  REM         PROGRAM  NAME: ENERGY 

30  REM         .................... 

40  REM  PROGRAM  TO  VERIFY  ZERO  STRAIN  ENERGY  PICTURE  OF  RIGID 
BODY  MOTION 

50  DIM  X5T<21) ,XND<21) ,Y5T(21) ,YND(21) ,X<4) ,Y<4) ,D<3.3) ,EC3> , 
GC3) ,U<24> ,B<3,24) , N<2,24) , FXC4) , FY<4) ,PX(4) ,PY<4) 

60  FOR  1=1  TO  4 

70  READ  X(I),V(I) 

80  NEXT  I 

90  DATA  0.0,160,0.160.160.0.160 

100  PI=3.1415927# 

110  CS=C0S(PI/4) 

120  SN=SIN(PI/4) 

130  PRINT 

140  INPUT-WHAT  KIND  OF  STRAIN  CONDITION  APPLIED  C   CUi»*<-).l 

&    Vi»*<->.1) type    1    <THETA= .00125) type    2 

CUi«*<-).l  S,  Vl  =  *(-).l  a  THETA=. 00125) type  3]";A 

150  PRINT 

160  INPUT"DO  YOU  WANT  U   or  V";AS 

170  IF  AS="N"  OR  AS="n"  THEN  210 

180  PRINT 

190   INPUT"WHAT   DO  YOU  WANT   U type  1  or   V type  2 

";US 

200  GOTO  270 

210  PRINT  :INPUT"DO  YOU  WANT  EACH  OF  STRAIN       Y/N";A3 

220  IF  AS="N"  OR  AS="n"  THEN  250 

230  PRINT  :INPUT"WHAT  DO  YOU  WANT (Epsilon  x(l)  Epsilon 
y(2)      Gamma  xy(3)      TYPE  N0.";AA 

240  GOTO  270 

250  PRINT  : INPUT"DO  YOU  WANT  THE  5TRAIN  ENERGY  PICTURE 
TYPE  Y  ";AS 

260  SE  =  1 

270  IF  A=2  THEN  360 

280  FOR  1=1  TO  24 

290  READ  U(I) 

300  NEXT  I 

310  DATA  .1,-. 1,0, 0,0,0 

320  DATA  .1, .1,0,0,0,0 

330  DATA  -.1,-1,0,0,0.0 

340  DATA  -.1,-. 1,0, 0,0,0 

350  IF  A=l  THEN  440 

360  THTA=. 00125 

370  FOR  1=1  TO  4 

380  P=6*(I-1>*1 

390  U(P+2)=0 

400  U<P*3)=THTA 

410  U(P*4)=-THTA 

420  U(P+5)=0 

430  NEXT  I 

440  E=60000! :M= . 17 

450  D1=E/(1-M»M) 
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460  D(1,1)=D1:D(1,2>=D1»M 

470  DC2,1)=DC1,2> :DC2,2>=D1:D<3.35=<1-M> »Dl/2 

480  DX=<X<2> -X(l> ) /20 

490  DY=(Y(4>-Y(1) )/20 

500  FOR  1=1  TO  21 

510  XST(I5=X(1> 

520  XNDCI>=XC2> 

530  YST<I)=Y(3) 

540  YNDCI)=Y<2) 

550  NEXT  I 

560  H=0:V=0 

570  CLS 

580  FOR  1=1  TO  4 

590  J»I*1:IF  J=5  THEN  J=l 

600  XP=20*X<I)*Y<I)»CS:YP=25+Y(I)«SN 

610  XQ  =  20  +  X(J)i-Y(J)«CS:YQ  =  25*Y(J)»SN 

620  LINE  <XP«2,220-YP)-(XQo2,220-YQ> 

630  NEXT  I 

640  FOR  1=1  TO  4 

650  X1=XCI) :Y1=Y(I) 

660  12=1*1 

670  IF  12=5  THEN  12=1 

680  X2=X<I2) :Y2=Y(I2) 

690  DPX=<X2-X1)/100:DPY=CY2-Y1)/100 

700  FOR  J=l  TO  101 

710  X=X1*(J-1)»DPX:Y=Y1»(J-1) «DPY 

720  GOSUB  1280 

730  NEXT  J 

740  NEXT  I 

750  FOR  1=1  TO  2 

760  FOR  XI=1  TO  21 

770  IF  1=2  THEN  820 

780  XO=XST<XI> 

790  XEND=XND(XI) 

800  DDX=<XEND-X0)/100 

810  Y=YMN*DY«(XI-1) : GOTO  860 

820  YO=YST<XI> 

830  YEND=YND(XI) 

840  DDY=<YEND-Y0)/100 

850  X=XMN*DX»(XI-1) 

860  FOR  J=l  TO  101 

870  IF  1=2  THEN  890 

880  X=X0+(J-1)»DDX:G0T0  900 

890  Y=Y0+(J-l;«DDY 

900  GOSUB  1280 

910  NEXT  J 

920  NEXT  XI 

930  NEXT  I 

940  IF  SE=1  THEN  1180 

950  IF  US=1  THEN  970 

960  IF  U5  =  2  THEN  990  ELSE  1010 

970  LOCATE  24,15:PRINT  "U" 

980  GOTO  1090 
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990  LOCATE  24.15: PRINT  "V" 

1000  GOTO  1090 

1010  IF  AA  =  1  THEN  1040 

1020  IF  AA  =  2  THEN  1060 

1030  IF  AA  =  3  THEN  1080 

1040  LOCATE  24,5:PRINT  "EDsilon  x" 

1050  GOTO  1090 

1060  LOCATE  24.5:PRINT  "EDsilon  y" 

1070  GOTO  1090 

1080  LOCATE  24,5:PRINT  "GAMMA  xy" 

1090  IF  A  =  l  THEN  1120 

1100  IF  A  =  2  THEN  1140 

1110  IF  A  =  3  THEN  1160 

1120   LOCATE  24,20:PRINT  "DUE  TO  Ul,2=.l  U3.4=-.l  &   VI. 4=-.: 

V2,3=.l" 
1130  GOTO  1260 

1140  LOCATE  24,20:PRINT  "DUE  TO  THETA= . 00125" 
1150  GOTO  1260 
1160   LOCATE   24,20:PRINT  "DUE  TO   Ul,2=.l   U3, 4  =  - . 1 . VI , 4= - . : 

V2,3=.l  £,  THETA=. 00125" 
1170  GOTO  1260 
1180  IF  A=l  THEN  1210 
1190  IF  A  =  2  THEN  1230 
1200  IF  A  =  3  THEN  1250 
1210  LOCATE  24.5  : PRINT"STRAIN  ENERGY  DUE  TO  U1.2=.l  U3.4=-.: 

S.  VI, 4  =  -. 1  V2,3=.l" 
1220  GOTO  1260 

1230  LOCATE  24.5  : PRINT"5TRAIN  ENERGY  DUE  TO  THETA= . 00125" 
1240  GOTO  1260 
1250   LOCATE  24.5  :PRINT"5TRAIN  ENERGY  DUE  TO  U1.2».l   U?.4=- 

.1,  VI, 4=-. 1  V2,3=.l  ,  and  THETA= . 00125" 
1260  AS=INKEYS:IF  AS  =  ""  THEN  1260 
1270  END 

1280  V1=X-X(1) :L1=XC2)-X<1) 
1290  V2=Y-Y<1) :L2=Y<4)-Y(1) 
1300  FOR  11=1  TO  2 
1310  V=V1:L=L1 
1320  IF  11=2  THEN  V  =  V2 
1330  IF  11=2  THEN  L=L2 
1340  LV=L-V 

1350  F1=LV/L*V»LV»(L-2«V)/L~3 
1360  F2=VyL*V»LV»<2»V-L)/L-3 
1370  F3=V»LV"2/L*2 
1380  F4=-LV»V~2/L"2 
1390  IF  US  =  1  OR  US  =  2  THEN  1440 
1400  F5=  C-6»V«L  +  6»V2)/L~3 
1410  F6= (6»V»L-6»V~2)/L"3 
1420  F7=(L-2-4»L»V"-3«V~2)/L~2 
1430  F8=(-2«L»V+3»V-2)/L~2 
1440  IF  11=2  THEN  1480 

1450  FX(1)=F1 :FX(2)=F2:FX(3)=F3:FX(4) =F4 
1460  PX(1)=F5:PX(2)=F6:PX(3>=F7:PX(4)=F8 
1470  GOTO  1500 
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1480  FYC1:>=F1:FYC2>=F2:FY<3'>=F3:FY<4)=F4 

1490  PY(1)=F5:PY(2)=F6:PY(3)=F7:PY(4)=F8 

1500  NEXT  II 

1510  FOR  11  =  1  TO  4 

1520  IF  11=1  THEN  C=0 

1530  IF  11=2  THEN  C  =  6 

1540  IF  11=3  THEN  C=12 

1550  IF  11=4  THEN  C=18 

1560  IF  11=1  OR  11=4  THEN  K=l  ELSE  K=2 

1570  IF  II<=2  THEN  KK=1  ELSE  KK=2 

1580  N(1,1+C)=FX(K)«FY(KK) 

1590  N(1,3*C>=FX(K*2)»FY(KK) 

1600  N(1,5*C)=FX(K> »FY(KK+2) 

1610  N(2,2*C)=N(1,1-C) 

1620  N<2.4*C)=N<1,3*C) 

1630  N(2,6*C)=N(1,5*C) 

1640  IF  US  =  1  OR  US=2  THEN  1710 

1650  B(1,1*C)=PX<K)«FY<KK) : B<2,2*C> =FX CK) »PY CKK> 

1660  B(1,3*C)=PX(K*2)»FY(KK) :B (2, 4+C> =FX <K  +  2) »PY (KX) 

1670  B(1,5*C>=PX(K)»FYCKK*2) : B <2, 6*C) =FX <K) »PY (KK+2) 

1680  B(3,1*C)=B(2,2*C) :B (3,2*C) =B < 1 » 1 +C) 

1690  B<3,3*C>=B<2,4+C> : B < 3,4+C> =B ( 1 , 3+C) 

1700  B(3,5-C)=B(2,6*C) :B(3,6»C)=B(1,5*C> 

1710  NEXT  II 

1720  IF  US  =  1  OR  US  =  2  THEN  1940 

1730  FOR  K=l  TO  3 

1740  E(K)=0 

1750  FOR  KK=1  TO  24 

1760  E(K)=E(K)-»B(K,KK)  »U(KK) 

1770  NEXT  KK,K 

1780  IF  SE=0  THEN  1900 

1790  FOR  K=l  TO  3 

1800  G(K)=0 

1810  FOR  KK=1  TO  3 

1820  G(K)=G(K)*D(K.KK)»E(KK> 

1830  NEXT  KK,K 

1840  FOR  K=l  TO  3 

1850  GG  =  GG-i-E<K)«G(K) 

1860  NEXT  K 

1870  GG=GG/2 

1880  2=GG:GG=0 

1890  S2=800:G0T0  2010 

1900  IF  AA=1  THEN  Z=E(1) 

1910  IF  AA=2  THEN  2=E<2) 

1920  IF  AA=3  THEN  2=E<3> 

1930  SZ=10000:G0T0  2010 

1940  UU=0 

1950  FOR  11=1  TO  24 

1960  IF  US=1  THEN  G=l  ELSE  G=2 

1970  UU=UU*U<II)«N(G,II) 

1980  NEXT  II 

1990  2=UU:S2=150 

2000  IF  A=2  THEN  SZ=500 
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2010  XP=20*X*Y*CS 

2020  YP  =  25h-Y»SN*SZ»2 

2030  IF  J>1  THEN  2050 

2040  LINE  <2«XP,220-YP>-<2«XP,220-YP) :G0T0  2060 

2050  LINE-<2»XP,220-YP> 

2060  RETURN 
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10  REM 

20  REM  PROGRAM  NAME:  HERMITIAN 

30  REM  **»***»#»* w * w »»#»*»»#» • 

40  DIM  XSTC21) ,XND(21) ,YSTC21) , YND(21> ,X(4) , Y(4) 

50  FOR  1=1  TO  4 

60  READ  X(I) ,Y(I) 

70  NEXT  I 

80  DATA  0,0,120,32,168,100,84.128 

90  PI=3. 1415927* 

100  CS=C0S(PI/4) 

110  SN=SIN<PI/4) 

120  INPUT"NODE  OF  INTEREST" ; ISN 

130   INPUT"DO   YOU  WANT  THE  GRAPH  WHEN  UNIT   NODAL   PARAMETER 

APPLIES  Y/N'^AS 
140  IF  AS="Y"  OR  AS="y"  THEN  270 

150  INPUT"D0  YOU  WANT  2nd  ORDER  £  Cx.y)  Y/N";AS5 
160  IF  AS="N"  OR  AS="n"  THEN  190 
170   INPUT"WHAT   DO  YOU  WANT  <W=1  W,x  =  2  W,y  =  3   W,xx=4   W.xy=5 

W,yy=6)  TYPE  N0.";F 
180  GOTO  290 

190  INPUT"D0  YOU  WANT  fCx.y)  Y/N";AS 
200  IF  ASs  =  "N"  OR  As3="n"  THEN  220 
210  F=7:G0T0  290 

220  INPUT"DO  YOU  WANT  df/dx  Y/N";AS 
230  IF  A«="N"  OR  AS="n"  THEN  250 
240  F=8:G0T0  290 

250  INPUT"DO  YOU  WANT  df/dy  Y/N";AS 
260  IF  AS="Y"  OR  AS="y"  THEN  F=9:G0T0  290 
270   INPUT"DO  YOU  WANT  THE  GRAPH  BY  (lst(l)  OR  2nd  <2>  ORDER) 

TYPE  NO.";0D 
280   INPUT"WHAT  GRAPH  DO  YOU  WANT(u=l  v=2    u.x  =  3   v,x  =  4   u.y  =  5 

v,y=6)  TYPE  NO . " ; S 
290  XM=l:YM=l:XN=l:YN=l:XMX=X(l) :YMX=YC1) :XMN=XC1> :YMN=YCI) 
300  FOR  K=l  TO  4 
310  IF  XMX>X(K)  THEN  330 
320  XMX=X(K) :XM=K 
330  IF  YMX>Y(K)  THEN  350 
340  YMX=Y<K> : YM=K 
350  IF  XMN<X(K)  THEN  370 
360  XMN=X<K) :XN  =  K 
370  IF  YMN<Y(K)  THEN  390 
380  YMN=Y(K) :YN=K 
390  NEXT  K 
400  DX=<XMX-XMN)/20 
410  DY=(YMX-YMN)/20 

420  Nl=YN:N2=YN*l:IF  N2=5  THEN  N2=l 
430  N3=N1*2:IF  N3=5  THEN  N3=l 
440  N4=N1-1:IF  N4=0  THEN  N4  =  4 
450  NNl=XN:NN2=XN+l:IF  NN2=5  THEN  NN2=1 
460  NN3=NNl+2: IF  NN3  =  5  THEN  NN3=1 
470  NN4=NNl-l:IF  NN4=0  THEN  NN4=4 
480  FOR  1=1  TO  21 
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490  Y  =  YMN*(I-1.1  »DY 

500  IF  Y<=Y<N4>  THEN  540 

510  K=CY-Y<N4> )/(Y(N3)-Y(N4> > 

520  XST(I)=X(N4) *K» <X <N3> -X <N4> ) 

530  GOTO  560 

540  K- (Y-Y(Nl) )/ (Y(N4>-YCN1) > 

550  XST(I)=X<N1)+K»<XCN4) -XCN1) ) 

560  IF  Y<=Y(N2>  THEN  620 

570  IF  YM  =  N4  AND  Y>Y(N3)  THEN  600 

5S0  K=CY-Y(N2)  ~>  J  <Y(N3)-Y(N2)  > 

590  XND<I)=X<N2)*K»CX(N3)-X(N2> ) :GOTO  640 

600  K=<Y-Y(N3) >/ <Y<N4>-Y<N3) ) 

610  XND<I>=X<N3)<-K*<X(N4)-X(N3>> : GOTO  640 

620  K=<Y-Y(N1) )/(Y(N2) -Y(N1) ) 

630  XNDCI)=XCNl>i-K*<XCN2)  -X(N1)  ) 

640  X=XMN*(I-1) »DX 

650  IF  X<=X<NN4>  THEN  690 

660  K=(X-X(NN4))/(X<NN3)-X<NN4> ) 

670  YST(I>=YCNN4> *K» C Y < NN3) - Y ( NN4 > ) 

680  GOTO  710 

690  K=<X-X<NN1)>/CX(NN4)-X<NN1) ) 

700  YST<I)=YCNN1)*K«(YCNN4) -Y(NN1) ) 

710  IF  X<=XCNN2)  THEN  740 

720  K=(X-X<NN2> ) / (X <NN3> -X ( NN2> ) 

730  YND(I)=Y(NN2)*K»(Y(NN3)-Y(NN2) ) : GOTO  760 

740  K=(X-X(NN1> >/<X<NN2)-X<NNl>) 

750  YND(I)=Y(NNl)->-K»<Y(NN2)  -Y(NN1)  ) 

760  NEXT  I 

770  H=0:V=0 

780  CLS 

790  FOR  1=1  TO  4 

800  J=I*1:IF  J=5  THEN  J  =  l 

810  XT=X(I) :YT=Y<I) 

S20  XP=20+X<IJ*Y(I)»CS:YP=25fY(I)»SN 

830  XT=XCJ) :YT=Y(J) 

840  XQ=20*X( J>  *Y( J>  *CS: YQ  =  25*Y(J) «SN 

850  IF  H<XP  THEN  H=XP 

860  IF  H<XO  THEN  H=XQ 

870  IF  V<YP  THEN  V  =  YP 

880  IF  V<YQ  THEN  V=YQ 

890  LINE  <XP»2,220-YP)-(XQ»2,220-YQ) 

900  NEXT  I 

910  FOR  1=1  TO  4 

920  X1=X<I) :Y1=Y(I) 

930  12=1*1 

940  IF  12=5  THEN  12=1 

950  X2=X(I2) :Y2=Y(I2) 

960  DPX=  <X2-X1)/100:DPY=  (Y2-YD/100 

970  FOR  J=l  TO  101 

980  X=X1*(J-1)*DPX:Y=Y1*(J-1) *DPY 

990  GOSUB  1740 

1000  NEXT  J 

1010  NEXT  I 
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1020  FOR  1=1  TO  2 

1030  FOR  XI=1  TO  21 

1040  IF  1=2  THEN  1090 

1050  XO=XST(XI) 

1060  XEND=XNDCXI) 

1070  DDX=(XEND-X0>/100 

1080  Y=YMN*DY»(XI-1) : GOTO  1130 

1090  YO=YST(XI) 

1100  YEND=YND(XI> 

1110  DDY=(YEND-Y0)/100 

1120  X=XMN+DX»(XI-1> 

1130  FOR  J=l  TO  101 

1140  IF  1=2  THEN  1160 

1150  X  =  XO-i-<J-l)*DDX:GOT0  1170 

1160  Y  =  Y0<-CJ-1)»DDY 

1170  GOSUB  1740 

1180  NEXT  J 

1190  NEXT  XI 

1200  NEXT  I 

1210  IF  F=>1  GOTO  1340 

1220  IF  S>=1  AND  0D=1  THEN  1240 

1230  IF  S>=1  AND  0D=2  THEN  1260 

1240   LOCATE   2,20:PRINT  "By  using  1st  order  Hermitian   Shaoe 

Function" 
1250  GOTO  1280 
1260   LOCATE  2.20:PRINT  "By  using  2nd  order   Hermitian   Shape 

Function" 
1270  GOTO  1280 
1280  IF  S=l  THEN  1430 
1290  IF  S=2  THEN  1450 
1300  IF  S=3  THEN  1470 
1310  IF  S=4  THEN  1490 
1320  IF  S=5  THEN  1510 
1330  IF  S=6  THEN  1530 
1340  IF  F=l  THEN  1550 
1350  IF  F=2  THEN  1570 
1360  IF  F=3  THEN  1590 
1370  IF  F=4  THEN  1610 
1380  IF  F  =  5  THEN  1630 
1390  IF  F=6  THEN  1650 
1400  IF  F=7  THEN  1670 
1410  IF  F=8  THEN  1690 
1420  IF  F=9  THEN  1710 
1430   LOCATE  24 ,20 : PRINT"Epailon  x  due  to  u=l  &  Gammma  xy  due 

to  v=l  at  node";ISN 
1440  GOTO  1720 
1450  LOCATE  24 ,20  :PRINT"Epsilon  y  due  to  v  =  l  &  Gammma  xy   due 

to  u=l  at  node";ISN 
1460  GOTO  1720 
1470   LOCATE   24 ,  20  :PRINT"Epsilon  x  due  to  u,x  =  l  &  Gammma   xy 

due  to  v,x=l  at  node";ISN 
1480  GOTO  1720 
1490   LOCATE  24,20:PRINT  "Epsilon  y  due  to  v,x  =  l  6.  Gammma  r.y 
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due  to  u,x  =  l  at  node" : I3N 
1500  GOTO  1720 
1510   LOCATE  24,20:PRINT  "Epsilon  x  due  to  u,y  =  l  &  Ganmma   xy 

due  to  v,y=l  at  node" : ISN 
1520  GOTO  1720 
1530   LOCATE  24.20:  PRINT  "Epsilon  y  due  to  v,y=l  &  Gammmn  xy 

due  to  u,y=l  at  node";ISN 
1540  GOTO  1720 
1550   LOCATE  24.5:PRINT  "2nd  Order  Haniiion   Shape   Function 

using  nodal  w  of  node"; ISN 
1560  GOTO  1720 
1570   LOCATE   24,5:PRINT  "2nd  Order  Hermitian  Shape   Function 

using  nodal  dw/dx  of  node" ; ISN 
1580  GOTO  1720 
1590   LOCATE  24,5:PRINT  "2nd  Order  Hermitian   Shape   Function 

using  nodal  dw/dy  of  node" ; ISN 
1600  GOTO  1720 
1610   LOCATE   24,5:PRINT  "2nd  Order  Hermitian  Shape   Function 

using  nodal  d<dw/dx)/dx  of  node" ; ISN 
1620  GOTO  1720 
1630   LOCATE  24,5:PRINT  "2nd  Order  Hermitian   Shape   Function 

using  nodal  d<dw/dx)/dy  of  node" ; ISN 
1640  GOTO  1720 
1650   LOCATE   24,5:PRINT  "2nd  Order  Hermitian  Shape   Function 

using  nodal  d(dw/dy)/dy  of  node";ISN 
1660  GOTO  1720 
1670   LOCATE  24,5:PRINT  "1st  Order  Hermitian   Shape   Function 

using  nodal  w  of  node"; ISN 
1680  GOTO  1720 
1690   LOCATE   24,5:PRINT  "1st  Order  Hermitian  Shape   Function 

using  nodal  dw/dx  of  node"; ISN 
1700  GOTO  1720 
1710   LOCATE  24, 5: PRINT  "1st  Order  Hermitian   Shape   Function 

using  nodal  dw/dy  of  node" ; ISN 
1720  AS=INKEYS  :  IF  AS=""  THEN  1720 
1730  END 

1740  K2=ISN*l:IF  K2>4  THEN  K2=K2-4 
1750  K3=ISN*2:IF  K3>4  THEN  K3=K3-4 
1760  K4=ISN-3:IF  K4>4  THEN  K4=K4-4 
1770  L1=<XCK2)-X)»(Y(K3>-Y>-<X<K3)-X)«<Y(K2)-Y> 
1780  L2=(X(K3)-X)«(Y<K4)-Y)-<X<K4)-X) »(Y(K3)-Y) 
1790   L3=<X<K2)-X<ISN> )*< Y<K3)-Y<ISN) >-CX<K3> -XCISN>  >"<Y(K2>- 

Y(ISN)) 
1800   L4=<X(K3)-X(ISN) ) »  <  Y  <K4  >  -  Y  <  ISN  > ) - CX < K4> -X < ISN ) ) * < Y < K3 ) - 

YCISN)) 
1810  L34=L3»L4 
1820  L=L1»L2/L34 
1830  XX=X-X(ISN) :YY=Y-Y(ISN> 
1840  0ML=1-L 

1850  IF  0D=2  OR  (F>=1  AND  F<7)  THEN  1880 
1860  H=1-3»0ML'2*2»0ML*3 
1870  PHL=6»0ML-6»0ML"2  : GOTO  1900 
1880  H=l-10-0ML-3*15«0ML-4-6»0ML-S 
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1S90  PHL=30«0ML-2-60»0ML~3*30«0ML' 

1900  IF  F>  =  7  THEN  2340 

1910  IF  F>»1  AND  F<7  THEN  2490 

1920  LX1=Y(K2)-Y<K3) 

1930  LY1=X(K3)-X(K2) 

1940  LX2=Y<K3>-Y<K4> 

1950  LY2=X(K4)-X<K3) 

1960  LX=CL1»LX2*L2*LX1)/L34 

1970  LY=(L1»LY2*L2»LY1)/L34 

1980  PHX=PHL»LX 

1990  PHY=PHL-LY 

2000  IF  5=2  THEN  2190 

2010  IF  S=3  THEN  2100 

2020  IF  S=4  THEN  2230 

2030  IF  S=5  THEN  2140 

2040  IF  S=6  THEN  2280 
2050  Z=PHX:S2=1000 

2060  IF  ISN=2  THEN  SZ=4000 

2070  IF  ISN=3  THEN  SZ=3000 

2080  IF  ISN=4  THEN  SZ=2000 

2090  GOTO  2780 

2100  2=H*XX»PHX:SZ=30 

2110  IF  ISN  =  2  THEN  S2=10 

2120  IF  ISN>=3  THEN  SZ  =  45 

2130  GOTO  2780 

2140  Z=YY«PHX:SZ=60 

2150  IF  ISN=2  THEN  SZ=30 

2160  IF  ISN=3  THEN  SZ=135 

2170  IF  ISN=4  THEN  SZ  =  90 

2180  GOTO  2780 

2190  Z=PHY:SZ=8000 

2200  IF  ISN=2  THEN  SZ=800 

2210  IF  ISN>=3  THEN  SZ=4000 

2220  GOTO  2780 

2230  Z=XX»PHY:SZ=180 

2240  IF  ISN  =  2  THEN  SZ=60 

2250  IF  ISN=3  THEN  S2=45 

2260  IF  ISN  =  4  THEN  SZ=30 

2270  GOTO  2780 

2280  Z=H*YY»PHY:SZ=100 

2290  IF  ISN=2  THEN  SZ=60 

2300  IF  ISN  =  3  THEN  SZ=35 

2310  IF  ISN=4  THEN  SZ=25 

2320  GOTO  2780 

2330  IF  F<=6  THEN  2490 

2340  IF  F  =  7  THEN  2370 

2350  IF  F=8  THEN  2400 

2360  IF  F  =  9  THEN  2450 

2370  Z=H:SZ=150 

23S0  IF  ISN>=3  THEN  SZ  =  75 

2390  GOTO  2780 

2400  Z=H»XX:SZ=6 

2410  IF  ISN=2  THEN  SZ=.24 
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2420  IF  ISN=3  THEN  SZ=.8 

2430  IF  ISN  =  4  THEN  SZ  =  2 

2440  GOTO  27S0 

2450  Z=H*YY:SZ=3 

2460  IF  I5N  =  2   THEN  5Z  =  1.1 

2470  IF  ISN  =  4  THEN  SZ=1.8 

2480  GOTO  2780 

2490  IF  F"l  THEN  2550 

2500  IF  F  =  2  THEN  2600 

2510  IF  F=3  THEN  2650 

2520  IF  F  =  4  THEN  2690 

2530  IF  F=5  THEN  2730 

2540  IF  F=6  THEN  2760 

2550  Z=H 

2560  IF  ISN  =  1  OR  I5N  =  2  THEN  SZ=150 

2570  IF  ISN=3  THEN  SZ=75 

2580  IF  ISN  =  4  THEN  SZ  =  50 

2590  GOTO  2780 

2600  Z=H»XX:SZ=4 

2610  IF  ISN=2  THEN  SZ=.4 

2620  IF  I5N=3  THEN  SZ=2 

2630  IF  ISN=4  THEN  SZ=1.6 

2640  GOTO  2780 

2650  Z=H»YY:SZ=3 

2660  IF  ISN  =  2  THEN  SZ=1.3 

2670  IF  ISN=4  THEN  SZ=1.5 

2680  GOTO  2780 

2690  Z=H«XX~2/2:SZ=.12 

2700  IF  ISN  =  2  OR  ISN=3  THEN  SZ=.07 

2710  IF  ISN  =  4  THEN  SZ=.24 

2720  GOTO  2780 

2730  Z=H«XX»YY:SZ=.05 

2740  IF  ISN=2  OR  ISN=3  THEN  SZ=.l 

2750  GOTO  2780 

2760  Z=H»YY"2/2:SZ= .08 

2770  IF  ISN=2  OR  ISN  =  3  THEN  SZ=.24 

2780  XP=20*X*Y*CS 

2790  YP=25*Y«SN*SZ«Z 

2800  IF  J>1  THEN  2820 

2810  LINE  <2»XP,220-YP)-(2»XP.220-YP> : GOTO  2830 

2820  LINE-(2*XP,220-YP) 

2830  RETURN 
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ABSTRACT 

The  state-of -the-art  of  stress  analysis  and  design  of 
shear  walls  were  reviewed.  The  mathematical  properties  of 
Hernitian  shape  functions  were  investigated. 

The  development  of  a  sequence  of  computer  programs 
based  on  the  use  of  finite  element  method  as  well  as  the  data 
management  and  reduction  for  the  Zenith  100  computer  was 
discussed  and  presented.  The  applications  of  these  programs 
were  illustrated  by  solving  a  model  shear  wall.  The  obtained 
results  were  displayed  by  computer  plots.  The  numerical 
values  were  compared  to  the  experimental  data.  The  findings 
through  the  investigation  were  summarized  in  the  conclusions. 


